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Abstract 



Effective field tfieory tecfiniques are used to describe tfie interaction of fieavy fiadrons 
in a model independent way. Predictability is obtained by exploiting the symmetries 
of QCD. Heavy hadron chiral perturbation theory is reviewed and used to describe 
D* decays. The phenomenologically important D*Dt: coupling is extracted from data 
working to first order in the chiral and heavy quark symmetry breaking parameters. 
A method is described for determining jl^bl from exclusive semileptonic B and D 
decays with 10% uncertainty. An effective field theory for two-nucleon systems is 
then discussed. The large S-wavc scattering lengths necessitate expanding around a 
non-trivial fixed point. A detailed discussion of the interplay between renormalization 
and the power counting is given. In power counting pion interactions with nucleons 
it is useful to consider three classes of pion: potential, radiation, and soft. A power 
counting for massive radiation is developed. Finally, it is shown that the leading terms 
in the effective theory for nucleon-nucleon interactions are invariant under Wigner's 
SU(4) spin-isospin symmetry in the infinite scattering length limit. 
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Chapter 1 Introduction 



The minimal standard model is an appealing theory which describes the strong, weak, 
and electromagnetic interactions in terms of 19 input parameters. This renormaliz- 
able quantum field theory gives quite an accurate description of nature, as shown 



through precision tests of QED p|, the Electro- Weak sector [jTO|, |TT|, and to a lesser 
extent QCD [0, |T3|, [1^. These tests examine observables for which a perturbative 
treatment of the couplings is applicable. For QCD, this is a valid approach for high 
energy processes due to asymptotic freedom p!5| , [1^. However, at low energy or large 
distance the coupling becomes strong and the quarks and gluons are confined into the 
observed mesons and baryons. At these energies a non-perturbative approach, such 
as lattice gauge theory, is necessary. Since lattice calculations are still fairly crude, 
it is reasonable to ask if the non-perturbative nature of QCD can be handled in an- 
other model independent fashion. In certain situations the answer is yes, because the 
symmetries and dynamics of QCD provide other expansion parameters besides the 
strong coupling. Expanding about a symmetry limit provides us with a means for 
describing non-perturbative effects by a series of low energy parameters (matrix ele- 
ments or effective couplings) which can be determined from experimental data. This 
approach is predicative since there are typically several observables that depend on 
a given parameter. In many ways this is complimentary to lattice QCD calculations, 
which can then concentrate on calculating these parameters. 

Effective field theory is a useful tool for implementing these ideas. We begin by 
writing down fields for the relevant degrees of freedom, and constructing an effective 
Lagrangian. The Lagrangian includes all possible terms that transform correctly 
under the symmetries, and is typically non-renormalizable with an infinite number 
of terms. These terms are organized in importance by power counting in a small 
parameter. Identifying a small expansion parameter usually depends on having scales 
which are widely separated. The low energy Lagrangian is a sum of terms of the 
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form L ~ C(/x)0(/x). The coefficients C(/i) and operators encode short and 

long distance physics respectively, and the renormalization point /i separates the 



two regimes. This is essentially a Wilsonian operator product expansion |T^. The 
effective field theory approach is important for several reasons. In an effective field 
theory different scales in the problem are separated, so that one can concentrate on the 
most interesting physics at a particular scale. Furthermore, the power counting gives 
us a way to estimate the uncertainty in working at a given order. Finally, calculations 
are often much simpler in the effective theory. Depending on the situation, effective 
field theories are used in two somewhat distinct ways, either from the toj) down or 
from the bottom up. 

In a top down approach the high energy theory is understood, but we find it 
useful or necessary to use a simpler theory at lower energies. Since the high energy 
theory is known, the C(/i) couplings can be calculated by performing a perturbative 
matching at the high scale. The theory is then run down to the desired low energy 
scale using the renormalization group. Solving the renormalization group equations 
for the running of the coefficients, we sum potentially large logarithms between the 
two scales. At the low scale, matrix elements of the operators are natural in size. 
A standard example is the calculation of QCD corrections to weak processes at mo- 



menta p ^ 90 GeV (see Ref. [|T8| for a review). Here integrating out the W and 
Z leaves four-fermion interactions and an expansion in p^ /rn^. A second example 
is non-relativistic QED (NRQED) |TP[, which is used in describing the electromag- 
netic interactions of non-relativistic leptons. NRQED is especially useful in describing 
Coulombic bound states such as positronium, where a pure coupling constant expan- 
sion is inappropriate. Instead, a dual expansion is performed in the electromagnetic 
coupling and the velocity of the non-relativistic leptons. This effective theory is es- 
pecially tractable since both the coefficients and matrix elements can be calculated. 
In QCD the quark masses are such that Aqcd ^ ''^c,h = ^q- In the limit mq oo, 
QCD exhibits additional flavor and spin symmetries, called heavy quark symmetry 



(HQS) pO| , [21[] . Heavy Quark Effective Theory (HQET) uses these symmetries and 
an expansion in Aqcd/^^q to make predictions for processes involving hadrons con- 
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taining one heavy quark. At high energies this effective field theory is matched onto 
QCD. In this case the matrix elements are typically not calculable, but are still re- 
lated by HQS. For systems with two heavy quarks, the appropriate effective theory 
is called non-relativistic QCD (NRQCD) [||]. 

A second approach to effective field theory is from the bottom up. In this case 
the high energy theory is either unknown or not calculable. A well known example 
is SU(3) chiral perturbation theory, which exploits the pattern of dynamical chiral 
symmetry breaking observed in QCD. In the limit rriq —>■ 0, QCD has additional 
chiral symmetries, giving enhanced predictive power. Approximate chiral symmetry 
is a result of the small light quark masses, rrig = rriu^^s ^ Aqcd- Phenomenologically, 
this approach is valid for energy and momenta ^ A where A ~ 1.2 GeV is the chiral 
symmetry breaking scale. The relevant degrees of freedom here are the pions, kaons, 
and eta which are the pseudo-Goldstone bosons of the SU{?>)l x SU{?>)r SU{?>)v 
breaking. At low energy, matching onto QCD is not possible, so the couplings C(/i) in 
this low energy theory must be determined from experimental data. However, because 
our fields correspond to the asymptotically observed particles, the matrix elements 
are calculable. Processes involving a single heavy hadron can also be incorporated 
in this approach by combining the power counting in HQET and chiral perturbation 
theory into heavy hadron chiral perturbation theory |2^, The effective 

field theory approach has also been extended to processes with two or more heavy 
particles, such as nucleon-nucleon interactions [|2^, |2^, ^ ^, The latter 
theory will be discussed in some detail. 

In the modern view, the standard model itself is a low energy effective theory. 
As an effective field theory it includes the usual Lagrangian as well as operators of 
dimension five and higher built out of standard model fields. Such operators are 
suppressed by powers of a scale A, where A is a measure of the energy at which the 
new physics becomes relevant. At energies ~ A the standard model effective field 
theory must be replaced by something more fundamental. The fact that the standard 
model is renormalizable is significant since it implies that the scale A is not generated 
by standard model interactions, and is therefore, in principal, unconstrained. Large 
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values of A then explain several of the beautiful features of the standard model, such 
as baryon and lepton number conservation and the absence of flavor changing neutral 
currents. 

In this thesis several applications of effective field theory are discussed. The focus 
will be on using chiral perturbation theory for processes with heavy particles. Much 
of this material has now been published [0, ^ 0, ^ |38|, but unpublished 
material appears in sections p.3.2| and chapter 8. Chapter 2 reviews the necessary 
theoretical tools and establishes notation. In chapter 3, D* decays are investigated. 
The D*^ ^ D*'^, and D* branching fractions are used to extract the D*Dtt and D*D'y 
couplings working to first order in the symmetry breaking parameters, rrig and l/rric- 
Important effects due to the heavy meson mass splittings and unknown order niq 
couplings are included. Predictions for the D* and B* widths are given. Chapter 4 
discusses a method for determining \Vub\ from exclusive B semileptonic decay. The 
calculable deviation from unity of the double ratio of form factors {f^^^^^/f^^~*^*^)/ 
{f(D^p)/f(^^^*)) is determined using chiral perturbation theory and is found to be 
small. It is shown that combining experimental data from B ^ p £ B K*£i and 
D ^ piui can lead to a model independent determination of \ Vub\ with an uncertainty 
from theory of about 10%. 

In chapter 5 an effective field theory for nucleon-nucleon interactions is discussed. 
The power counting in this theory is controlled by the presence of a non-trivial ultra- 
violet fixed point or, equivalently, a bound state near threshold. Two renormalization 
schemes which have manifest power counting are discussed in detail, the Power Diver- 
gence Subtraction scheme (PDS)|31[] and an off-shell momentum subtraction scheme 
which we call the OS scheme. Comparing results in these schemes gives us a method 
for determining if a statement about the behavior of the theory is scheme depen- 
dent. The effect of low energy poles on the organization of the perturbation series is 
explained. Comments are also made regarding the constraints that ultraviolet diver- 
gences make on the power counting. Theoretical and empirical arguments are then 
given about the range of this theory. 

In chapter 6, radiative pion effects are discussed. It is shown that for the purpose of 
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power counting the pion interactions should be divided into three classes: potential, 
radiation, and soft. A power counting is introduced for systematically including 
radiation pion effects. The leading order radiation pion graphs for nucleon-nucleon 
scattering are evaluated. The power counting for soft pions is also discussed. 

Chapter 7 discusses the symmetries of the lowest order nucleon effective field 
theory. It is shown that in the limit where the NN ^Sq and ^Si scattering lengths, 
qC^So) g^j^^j a^^^'^\ go to infinity, the leading terms in the effective field theory for 
strong NN interactions are invariant under Wigner's SU(4) spin-isospin symmetry. 
This explains why the leading effects of radiation pions on the S-wave A^'A'" scattering 
amplitudes vanish as a^^'^°) and a^^^^^ go to infinity. Implications of this symmetry 
are also discussed for A"A^ — > A^A^axion and ^ np. 

A brief discussion of predictions for the ^5*1 — ^ Di mixing parameter ei is given 
in chapter 8. Working in the theory with pions at NNLO gives a one parameter 
prediction for ei(p). The accuracy of this prediction is compared to results in the 
theory without pions. 

Chapter 9 contains concluding remarks. 
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Chapter 2 Theoretical Background 

This chapter introduces the chiral perturbation theory formahsm for theories with 
zero, one or two heavy particles. We begin with the QCD Lagrangian, 

Cqcd = -\g%G^^'' + q{ilf) - m,)q + Q{ilp - mQ)Q + g.f. + c.t. , (2.1) 

where G^^, is the field strength for the gluon field A^, = + igA^T^ is the color 
covariant derivative, g.f. stands for gauge-fixing and ghost terms, and c.t. stands 
for counterterms. The field q includes the three light quark fields u, d, s with masses 
"^q = TT^u,d,si while Q includes the three heavy quarks c, 6, t with masses uiq = rric^h^t. 
The quark-gluon interaction is fiavor blind so in QCD only the masses distinguish the 
quarks. The Lagrangian in Eq. is renormalizable, Lorentz invariant, and is also 
invariant under parity, charge conjugation, and time reversal^ The quantum theory 
of QCD depends on another dynamically generated scale, Aqcd ~ 250 MeV, where 

, , qiu)"^ 4:71 
a, /X = ^ = +.... 2.2 

Here /i is the renormalization point, and /5o = llNc/3 — 2nf/3 is the lowest order 
coefficient of the QCD beta function for nj fiavors and Nc{= 3) colors. QCD is 
asymptotically free [|15[ [l^, as{fi oo) 0, making perturbation theory valid at 
large energies. At low energy as{fi) becomes large, and the quarks and gluons become 
confined. Confinement is a non-perturbative phenomenum and a direct proof from 
QCD has not been given. 

In the limit rrig the light quark term in Eq. ( |2.1|) is invariant under the chiral 



^Motivated by instanton configurations, a term of the form 0/(647r^) G^^G^^" can be added to 
'Cqcd- This term violates parity and time-reversal invariance. In nature is tiny, limits on the 
neutron electric dipole moment give 9 < 10~^. The occurrence of this unnaturally small value 
is known as the strong CP problem. 



symmetry transformation 
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qi^LqL, qR^Rqn, where L e SU (3) l and R e SU (3) r. (2.3) 

Since rriq ^ Aqcd this is an approximate symmetry of QCD. This symmetry is 
spontaneously broken, SU{3)l x SU{3)fi SU{3)v, by the vacuum expectation 
value 

{0\q-q'^\0)=v S^' , where v ~ A^qcd • (2-4) 

The breaking of chiral symmetry is another non-perturbative effect, and occurs at a 
scale ~ 1 GeV. The up and down quarks are much lighter than the strange quark, 
so SU{2)l X SU{2)r is an even better symmetry. In this case the unbroken SU{2)y 
subgroup is isospin. Chiral symmetry has important implications for the interaction 
of pions, kaons, and the eta with each other as well as with the heavier hadrons. 

In the limit itlq —>■ oo the heavy quark sector in Eq. (|2.1| ) also exhibits additional 
symmetries. Consider a heavy quark with momentum p = rnqv which interacts 
with a gluon with momentum /c, so that the final momentum of the heavy quark is 
p' = rnqv + k. For k ~ Aqcd the velocity v of the heavy quark is conserved up to 
small terms of order Aqcd/'"^q and becomes a useful label for the heavy quark field. 
To construct a Lagrangian with a good mq — > oo limit we set 

Q[x) = e-*™«"-"/i„(x) + . . . , (2.5) 

where + 'f)) = and = 1. Momenta of order Aqcd cannot produce a heavy 
anti-quark, so in the heavy quark sector the anti-particles can be integrated out as 
indicated by the ellipsis in Eq. ( p.5|) . The number of heavy quarks is then conserved 
in the effective theory. After some straightforward algebra, the Lagrangian for heavy 



quarks with velocity v becomes] 40 



C^ = Kiv- DK +0{l/mq). (2.6) 
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At leading order this Lagrangian is independent of the spin and flavor of the heavy 
quark. For heavy quarks we have an SU {2N)^ symmetry, known as heavy quark 
symmetry |2l]]. This symmetry has important implications for exclusive and inclu- 
sive decays involving hadrons with a heavy quark, such as the D,D*,B,B*,Ab,Ac, 
etc. For a more detailed discussion of the implications of heavy quark symmetry see 



Ref. [41 



The next few sections explain the implication of these symmetries in the formula- 
tion of the low energy effective field theories for interaction of the light pseudoscalars 
TTjKjf], with heavy hadrons including the charm and bottom mesons D,D*,B,B*, 
vector mesons p,K*,(j),uj, and nucleons N = p,n. For reviews of the formalism for 
zero and one heavy particle see [0, For effective field theory with two heavy 

particles see 



2.1 Chiral symmetry and chiral perturbation the- 
ory 

In this section SU{3) chiral perturbation theory is reviewed. The formalism for the 
more accurate but less predictive SU{2) chiral perturbation theory follows in a sim- 
ilar manner. The eight pseudo-Goldstone bosons vr* that arise from the breaking 
SU{3)l X SU{3)r SU{3)v are identified with the observed light pseudoscalar 
mesons, , n'^ , K"^ , K'^ , K'^ , and t]. These will be encoded in the exponential repre- 
sentation 



S = ^2 = exp 



(2.7) 



where E^E = C,^^ = 1 and the A* are 3x3 matrices such that 



n = Ti'X' 



TT 



TT 



K- 



-TX 



{21 
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In Eq. ( |2.7| ) / ~ = 131 MeV, where is the pion decay constant, 



{0\ur75d\n-ip)) = -tf^pr (2.9) 
Note that 11 transforms as an octet under the unbroken SU{3)v- The ^ field in 



Eq. ( p.7| ) will become useful when we add heavy matter fields in section |2.2.1| . Under 
an SU{3)l x SU{3)ji transformation 



^ ^ L^f/t = ;7^i?t ^ ^here U = V^RWRJ:^L^ . (2.10) 

The non-zero quark masses rUg = diag(m„, m^, m^) break the chiral symmetry. 
To include rriq in our low energy Lagrangian we need to write a term that transforms 
in the same way as the light quark mass term in £qcd in Eq- ( |2.1| ). To do this we 
pretend that rrig LirtqR'^ under a SU{3)l x SU{3)r transformation, and then form 
invariants with m^. The Lagrangian with the fewest derivatives and powers of 
that satisfies the symmetry constraints is 



£2 £2tj 

= ^Tra^Sa^St + ^^Tr(m,S + mgEt), (2.11) 



where J^Bq/A = f in Eq. ( |2.4|) . Note that expanding the E fields in terms of 11 gives a 
canonically normalized kinetic term for 11 plus an infinite number of interaction terms 
with determined coefficients. In the interaction terms the pseudo-scalar fields are 
derivatively coupled, which is a general feature of chirally invariant couplings involving 
n. This follows from the fact that constant goldstone boson fields vTj are a rotation of 
S, and correspond to an equivalent vacuum for the spontaneous symmetry breaking. 
When neglecting isospin violation it is conventional to define m = (m^ + mrf)/2. In 
this case the meson masses are 

ml = 2Bom , m|- = Bo{m + m^) , = ^^(2ms + m) . (2.12) 
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Therefore, one power of a quark mass corresponds to two powers of a meson mass. 
This is the most general leading order behavior, given that these squared masses have 
a Taylor series in rrig, and that a constant term is forbidden by the fact that in the 
limit rrig the mesons are massless Goldstone bosons. 

The Lagrangian in Eq. ( p.ll| ) is not the most general one that is invariant under 
the desired symmetries. In particular, we can add an operator with dimension m 
that involves more derivatives or powers of rrig and a coupling of dimension 4 — m. 
After using the equations of motion^ there are 10 linearly independent terms with 



dimension coefficients |5y, |5^. For example 



4^) = ai [Tra^S9^St]' + ... . (2.13) 

Couplings like ai encode information about the short distance physics which was 
integrated out, so their scale is set by short distance scales like the chiral symmetry 
breaking scale A^. If p is a typical momentum, then higher dimension operators 
are suppressed by powers of p^/A^, m^/A^, and m|^/A^. This is the chiral power 
counting. It is convenient to consider ~ ~ mj^ ~ niq and then call Eq. ( p.ll|) 



the 0{p'^) Lagrangian. Since the particles in are relativistic, E"^ = p"^ + m^, and 
counting powers of the energy and powers of momenta are equivalent here. 

Along with higher dimension operators we must also consider loop corrections. 
These corrections are necessary, for instance, to restore unitarity to the S-matrix. 
The chiral power counting can also be applied to loop diagrams. Consider a graph 
with L loops, and vertices that are Weinberg |^ proved that this diagram 

is 0{p^) where 

D = 2(L + 1) + ^(m - 2)r2„ > 2. (2.14) 

m 

Each additional loop adds two powers of p. Instead of remembering the formula 
in Eq. ( p.l4|) we can power count an arbitrary loop graph by assigning appropriate 



^Note that after using the equations of motion the remaining Lagrangian can also be used in loop 
calculations ||^, ||, ||. 
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powers of p to the vertices, a for each loop integration, and factors of l/p^ for 
each propagator. If loop integrals are regulated in a mass independent way, then 
this counting is not affected by the divergences. It is convenient to use dimensional 
regularization where we continue the dimension of space time to c? = 4 — 2e. At 
O(p^) we must include the 0{p'^) vertices at tree level and the 0{p'^) vertices at 
one loop. Because of renormalization, these two contributions can not be separately 
specified in a unique manner. A logarithmic divergence in a loop integral induces a 
l/e + ln(/i^/p^) dependence in the result. The 1/e pole is subtracted or absorbed into 
a 0{p'^) coefficient, so in this sense these couplings act as counterterms. The fact 
that divergences are polynomial in the mass or momentum squared |Q along with 



the chiral power counting implies that divergences can always be absorbed in this 
way. This theory is said to be renormalizable order-by-order in the power counting. 
The finite (e independent) part of the 0{p'^) coupling depends on /x in such a way 
that it cancels the dependence from the loop. Changing /i changes the value of the 
loop with a compensating change in the value of the 0{p^) coupling. It might seem 
strange that the /z dependence exactly cancels (which is different than the situation 
in perturbative QCD where the /i dependence only cancels to a given order in ^^(/i)). 
However, this is nothing more than the statement that if we could calculate the full 
amplitude then it would be independent of the renormalization point. Therefore, 
expanding this amplitude in a power series in p^ and iriq gives coefficients which are 
termwise independent of 

The /i dependence of the loops and counterterms gives us a method for determining 
the size of A-^, called naive dimensional analysis [Q. Consider the graphs for tttt 
scattering at 0{p'^). For simplicity we use SU(2) chiral perturbation theory. Setting 
constants of order unity equal to 1, the amplitude takes the form 

m 



P + /2(47r/)2 



(^)+^]+7l"i(/^) + ---' (2-15) 



where is a number and p is the center of mass momentum. At tree level the first 
term in Eq. (|2.11|) gives an order p^ contribution which is the first term in Eq. (|2.15| ). 
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The second term is the contribution from the loop graphs, and the factor of l/(47r)^ 
arises from the loop integration. The third term in Eq. ( p.l5| ) is the contribution 
of the operator in Eq. (|2.13|) . For simplicity, terms with dependence have been 
left out as indicated by the elhpses. If the value of /x is changed then the second 
and third terms change in size while the sum stays the same. If there is no fine- 
tuning of parameters, then Q;i(/i) must be at least as large as the change to the loop 



graph induced by rescaling /i by an amount of order 1 ||5^. Thus, naive dimensional 
analysis implies that the second and third terms will be roughly the same size, and 
ai{fi) ~ l/(47r)^ ~ 0.006. It also implies that a natural size for the chiral symmetry 
breaking scale isQ ~ 47r/. When the 0{p'^) couplings are fit to data, they are found 
to be ~ 10^^ with /i = rup, which agrees with the dimensional analysis argument. 
The choice of /i reflects the fact that ai(yu) knows only about short distance scales 
(and in particular is independent of m^) so 

«i(/^)--(^[ln(^)+i^'], (2.16) 

where K' is a constant and A ~ rrip or A^. To avoid large logarithms in the coefficients 
we pick yU ~ A. This leaves potentially large logarithms in the matrix elements. 
In (/i^/m^). However since the theory is finite in the chiral limit — > these come 
multiplied by a power of m^, so although they are enhanced relative to other (9(m^) 
terms, they are not particularly large. 



2.2 Dynamics with one heavy particle 
2.2.1 The D''*^ and B^*'> and heavy quark symmetry 



The use of heavy quark symmetry] 20, 21] results in a dramatic improvement in our 



understanding of the spectroscopy of hadrons containing a single heavy quark. In 

the limit where the heavy quark mass goes to infinity, mq oo, such hadrons are 

•^Some authors use = 93MeV rather than / ~ 131 MeV, so = At:F.^. Dimensional analysis 
can not tell the difference. 
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classified not only by their total spin J, but also by the spin of their light degrees 
of freedom (i.e., light quarks and gluons), si [Q. In this limit, hadrons containing 
a single heavy quark come in degenerate doublets with total spin, J± = si ± ^, 
coming from combining the spin of the light degrees of freedom with the spin of 
the heavy quark, sq = |. (An exception occurs for baryons with si = 0, where 
there is only a single state with J = |.) The ground state mesons with Qq flavor 
quantum numbers contain light degrees of freedom with spin-parity sj'- = ^ , yielding 
a doublet containing a spin zero and spin one meson. For Q = c these mesons are 
the D and D*, while Q = b gives the B and B* mesons. The observed doublets are 
indeed very close in mass, m^i = 1.867 GeV, mo* = 2.008 GeV, uib = 5.279 GeV, 
and TUB* = 5.325 GeV|^. The heavy quark flavor symmetry gives further relations 
between the D^*^ and B^*\ 

The heavy mesons come in triplets under the SU{3)v symmetry, {D^, D^, Ds), 
p*o^ D*), {B-, fiO, B,), and (5*", B*). We will use the dimension 

3/2 HQET velocity dependent fields pP{v) and Pa^'^^^\v) (a=l,2,3), where p(^)(i;) 
destroys a D with velocity v, etc. It is convenient to include Pa and in a 4 x 4 
matrix 



p<Q),^^ - PP 75 



(2.17) 



2 

ni'^^ transforms linearly under both a heavy quark spin transformation -D(-R), and 
under a heavy quark flavor transformation U G SU (2) 



H(Q) D{R)Hfi^ , ^ UijH'f^^'^ , (2.18) 

and satisfies ^if^ = Ha = —Hai^. The conjugate field is defined as Ha = 7°-ff]7°. It 
is also convenient to define vector and axial vector currents, 

^(^t^M^ + ^^M^t)^ and A^='-ied''^-^d^e), (2.19) 



which contain an even and odd number of 11 fields respectively. Under a SU{3)l x 



SU{3)ji transformation 
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Ha^HbUl, V ^ UV^'W + zUd''U^ , and A^" ^ UA^^UK (2.20) 

The lowest order Lagrangian invariant under these symmetries is 

= -TTHatv-D,aH, + gT:THaH,^^^,Al, (2.21) 

where the chiral covariant derivative is = Sabd'^ — V^'^. The Ha propagator 
derived from the kinetic term in Eq. (p.21|) is often referred to as static, since in the 



rest frame v = (1,0) the equations of motion give zero energy for an onshell particle. 
(Recall that analagous to Eq. ( |2.5| ), the v dependent fields already have a factor of 
rriH subtracted from their energy.) Like in Eq. (|2.11 ), Ch is organized by an 



expansion in derivatives and powers of nig. Ch also involves an expansion in powers 
of l/rriQ, where terms at order l/mq break heavy quark symmetry. 

Since ~ A^' ~ p the Lagrangian in Eq. ( p.21| ) is 0{p). It contains one coupling 
g for p(Q)*p(Q)n and p(Q)*p(Q)*n. This coupling will be discussed in greater detail 
in Chapter 3. The propagators for the heavy pseudo-scalar and vector mesons, 

^Sab -^6a, {g -v v.) ^ ^^.22) 



2{y-k + ie)' 2{v k + ie 

are 0{l/p). When power counting graphs, the meson propagators give l/p^ and the 
loop measure gives a as before. Because of the form of the heavy propagators and 
couplings, the power counting involves powers of p and 171^^. Higher order corrections 
to the Lagrangian in Eq. ( p.21| ) are discussed in chapter 3. 

2.2.2 The heavy- vector meson chiral Lagrangian 



In this section we extend the heavy matter formalism to the vector mesons, p^, p°, 
0, uj, K*^, K*^, and K*^, following the presentation in Ref. [^. One might ask if 



these mesons should be treated relativistically; the lightest has mass nip = 770 MeV 
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which is starting to approach the low momentum regime of interest. However, if these 
particles are not treated as heavy then predictive power is lost. Unlike the pion, the 
vector mesons are not pseudo-Goldstone bosons, so they do not have to be derivatively 
coupled and their self-interactions are not constrained by chiral symmetry. When 
they are treated as heavy, the interaction terms in the Lagrangian can be expanded 
in derivatives giving an expansion in powers of p/mp. 

The vector meson fields are introduced as a 3 x 3 octet matrix and a singlet 



5. 



k: 



\ 



/,(o) 



(2.23) 



where v ■ O = v ■ S = 0. The dependence of these fields on the fixed four- velocity v 
has been suppressed. Under SU{3)l x SU{3)r the fields in Eq. ( ^.23| ) transform as 



(2.24) 



The 0{p) Lagrangian is ||58[ 



C 



V 



- SUv d - TiOl ivDO" 

+zg. Si Tr(aAA) v^e^""'^ + h.c. + ig^ Tr({Oi, 0,}Ax) t;.e^'^'^ 



(2.25) 



where the chiral covariant derivative is D'^O^ = d'^O^ + \y'^ , O^] and and are 
given in Eq. ( p.l9| ). The octet and singlet originally have masses /io and /is- When 
the velocity dependent fields are constructed we rescale both and by a common 



factor, A/2/i8 e'^^^'"'^. This leaves a term involving the mass difference, A/i = /^q — /is < 
200 MeV which may be treated as order rriq. Corrections to Eq. (p.25|) involving the 
quark mass matrix rriq induce mass differences between the vector mesons. The mass 
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eigenstates of the (j)^^^ — (p^^^ mass matrix are 



sin^ - cos^ , \uj) = cos^ + sin^ 



(2.26) 



where the SU{3)v prediction for the mixing angle is tan^ ~ ±0.76. 



Further predictive power can be obtained by considering the hmit of large A^^^ [58 
In this limit Afi = 0, tan^^ = l/\/2 and the octet and singlet mesons can be combined 
into a single nonet matrix 



N„ 



V 



pI/V2 + uj,/V2 
K*- 



K*+ 


\ 



















(2.27) 



At leading order in A^^ the Lagrangian in Eq. ( |2.25|) becomes 

Cv = -T:rNUvDm + tg2TTi{NlN,}Ax)v^e''''^'' . 



{2.2t 



Chiral and heavy quark symmetries can be used to relate the form factors describ- 
ing the semileptonic decays: D K*ii'£, D ^ piui, B ^ K*ii, and B pii^e- 
In chapter 4 the heavy vector meson formalism described in this section will be used 
to estimate symmetry breaking corrections to the heavy quark and chiral symmetry 
relations between these form factors. 



2.3 Dynamics with two heavy particles 
2.3.1 Two nucleon effective field theory 

This section considers an effective field theory for two heavy particles. The application 
in chapter 5 involves nucleon-nucleon scattering, so the particles are taken to be 
nucleons. For processes with one nucleon a formalism similar to that in Section |2.2.1 



may be used. For simplicity, pion-nucleon interactions will not be considered in this 
section, but will be considered in chapter 5. This simplification will allow us to 
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emphasize the quahtatively new features of the two-nucleon theory. Without pions 
the effective field theory is a valid description of nucleon interactions for momenta 
p <^ m^. 

Below the scale m^r, the pion can be integrated out, leaving a theory of non- 
relativistic nucleons interacting via contact interactions. The nucleon field is a 
doublet under isospin. The full Lagrangian in the two nucleon sector is given by: 



NN 



.a, + ^V(2M) + ...jiV-5^5^c£o£l + ... 



(2.29) 



s m=0 



where M is the nucleon mass, and the ellipsis refers to relativistic corrections. The 
transformation to non-relativistic fields is analogous to Eq. ( |2.5| ) where here it is 
convenient to choose v = (1,0). In Eq. ( |2.29| ), is an operator with 2m spatial 
derivatives and four-nucleon fields. We will work in a basis in which these operators 
mediate transitions between ingoing and outgoing two-nucleon states of definite total 
angular momentum. The superscript s will give the angular momentum quantum 
numbers of these states in the standard spectroscopic notation, '^^'^^Lj. States with 
(— l)'^"'"'^ even are isospin triplets, while those with (—1)"^^^ odd are isosinglets. If 
we denote the incoming and outgoing orbital angular momentum by L and L', then 
any operator mediating a transition between these states must contain at least L + L' 
derivatives. For states with 5* = 0, |L — L'\ = 0, while for states with S = 1, 
\L -L'\=0 or 2. 

In this section only S-wave transitions {L = L' = 0) will be discussed. For s = ^So 
or ^Si the first two terms in the series are 



E 



C'W(jyrp,(.)^)t(iv^p.(-)iV) 



(iV^p.(^)iV)t(iV^p.(^) N) + h.c. 



(2.30) 



+ 



where the matrices P/^^ project onto the correct spin and isospin states 



Pf''^^ = 4= i^^2) (^r^r,) , pf^^^ = (ta.a,) {ir,) . (2.31) 



The Gahlean invariant derivative in Eq. (|2.30|) is = V — 2V • ^ + ^ , and the 



elhpsis denote contributions with more derivatives. The normahzation in Eq. ( ^.30[ ) 
imphes that between S-wave states the Feynman rules are 

Co 

-tCo, (2.32) 



[{Pi - P2Y + (P3 - Pif] = -i C2P^ 



where the last equality is true when the nucleons are onshell in the center of mass 
frame, and p is the magnitude of the center of mass momentum. 

In a theory with two heavy particles using a static propagator is problematic. The 
loop graph 

=i-^Cor[-^— (2.33) 

has a pinch-singularity in the go integration at small go. This infrared singularity 
indicates that the static propagator is missing some essential physics. In theories 
with two heavy particles the kinetic energy term in Eq. ( p.29| ) becomes a relevant 
operator of the same order as the dt term. Including this term removes the singularity. 
The equations of motion for external nucleons are then po = P^/(2M), so for power 
counting, the nucleon energies and momentum are not equal in size. In dimensional 
regularization the loop graph in Eq. ( p.33|) is finite 
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. ,2 MV-ME - te 
= -^(^o) = -r{C,) [-^). 

Here E = p'^/M is the center of mass energy. In the second hne of Eq. (|2.34|) the qo 
integration was done by contour integration. To power count this graph note that 
~ p'^/M, g ~ so d^q ~ p^/M, and each of the nucleon propagators gives an 
M/p'^. The graph is therefore order p in agreement with Eq. ( p. 341) . 



The resuh in Eq. ( |2.34| ) has a factor of the nucleon mass M in the numerator. Since 



each loop with two nucleons gives an additional factor of M one might worry that these 
large factors will spoil the power counting. The reason for using a non-relativistic 
expansion in the first place was that each graph scales as a definite power of M, so we 
can keep track of these large factors^. From Eq. ( p.29| ) the coupling Cq has dimension 
—2. To count factors of M we must determine how the couplings scal^ with M. 
To determine this, rescale all energies, g° —>■ q'^/M, and time coordinates, t —>■ Mt, 
so that dimensionful quantities have the same size (ie., are measured in units of p). 
If we demand that the action is independent of M, then since the measure d'^x ~ M, 
the Lagrange density C ~ 1/M. The kinetic term determines that our nucleon fields 
scale as N{x) ~ M°, so from Eq. ( p.29| ) the coupling 

C2m-^l/M. (2.35) 

With the M scaling for the couplings determined, the scaling of any Feynman graph 
can be found. A nucleon propagator gives one power of M , and each momentum 
space loop integration gives a 1/M. For bubble graphs that have insertions of the 
four- nucleon operators, Np = Nl + Ny — 1, where Np, Nl, Ny are the number of 
propagators, loops and vertices. Thus, when relativistic corrections are neglected any 



^What we are keeping track of is the explicit M dependence in the Lagrangian. The Lagrangian 
does have further imphcit M dependence since M is a function of Aqcd and the scales that appear 
in the short distance coupHngs depend on Aqcd as weh. For this reason saying we know the M 
dependence of an ampUtude is not as strong a statement as saying we know the m,r dependence. 

^If the contact interactions were replaced with Coulombic photon exchange then the interaction 



would not involve any powers of M. In this case the graph in Eq. (2.34) would scale as a^M/p^ ~ 
/{p^v) where a — /An is the fine structure constant. This is a factor of a/v times a single 
photon exchange. For a Coulombic bound state a/v '-^ 1. Summing the most singular a/v terms is 
equivalent to solving the Schroedinger equation in a Coulomb potential. 
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graph built out of the interactions in Eq. ( p.29|) scales as since Np — — Ny = 
— 1. Therefore, the 2^2 scattering amplitude A ~ 1/M. With the definition of 
A used here, this scaling gives a finite cross-section in the M — > cxd limit which is 
physically sensible. (Note that there is a limit of QCD where M — * oo but Aqcd is 
finite, namely large N^. with Ncag held fixed |]5^.) Since all graphs scale the same 



way with M, M is irrelevant to the power counting. Relativistic corrections are 
included perturbatively0, and are generally suppressed by p^/M^ relative to the 
leading contribution to an observable. 

Applying dimensional analysis to the short distance coupling constants now gives 

11 1 

where A denotes the scale of short distance physics that was not included exphcitly 
(and we are assuming A < M). Treating the C2m couplings perturbatively gives an 
expansion in p/A. In the current case we expect that A ~ m^r. 

In nature, however, the dimensional analysis in Eq. ( |2.36| ) fails. This is because the 



nucleon-nucleon system is fine tuned to have bound states near threshold 113, H, p. 
In the ^5*1 channel this bound state is the deuteron with binding energy B = 2.2245 ± 



0.0002 MeV. This energy corresponds to the momentum 7 = V MB = 45.7 MeV ^ 
Aqcd- This bound state gives a pole in the two-to-two scattering amplitude which 
can not be reproduced perturbatively and therefore limits the range of the momentum 
expansion to p < 7. In the ^Sq channel the situation is even worse. In this case there 
is a pseudo-bound state sitting 8 MeV above threshold which limits the momentum 
expansion to p < 8 MeV. These bound states are related to the occurrence of unnat- 
urally large scattering lengths. Recall that low energy scattering can be described by 
an effective range expansion 

pcot<5(^) = -4y + ^r(^^p' + -- - , (2.37) 
where (5'-''^ is the phase shift, a*-**^ is the scattering length and Tq*'' is the effective range. 
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M 



In quantum mechanics, with a well-behaved potential, it can be shown that Tq is of 
order the range of the interaction, tq ~ 1/A, but a*^*-* can differ from A by orders of 
magnitude ||62|. From the formula for the amplitudes^ 



•^^^^ = S ^r^— ^ (2-38) 

Mpcot5W-«p ^ ^ 

it is straightforward to show that large negative or positive scattering lengths corre- 
spond to bound states just above or below threshold. Experimentally [|63| , 



^CSo) ^ -23.714 ± 0.013 fm, and a^'^^) = 5.425 ± 0.001 fm, (2.39) 

or = —8.3 MeV and = 36 MeV. These values of 1/a are small relative 

to and Aqcd- 

If the scale for the bound states or scattering lengths were set by A then the 
scaling in Eq. ( p.36| ) would be correct. The next section will discuss how the C2m 
coefficients scale for a theory with large scattering lengths. This power counting was 
worked out by Kaplan, Savage, and Wise [^, |[ (KSW). 



2.3.2 Power counting and ultraviolet fixed points 

In this section we will examine how unnaturally large scattering lengths affect the 
importance of four-nucleon operators. We also explain how linear ultraviolet diver- 
gences play a role in determining the KSW power counting, and why it is useful to 
recast this in the framework of the renormalization group. 

Treating the effective range term and higher powers of p"^ in Eq. ( ^.37| ) as pertur- 
bations, the amplitude is 

A^'^ = -^ ,, . , (2.40) 



^Strictly speaking Eq. ( 2.38| ) only holds in the channel. The ^Si channel is more complicated 
because of '^S'l — ^ Di mixing, but this mixing is a small effect which we will ignore for the time 
being, but discuss it in Chapter 8. 
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Figure 2.1: The leading order contribution to the A^A^ scattering amphtude. 

and has a good a ^ oo hmit. The pole at p = i/a^^^'^'^ corresponds to the deuteron 
bound state. The amplitude in Eq. ( p.40| ) is reproduced in the effective field theory 
by summing the loop graphs in Fig. p.l| using the result in Eq. (|2.34| ) 



oo 



a = -,c„V(-'P^g°)-= (2.41) 



Matching onto the effective range expansion then gives the values 

- is) 47ra(^) 

Thus, the Co coupling must be fine tuned to a large value to reproduce the observed 
scattering length [ 28 1. This result depends on our definition for the renormalized 
coupling constant Co, or in other words our choice of renormalization scheme (which 
is distinct from choosing a regularization method). The bars in Eqs. ( p.41|) and 
( p.42| ) indicate that only divergent terms have been subtracted, which is the Minimal 
Subtraction (MS) scheme. The sum in Eq. ( |2.4lD only converges for p < 1/a. Each 
loop graph in the sum grows with p, and the radius of convergence is 1/a. For p > 1/a 
we analytically continue and use the result on the right hand side of Eq. ( p.41| ). 
Looking back at Eq. (|2.34|), this one-loop graph is linearly divergent with d = 4, 



however in dimensional regularization this power divergence is not present. If we had 
calculated this loop graph with a finite momentum cutoff L then we would have 

XX ='(^») i X?i^^¥^rXX{'^^-^-)' 

where the ellipses denote terms that vanish as L oo. Defining a counterterm 
to cancel the term proportional to L, and taking L ^ oo gives back the result 
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in Eq. ( |2.34| ). Note that the hmit p oo and the integral over q do not commute 
because the integral is divergent. The growth of Eq. ( |2.43| ) with p is tied to this linear 
ultraviolet divergence. It is useful to use our freedom in defining the renormalization 
scheme to keep track of this. In the power divergence subtraction scheme (PDS)pT| 
additional finite subtractions are made in one-to-one correspondence with the linear 
divergences. In this scheme the coupling constants depend on the renormalization 
point hr and the value of the loop in Eq. (|2.43| ) becomes^ 
Co Co 



iM 



ip + fiR 



(2.44) 



In fact, the same result is obtained in a more physical scheme where the renormalized 
coupling is defined to be the value of the four point function evaluated at p = ifiR 
p7| , |6^ , |35| , |36| . These renormalization schemes will be discussed in greater detail in 
Chapter 5. 



In PDS, the renormalized coupling is [|3T 



R) 



M fiR- 1/aW 



(2.45) 



In this scheme the fine tuning as a — > oo, is that Co(/i/j) gets closer to its fiR ^ oo 
value. Cq^hr) has the beta function 



a M^R 2 

PO = /^i?^ Co(/iij) = — — Cq[iir) 

dfXR Arc 



(2.46) 



The renormalization group scaling gives us information about the behavior of the 
theory at the scale fiR. For fiR p <ti l/a, Cq^hr) behaves like a constant and 
pCo{fiR) can be treated perturbatively. For fiR p > 1/a, pCo(yUi?,) ~ 1 and the sum 
in Eq. (|2.41|) must be done. The factors of fiR, make each term in the sum roughly the 
same size which is good from the point of view of power counting. For fiR ~ p ^ 1/a 
the summation of Cq bubble graphs is always necessary and should be considered to 



''The notation fiji is used for the renormaUzation point in this section to agree with the notation 
used in chapter 5. 
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be infrared physics that has been built into the theory. In Ref. it was shown that 
the scahng in Eq. ( |2.45| ) is also reproduced with a Wilsonian renormalization group 
approach. 



The significance of the beta function in Eq. (|2.46| ) is shown more clearly in 



Fig. |2.2|°| . To make the beta function dimensionless it has been rescaled by MA/(47r) 



- MA ^ ^ a llr 

The prefactor a A can be safely ignored. The important dependence is in a fi^ since 
this factor measures the scale of interest relative to the scattering length. To study 
values of a fi^ from — oo to oo, a variable x is defined to map this range onto a compact 
interval. 



afiR = tan(^). (2.48) 



The dependence of f3{Co) on x is plotted in Fig. |2]^. Consider fixing the value 
of fiR > and varying the value of a. The points a = — oo, 0, and oo are fixed 
points of the beta function. Classically this makes sense since the scattering length 
is a measure of the interaction size. For a ~ or ±oo the size is so small or big 
that the interaction looks the same at all scales. In fact a = is a trivial non- 
interacting fixed point, whereas a = ±oo are non-trivial interacting fixed points 



where the theory is scale invariant at lowest order. Another feature in Fig. |2.2| is that 
/3 — ^ oo for fiR = 1/a > 0. This corresponds to the deuteron bound state. Performing 
perturbation theory about a = we can never describe this bound state, so we are 
limited to describing the region fiR < 1/a. If perturbation theory is performed near 
a = oo, then the deuteron is a physical state in the spectrum of the theory. If we are 
interested in physics at fiR ^ then the observed ^Sq and scattering lengths 
place us at the location of the stars in Fig. |2.2| . Looking at the distance along the x 
axis we are much closer to the a = ±oo fixed points than to a = 0. 



^This figure was inspired by a similar plot shown by David B. Kaplan in a physics colloquium. 
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For 1/a powers Q ^ p ^ fiR are counted instead of just powers of momenta. 
The graphs in Fig. ( |2.1| ) are leading order and all scale as l/Q. In the PDS renor- 
malization scheme this power counting is manifest. By solving the beta functions for 
the coefficients of operators with more derivatives it can be shown that 



for p ~ ^R ^ 1/a. Since insertions of these operators scale with non-negative pow- 
ers of Q, they may be treated perturbatively. This renormalization group scaling is 
also important in determining the power counting of operators that involve four nu- 



corresponds to the pole in our scattering amplitude, deuteron properties can be sys- 



When pion interactions are included they enter at order Q^, which is one higher 
order than insertions of Cq. Therefore, the discussion in this section is also relevant to 
the leading order theory with pions. In the theory with pions the physics encoded in 
the scale A in Eq. (|2.2|) changes. However, the power counting for the coefficients 
remains the same because pion effects are subleading corrections to the running of 
these operators. Pion interactions in the two-nucleon theory will be discussed in more 
detail in chapter 5. 




(2.49) 



cleon fields and fields like the photon |3^, |6^ and pion. Note that since the deuteron 



tematically calculated in this field theory P^. 
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Figure 2.2: Fixed point structure of the beta function for Cq{iir) in the PDS renor- 
mahzation scheme. 
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Chapter 3 Extraction of the D^Dir 
couphng from D* decays 



Combining chiral perturbation theory with heavy quark effective theory (HQET) 
gives a good description of the low energy strong interactions between the pseudo- 
Goldstone bosons and mesons containing a single heavy quark. In this chapter we 
extend the formalism in Chapter 2 to describe D* and B* decays at subleading order. 
Due to heavy quark symmetry there is one coupling, g, for D*D7i, 0*0*71, B*B7t, 
and B*B*'7i, and one coupling, (3, for D*Dj, D*D*'y, B*B'^, and B*B*'y at leading 
orderQ. The value of the coupling g is important, since it appears in the expressions 
for many measurable quantities at low energy. These include the rate B — > D^*'^TTih'i 
| p7| , |7D|, form factors for weak transitions between heavy and light pseudo- 

scalars [^, |2j, ^ |7^, ^ |7^, ^ |7|, decay constants for the heavy mesons 



80, Bl, B2l, weak transitions to vector mesons [83 , form factors for B 



D'^*Hh'i |8i, 85, 36|, and heavy meson mass splittings [3?, BS,|8^ (for a review see |90]). 
However, the value of g has remained somewhat elusive, with values in the literature 
ranging from ~ 0.2 to 1.0. Recently, a CLEO measurement of D*^ D'^j 
brought the experimental uncertainties down to a level where a model independent 
extraction of g is possible from D* decays. 

As a consequence of HQS, the mass splitting between D* and D mesons is small 
(of order AqQ^/mc), leaving only a small amount of phase space for D* decays. In 
the dominant modes, D* Dtt, and D* D'j, the outgoing pion and photon are 
soft, making the chiral expansion a valid framework. The branching ratios for D*~^ 
decay are n+ (67.6%), D+ 7r° (30.7%) and D+ 7 (1.7%) 0. A can only decay 



into D°7rO (61.9%) and (38.1^ 



since there is not enough phase space for 



"'^Where it is meaningful we use tt to denote any member of the pseudo-Goldstone boson SU(3) 
octet, and D* and D for any member of the triplets {D*'-\ D*) and {D'^ , Dg) with a similar 
notation for B* and B. 
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71^ . The D* decays predominantly to -0^7 (94.2%) with a small amount going 
into the isospin violating mode Dg 7r° (5.8%) Since a measurement of the widths 
of the D* mesons has not yet been made, it is only possible to compare the ratios of 
branching fractions with theoretical predictions. The ratio 

R+ = B{D*+ ^ D\+)/B{D*+ ^ D+7r°) (3.1) 

is fixed by isospin to be = 2\k^+\'^ /\k^o\^ = 2.199 ± 0.064 (where k^+,o are 
three momenta for the outgoing pions in the D* rest frame). This constraint is often 
used in experimental extractions of the branching ratios to reduce systematic errors. 

It is interesting to note that the quark model predictions] 92, 93| for D*^ and D*~^ 
decays agree qualitatively with the data. One can understand, for instance, why the 
branching ratio B{D*^ D^l) is small compared to B{D*^ D^'-y). In the quark 
model the photon couples to the meson with a strength proportional to the sum of the 
magnetic moments of the two quarks, /i2 = 2/(3 mc) — l/{'imd) for D*^ D'^'y and 
yUi =2/(3 nic) + 2/(3 m„) for D*^ D^'j. The rate for the former is then suppressed 
by a factor 



/ii 



^ / , n9 ('^f^/'^r — 1/2)^ 

= m./m, ^ y/ / \ ^ ^ 0.04 , 3.2 



where mass ratios appropriate for constituent quarks have been used, mu/m^ ^ 1, 
iTLd/rric — rriu/mc ~ 1/4. This suppression results from the opposite signs in fii and 
/i2, which in turn follow from the (quark) charge assignments and spin wavefunctions 
for the heavy mesons. 

In the quark model g = 1 and (3 ~ 3GeV~^, while for the chiral quark model 
g = 0.75 [^]. Relativistic quark models tend to give smaller values, g ~ 0.4 95 
as do QCD sum rules, c/ ~ 0.2 - 0.4 @, 0, [53 



Our purpose here is to use heavy meson chiral perturbation theory at one-loop to 
extract the couplings g and jS from D* decays. In other words, we wish to examine 
the sensitivity of a model independent extraction of g and (3 to higher order correc- 
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tions. For D* —>■ D'-f, analyses beyond leading order have included the heavy quark's 
magnetic moment which arises at 1/ rric [0, |100|| , and the leading non-analytic effects 



from chiral loops proportional to |9^ . ^/fn^ terms proportional to both ttik and 
m^r were found to be important. These effects do not introduce any new unknown 
quantities into the calculation of the decay rates. For D* D'y and the isospin 
conserving D* Dn decays the effect of chiral logarithms, mqln{fi/mq), have also 
been considered ||101|| . These are formally enhanced over other rriq corrections in the 
chiral limit, rriq 0, however, the choice of the scale fi leads to some ambiguity in 
their contribution. (This scale dependence is cancelled by unknown couplings which 
arise at order rriq in the chiral Lagrangian.) The isospin violating decay D* DgH^ 
has only been considered at leading order, where it occurs through rj — n^ mixing [|102| ]. 

Here the investigation of all D* decays is extended to one-loop, including symmetry 
breaking corrections to order rUq and l/rric. Further l/rric and rUq contributions 
considered here include the effect of nonzero D*-D and Ds-D^ mass splittings, and 
the exact kinematics corresponding to nonzero outgoing pion or photon energy in the 
loop diagrams. (Their inclusion is motivated numerically since m^o ~ m/j* — m^) ~ 
m/)^ — m/), and the decay ^ Z^tt'^ only occurs if m^)* —m/) > m^o.) To simplify the 
organization of the calculation these splittings will be included as residual mass terms 
in our heavy meson propagators. This gives new non-analytic contributions to the 
D* — s> Dti^ and D* — > D'y decay rates. (To treat the mass splittings as perturbations 
one can simply expand these non-analytic functions.) At order rriq there are also 
analytic contributions due to new unknown couplings which are discussed. These 
new couplings can, in principle, be fixed using other observables. We estimate the 
effect these unknown couplings have on the extraction of g and (5. 

The calculation of the decay rates to order rriq and l/rric is taken up in sec- 



tions and |3.2| . In section ^]3| we compare the theoretical partial rates with the 
data to extract the D*Dtt and D*D'~f couplings and discuss the uncertainty involved. 
Predictions for the widths of the D* and B* mesons are also given. Conclusions can 
be found in section O, 
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3.1 Da 71 decays 

In this section we construct the effective chiral Lagrangian that describes the decays 
D* Dn to first order in the symmetry breaking parameters rrig and l/rric. Going 
beyond leading order also involves including loops with the pseudo-Goldstone bosons. 
From Chapter 2 section 2 recall that the lowest order Lagrangian is 

£2 £2td 

Co = — Tr9'^S9^St + ^Tr(mgS + mgS^) 

8 4 

-Tr H,iv ■ DbaHt + gTi H^H^j^j^Al . (3.3) 

The last term in Eq. couples P*P'k and P*P*'k with strength g (where P = D,B) 
and determines the decay rate D* Dn at lowest order. At order nig ~ l/rric the 
following mass correction terms appear 

Cm = -^Tr Haa^'HaCJ^, + 2AiTr HaHmL + SA^Tr HaHaml, (3.4) 

where = ^{C.mq^'^ + ^'^rUqC,). The A'^ term can be absorbed into the definition 
of rriH by a phase redefinition of H. The A2 term is responsible for the D*-D mass 
splitting at this order, 

A = TTiD* — = —2\2/mc . (3.5) 

The term involving Ai splits the mass of the triplets of D and D* states. Ignoring 
isospin violation this splitting is characterized by 

6 = rriD* — niD* = uid^ — niD = 2Xi{ms — m) (3.6) 

where m = rriu = m^. For the purpose of our power counting 6 ~ rUg ~ l/m^ ~ A. 
The effect of these mass splitting terms can be taken into account by including a 
residual mass term in each heavy meson propagator. Since we are interested in decay 
rates we choose the phase redefinition for our heavy fields to scale out the decaying 
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particle's mass. For D*^ and D*^ decays the denominator of our propagators are: 
2v-kfoT D*^ and D*+, 2{v-k~6) for D*, 2{vk + A) for and D+, and2{v-k + A-6) 
for Dg. For the D* decays the denominators are the above factors plus 26. (If we had 
scaled out a different mass then the calculation in the rest frame of the initial particle 
would involve a residual 'momentum' for the initial particle, but would yield the same 
results.) This results in additional non- analytic contributions from one- loop diagrams 
which are functions of the quantities A/m.^- and S/rriT^., where m.^- G {m^, m^, 
Formally, m^. ~ rrig ~ A ~ 5 and one can expand these contributions to get back the 
result of treating the terms in Eq. (|3.4| ) as perturbative mass insertions. 

Another type of 1/ corrections are those whose coefficients are fixed by velocity 



reparameterization invariance [103, 82 



SCv = -77— Tr Ha{iD)l^Hb + —TlHci^lD^ac"^ ■ Aba - iv ■ Aac'^ba)Hb1fil5 ' 

2mQ niQ 



(3.7) 



The first term here is analogous to the HQET kinetic operator, Okin = q^^v {^DY hy, 



but written in terms of the interpolating fields Pa and P*^. In Eq. (|3.7|) the derivatives 
give powers of the heavy meson's momentum. There are also contributions from Okin 
that break the fiavor symmetry where the derivatives are order Aqcd- In conjunction 
with the HQET chromomagnetic operator, Omag = 2m^^^ % ^ai^G^^ h^, these contri- 
butions to the Lagrangian modify the dynamics of the heavy meson states. They 
give Aqcd/itT'Q corrections in the form of time ordered products with the leading 
order current [ |104|] , which induce spin and flavor symmetry violating corrections to 



the form of the 0*071 coupling. We account for these corrections by introducing the 
couplings gi and g2 in Eq. ( |3.8| ) below. The last term in Eq. ( pI7| ) contributes at 
higher order in our power counting since it is suppressed by both a derivative and a 
power of 1/ rric. 

Further terms that correct the Lagrangian in Eq. ( |3.3| ) at order ruq ~ 1/ rric include 
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HQ 




x X 




X X 




+ ^Tr HaH,^,^,Al + ^Tr Hai^i^MbAl + . . . , 




(3.8) 



where D^^A^^ = d^A^^ + [l^", A^jba and = 47r/. The eUipses here denote terms 



pion, as well as terms with {iv ■ D) acting on an H. For processes with at most one 
pion and H on-shell the latter terms can be eliminated at this order, regardless of 
their chiral indices, by using the equations of motion for H . The Kj coefficients contain 
infinite and scale dependent pieces which cancel the corresponding contributions from 
the one-loop D* Dtt diagrams. For the Ki and k'i terms only the combination 
will enter in an isospin conserving manner here. (The combination 
Ki — k[ will contribute an isospin violating correction to R^-) At a given scale fi, the 
finite part of k,^ can be absorbed into the definition of g. The decays D* —>■ Dn have 
analytic contributions from ki and K5 at order rrig. 

For rriQ = rric the term in Eq. ( |3.8| ) involving gi can be absorbed into g (this term 
only enters into a comparison with B* decays). The term g2 breaks the equality of 
the D*D7T and D*D*n couplings. Since we only need the coupling D*D*7i in loops 
we can also absorb (72 into the definition of g. Thus, our g is defined as the D*D7t 
coupling with l/mg corrections arising in relating it to the couplings for D*D*7i and 

The terms in Eq. ( ^.8] ) involving 62 and ^3 contribute to D* —>■ Dtt'^, entering in a 

fixed linear combination with the tree level coupling g of the form g — {62 + Ss)v -k/A^. 

These are ~ 10% corrections for the decays D* Dir. The energy of the outgoing 

pion is roughly the same for all three decays, f ■ /c ~ .144 GeV. Therefore, it is 
^The k'i term was not present in p3]. The factor Bq/A^ is introduced here for later convenience. 



linear in m 



|(,^mg^^ — ^^mg,^) which contribute to processes with more than one 
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impossible to disentangle the contribution of ^2,3 from that of g for these decays, and 
the extraction of g presented here will implicitly include their contribution. For other 
processes involving pions with different v ■ k these counterterms can give a different 
contribution. This should be kept in mind when this value of g is used in a different 
context. 

Techniques for one-loop calculations in heavy hadron chiral perturbation theory 
are well known and will not be discussed here. Dimensional regularization is used 
and the renormalized counterterms are defined by subtracting the pole terms 1/e — 
7 + log (47r). In doing this type of one- loop calculation an important integral is 



d^-''q -i b 



(27r)4-2^ (g2 -m'^ + i£)2{q-v -b + ie) {Any 



1 , , u\ ,Tn^ 
^ + ln^ +2-2F - 
e b 



(3.9) 

where l/e = 1/e — 7 + log (47r). F is needed for both positive and negative 6, so 

\x\ < 1 



X I I Va-^ 



f-tan-1^^ 



— In (x -|- Va^^ — 1 ) > 1 



(3.10) 



For 6 > the function F was derived in ||105| , |77| and agrees with the above formula^ 
For X = bjm < —\ the logarithm in Eq. ( p.lO|) has an imaginary part. This corre- 



sponds to the physical intermediate state where a heavy meson of mass vtin produces 
particles of mass mn + b and m. For the calculation here the imaginary part only 
contributes from F im.,^ j [d^ — A)), and was found to always be numerically insignifi- 
cant. Note that the real part of x F{\jx) is continuous everywhere, and differentiable 
everywhere except x = —1. Also F{1) = F{—1) = 0. 

The decays D°7r° and D*+ D+7r°, and D* D,7r° have decay rates F^, 



3t 



Eq. i ^J^ for F disa grees with |82| for x < 0. Their F{l/x) is even under x — > —x making 



Eq. (3.9) discontinuous at A = 0. Furthermore, their F has no imaginary part corresponding to the 



physical intermediate state. 
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b) 



c) 
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Figure 3.1: D and D* wavefunction renormalization graphs. The dashed hne repre- 
sents a pseudo-Goldstone boson. 

r^, and given by 



pa _ 9 



2 



12 71/2 



2 



K\'- (3.11) 



— » 

Here fc" is the three momentum of the outgoing pion, contains the vertex correc- 
tions, and = a/ 2'!,, contains the wavefunction renormahzation for the D* and 
D. When the ratio of to the D* D'-f rate is taken Z^^ will cancel out. However, 
Z^j does contribute to our predictions for the D* widths, where the ratio Z^^/Z^ will 
be kept to order g"^. The graphs in Fig. |3.1| give 



Z^D = l + ^^(Al,Ajj{[3m^-6(A + rfoniog(^)+ 3GiK,A + rfo)}, 

2 2 

Z%* = l + ^^(Aj,,Ajj{[3m^-4rf2-2(rfo-Aniog(£^) + 2GiK,rf^ 

+ 6-1^,^0- A)}, (3.12) 

where rrii is the mass of vr*, do = S^^6 for D*^ and decays and do = {6^^ — 1)6 
for decays. The notation in Eq. (|3.12|) assumes that we sum over b = 1,2,3 and 
i = 1, . . . ,8. The logarithms agree with ||101|| , except that we have kept terms of order 
A^ ~ c/q in the prefactor since these terms are enhanced for niq 0. Analytic terms 
of order A^ ~ d"^ are neglected since they are higher order in our power counting. 
The function Gi{a,b) in Eq. ( p.l2| ) has mass dimension 2, 



Gi(a,6) = -a^ + Ub^ ~-a^)F(a/b) + -(a^ -b^)^F'(a/b). (3.13) 
3 3 3 6 
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D 



Figure 3.2: Nonzero one-loop vertex corrections for the decays D*^ D^n^ and 
D*^ D^TT^ (a,b,c) and the pseudo-Goldstone boson wave function renormahzation 
graph (d). 

It contains an analytic part proportional to a^, and a non-analytic part which is 
a function of the ratio b/a. In the limit A — Eq. ( p.l2| ) gives Zo = Zo* in 
agreement with HQS. To obtain HQS in the finite part of the dimensionally regularized 
calculation of the graphs in Fig. ^]T] it was necessary to continue the D* fields to d = 
4 — 2e dimensions (so the D* polarization vector ea = (1 — ^)ia where ^ e*e" = —3). 

For a = 1,2 the decay proceeds directly so that at tree level Zl^^/Z^''^ = 1. At 
one loop we have non-zero vertex corrections from the graphs in Fig. |3.2| a,b,c. As 
noted in ||101|| , the two one-loop graphs that contain a D^*^D*nn vertex (not shown) 



vanish, and the graph in Fig. |3.2| c cancels with the vr*^ wavefunction renormahzation 
in Fig. 13d (this is also true for D*+ L>°7r+ and D* Z),7r°). For a = 1,2 the 
vertex corrections are 



+2Fi(mi,di,rf2) -4Fi(m„rfi,do) \ + qI{ki, k.,) , (3.14) 



where here d^ = S^^S, di = k - v + do, d2 = —A + do, and k is the outgoing momentum 
of the 7r°. The coefficient of the log (fi^/mf) term agrees with ||101|| . The function 
Fi in Eq. (p.l4|) has mass dimension 2 and contains both analytic and non-analytic 



parts, g'^^f contains the dependence of the rate on the (renormalized) counterterms 
ki{fi) and K^{fi). 

Fi{a,h,c) = J-a^ + -^^^[h{a'-h')F{a/h)-c{a'-c')F{a/c)\, 
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Figure 3.3: Nonzero vertex corrections for the decay D* Dg^^ which involve tt^ — r] 
mixing. The cross denotes leading order mixing while the triangle denotes mixing at 
next to leading order. 

a=l,2 _ < k.ifl) 

" (47r/)2 2 • ^^-^^^ 

We have ignored isospin violating counterterm corrections in In this case ^^'^ do 
not depend on K5, and furthermore are proportional to m^/(47r/)^, so these countert- 
erms are small. 

The decay D* DsTc^ is isospin violating, and the leading contribution occurs 
through f] — TT^ mixing [|102|] . To first order in the isospin violation the decay is sup- 



pressed at tree level by the mixing angle 6 = (1.00 ± 0.05) x 10 



1. = - "^^) = L (3 IQ) 

Zl 2(m,-m) 73 " 87.0' ^ ' ^ 

Beyond tree level we have corrections to the rj — mixing angle parameterized by 
Smix = 0.11 (Fig. |373|a), loop corrections to the 77— 7r° mixing graph (Figs. p3|b,c,d), 
as well as loop graphs with decay directly to that occur in an isospin violating 
combination (Figs. |3.2| a,b). The contribution of Fig. p.3| d is again cancelled by the 
pseudo-Goldstone boson wave function renormalization graph (Fig. p.2| d). Note that 
the decay D* — > Dgir^ cannot occur via a single virtual photon in the effective theory. 
In the quark model, decay to the spin and color singlet 7r° can occur if the single 
photon is accompanied by at least two gluons (with suppression a/Tr ~ 1/430 ||102|| ). 



We will neglect the possibility of such a single photon mediated transition here. Thus, 



1 (m„ - TTid) 



2 (nis — rh) 



n2 \i \8 \*t r ,,2 



l+6rm.+ log(^) mi + -{-di + d,d2 



2 



3 
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+dl - 2di do-2dl 



2 Fi(mi, di, d2) - A FAnii, di, d^ 



r(2 X* xi x*'!" 



fx 



log(-r^) + 2 Fi{mi, di, d2) - AFi{mi,di,do) 



mt 



(3.17) 



where for D* decay do = (5*^ — 1)S, di = k ■ v + do, and d2 = —A + d^. The tilde on 
the mass, rhi, indicates that isospin violation is taken into account. Also note that 



1 \«t™2 



b3"''i 



. In Eq. ( p. 171 ), the function gl depends on Ki 



and ^5, and at leading order in the isospin violation is 



[niu - rrid) 



2 ^ 



[m 



K5 



(3.18) 



In deriving this equation use has been made of = 23^171^ = 2i?ow^d = 2BQrh, 



m 



\ - ml/ 2 = Bonis, and m^± 



m 



(m„ - md)Bo. 



3.2 d: 



Da 7 decays 



To describe D* —>■ D'j, electromagnetic effects must be included, so the Lagrangian 
in Eq. ( pl3|) is gauged with a U{1) photon field B^. With octet and singlet charges, 
Q = diag(|, — |, — |) and Q' = | (for the c), the covariant derivative is ||106 



V^H = d^H + ieB^{Q'H-HQ)-V^H, 



where the vector and axial vector currents are now 



(3.19) 



(3.20) 
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However, this procedure does not induce a coupling between D* , D and without 
additional pions. Gauge invariant contact terms should also be included, and it is 
one of these that gives rise to the D*D'-f coupling (and a D*D*'~f coupling) 

Cp = ^Tr H^H^a'^-'F^^Ql. (3.21) 

Here P has mass dimension —1,Q^ = ^{^'^Q^ + ^Q^^), and F^j,^ = dfj,B^ — The 
terms which correct this Lagrangian at order rriq ~ l/rric have a similar form to those 
in Eq. (KEi) 



X X 
X X 

+ ^Tr HaH^^a^^Qltv ■ D,aF^, + ^Tr HaH^a^'QliD'^^^v^F,^ 

X X 



- -^Q' Tr Hacr^^'^HaF^, + ... . (3.22) 

The ellipses denote terms that do not contribute for processes without additional pions 
and/or can be eliminated using the equations of motion for H. For our purposes 
and in Eq. ( 3.22 ) are diagonal so only a = ai+a[ contributes. The finite part of as 



will be absorbed into the definition of (3. For itlq = rric, the f3i term can be absorbed, 
and we absorb the part of the f32 term that contributes to D*D^ since D*D*'y only 
contributes in loops for us. Thus, P is defined to be the D*D'y coupling at order 
l/rric. The last term in Eq. ( p.22|) is the contribution from the photon coupling to 



the c quark and has a coefficient which is fixed by heavy quark symmetry [ |10U| , 
This term is numerically important. However, here the leading order contribution to 
Daj is taken to be the ijiq ^ oo, nig ^ effect from Eq. (|3.21| ), so this l/rric 



term is part of the first order corrections. The Ti^2 terms are similar to the ^2,3 terms 
in Eq. ( p.8|) , and appear with (3 in the combination /3 — (ti + T2)v ■ k/A"^. Here Ti + T2 
will have an infinite part necessary for the one-loop renormalization. Again it is not 
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Figure 3.4: Nonzero vertex corrections for the decays D* D'y. 

possible to isolate the finite part of the (ri + T2) contribution from that of (3, so the 
extraction at this order includes the renormalized ri^2 with v ■ k 0.137 GeV. 
The decays D*^ D°7, D*+ D+-f, and D*^ 



Dsj have decay rates Tl^, F^, 



and given by 



"pa 

7 



a 
3 



Q 



aa I Q 



TTlr 



(3.23) 



where a ~ 1/137, is the three momentum of the outgoing photon, and the wave- 
function renormalization, Z^j, is given by Eq. (|3TT^) . To predict the D* widths. 



Z^f/Z" is kept to order g"^ and we take Z^^ x l/rric = l/rric. The vertex correction 



factor Z^ has nonzero contributions from the graphs in Fig. |3.4 Note that the two 
one-loop graphs that contain a D^*^D*n-]' vertex (not shown) do not contribute ||101|| . 



Furthermore, the graph in Fig. |3.4| b has no contribution from the D*D*-]' coupling 
which arises from gauging the lowest order Lagrangian in Eq. (|3.3|). Thus 



1 



.it 



1+ ^ KbQbbKa\ J ./^ 



2 r 2 

2 , ^/ j2 



mj H — [—di + did2 + d2 — 2di do — 2 do) 
3 



1 [X'\[Q,y 

(47r/)2 2Qaa 

(47r/)2 /5Q„, 
+e'"(ai,a5) , 



+2 Fi (mi, 6^1,^2) - 4:Fi{mi,di,do) 



log ( — 2) + ""^^ 



log (^)(c/o + — -) + F2{mi, do, k ■ v) 
mf 2 



(3.24) 



where is the charge of meson tt*, k is now the outgoing photon momentum, and 
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the di are as above (again they differ depending on whether it is D* or one of D*^ 
D*^ that is decaying). The coefficients of the mf log (fi'^/mf) terms agree with ||101 | 
The new function F2 has mass dimension 1, 



^2(0,6, c) 



2a2 r/(^+'=) , Fit) 

-26 -c / dt^^ 

c Ja/b t-^ 



-2b -c 



2a' 



Fix) 



1 F{x) 
4x^ ~ 2x^ ~ 4(x2 - 1) 



x=a/ (fc+c) 
x=a/b 



(3.25) 



It contains an analytic part proportional to 2dQ + v ■ k, and a non-analytic part which 
is a function of 6/mi and v ■ k/rrii. In Eq. ( p.24| ) g'^ contains the dependence of the 
rate on the (renormalized) counterterms ai(/i) and a^{fi). Assuming isospin to be 
conserved we have 



a=l,2 



ml ai 
2iAnfy 



{2m\ - ml) ai 
2 (Anf)^ 



{m\ - ml) as 

3gaa(47r/)2 ' 

{m\ - ml) as 



+ 



(4vr/)^ 



(3.26) 



By examining Eqs. ( |3.12|) , ( p.l4| ), (|3.17| ), and ( |3.24| ) we can get an idea of the 
size of the various one- loop corrections to and F^^. With our power counting 
A ~ (5 ~ f ■ A; ~ mq ~ so we can consider expanding in /S./mi, 5 /mi, and v ■ k/mi 



Gi{mi,b) 
Fi{mi,b, c) 
F2{mi,do,k- v) 



3 



6776 1662 
1 + ^ 



mt 



mt 



njb + c) ^ 16{b' + bc + c') ^ 

mi 9mf 
3dl + 3dok ■v + {k-v)'^ 



6mf 



+ 



(3.27) 



The leading terms in Gi and Fi are corrections to the rates. The second terms are 
order m'^^'^ and y/m^/ mc, and can be kept since they are unambiguously determined at 
the order we are working. The third and remaining terms in Gi and Fi are subleading 
in our power counting. The term -Tcmi in F2 is the formally enhanced contribution 
discovered in Note that there are no contributions to F2 proportional to 6 or 
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k ■ V. The second term in F2 in Eq. ( p.27| ) has contributions of order mq^"^, y/fn^k ■ v, 
and {k ■ vY / ^JffTq which again can be kept since they are unambiguously determined. 

The above power counting is sensible when mj is vhk or m^. We know that 
numerically ~A~5~/c-f, so for = the series in Eq. ( p.27| ) are not 
sensible. In the term — Trm^r in was found to be important, so we want to keep 
corrections with m^r dependence. Therefore, instead of expanding the non-analytic 
functions we choose to keep them in the non-analytic forms given in the Appendix. 
Numerically the one-loop corrections to V\ and Y\ are very small; with g = 1 they are 
of order ~ 2%. For F^, 6mix is a 11% correction to the tree level result in Eq. ( |3.16| ). 
Individually the terms proportional to g'^Fi and g'^ log (/i/mg) in Eq. (p.l7|) are ~ 10% 
corrections for g = 1. However, the loops graphs with rj — tt^ mixing tend to cancel 
those without 77 — 7?" mixing leaving a ~ 2% correction. The one-loop corrections to 



F" are larger, for instance the graph in Fig. |3.4| c gives sizeable corrections that are 
not suppressed by g^. Corrections to the coefficient of the leading g'^/jS term range 
from ~ 3% for D* and ~ 20% for D*° decay, to ~ 50% for the D*+. (The latter 
percentage is large because the only contribution for this decay comes from a charged 
pion in the loop of Fig. |3.4| d.) Corrections proportional to g"^ are only sizeable for 
D* — s> Dg^ where they are ~ 10% for g = I. 

3.3 Extraction of the couplings g and f3 

Using the calculation of the decay rates from the previous section, the couplings g 
and f3 can be extracted from a fit to the experimental data. Input parameters include 
rric = 1.4 GeV [|107| , |108|| , the meson masses from ||5^, A = rriD* — itld = 0.142 GeV, 



6 = m {,) —rrij^i*) = 0.100 GeV, and v-k which is determined from the masses. When 
isospin is assumed we use mx = 0.4957 GeV and = 0.1373 GeV. / is extracted 
from 7r~ decays. At tree level we use f = f\ = 0.131 GeV |^, while when loop 
contributions are included we use the one-loop relation between / and /tt ||52[ to get 
/ = 0.120 GeV. The ratio of the decay rates F" and F° are fit to the experimental 



numbers 
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^ D^-f)/B{D*^ -> D°7r°) = 0.616 ± 0.076 [57] 
B{D*+ D+^)/B{D*+ D+7r°) = 0.055 ± 0.017 [0 



B{Dl ^Ds'K^)/B{Dl ^ Dsi) = 0.062 ± 0.029 0, (3.28 



where the errors combine both statistical and systematic. Using the masses m^.o, 



mD*+, rriD*, and mass sphttings m£).o — m£,o, mD*+ — ?^d+5 "^dj — t^Ds from |^| 
gives the momentum ratios that appear in F^/F": 

\Ul\3 \t2\3 |p|3 



32.65 ±0.44, ^^ = 45.2 ±1.0, ^^ = 24.4 ±1.5. (3.29) 

\ki\' \kl\' \kl\' 



The errors here are clearly dominated by those in Eq. ( 3. 28 ). Equating the numbers in 



Eq. ( p.28|) to the ratio of rates from Eqs. ( p.ll| ) and ( |3.23|) gives a set of three nonlinear 



equations for g and f3 (where we ignore for the moment the unknown counterterms) . 
In general any pair of these equations will have several possible solutions. To find 
the best solution we take the error from Eq. (|3.28| ) and minimize the fo^^ the fit to 



the three measurements. We will restrict ourselves to the interesting range of values, 
< g < 1 and < /? < 6, discarding any solutions that lie outside this range. (The 
sign of g will not be determined here since it only appears quadratically in F" and 
F".) 

To test the consistency of the chiral expansion we will first check how the ex- 



traction of g and /5 differs at various orders. The results are given in Table |0|. At 
tree level only the ratio (3/g is determined, and the is rather large. We might 
next consider adding the contribution from the chiral loop corrections to D* Dj 
which go as y/rfi^. However, this does not lead to a consistent solution between the 
three data points unless jS is negative. This signals the importance of the Q'/ rric con- 
tribution in Eq. ( ^.23| ) corresponding to a nonzero heavy quark magnetic moment. 
Adding this contribution gives the results in the second row of Table p.l| , where there 
are now two solutions with similar in the region of interest. Adding the chiral 
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Order 


9 


/? 




9 


/9 




tree level 


^9 = 


3.6 


30. 








+Q' /rric + one-loop with y/fi\ 


0.23 


0.89 


4.3 


0.45 


2.8 


3.7 


+ cliiral logs 


0.25 


0.78 


4.1 


0.56 


3.2 


1.4 


one-loop with nonzero A,5,f • fc, 


0.25 


0.86 


3.9 


0.83 


6.0 


2.5 


without analytic rUq terms 














order irtq ~ 1/ irtc with 


0.265 


0.85 


3.0 


0.756 


4.9 


3.9 



Table 3.1: Solutions for g and /5(GeV^^) which minimize the associated with a fit 
to the three ratios in Eq. (|3.28| ). There are two solutions in the region of interest. 

logarithms, log (yu/mg), at scale /i = 1 GeV gives the solutions in the third row. 
Taking nonzero 5, A, and v ■ k m the non-analytic functions Fi and F2 gives the 



solutions in the fourth row of Table |3.1| , where the value of g in the second solution 
has increased by ~ 50%. For these two solutions only the analytic dependence 
has been neglected. Finally, the solutions in row five include the analytic depen- 
dence with the counterterms set to zero (at = 1 GeV). The uncertainty associated 
with these counterterms will be investigated below. It is interesting to note that the 
extracted value of g in the second column of Table ^.1| changes very little with the 
addition of the various corrections. 



One can see more clearly how these solutions are determined by looking at Fig. p75 
The central value for each ratio of decay rates in Eq. ( |3.28| ) gives a possible contour in 
the g-(3 plane, as shown by the solid (-D*'^), dashed (-D*+), and dotted (-D*) lines. An 
exact solution for two of the ratios occurs at the intersection of two of these contour 
lines. However, a good solution for all three ratios requires a point that is close to all 



three lines. The solutions in the fifth row of Table ^3] are indicated by stars in Fig. p~5 
The size of the experimental uncertainties can be seen in the 68% confidence level 
ellipses which are shown as shaded regions in the figure (for two degrees of freedom 
they correspond to ^ Xmm + 2-3). These regions are centered on the solid line 
since the D**^ ratio has the smallest experimental error. The errors in Eq. ( p.28| ) give 
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> 
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CO. 




0.4 0.6 



Figure 3.5: Solution contours in the g-[5 plane for the situation in row 5 of Table pTl 
The solid, dashed, and dotted lines correspond to solution lines for the D*^ , D*~^, and 
D* decay rate ratios respectively. The stars correspond to the minimal solutions 
and the shaded regions correspond to the 68% confidence level of experimental error in 
the fit. The hatched region is excluded by the experimental limit r{D*^) < 0.13 MeV 

m. 



the following one sigma errors on the two solutions 



9 
9 



0.2651;°?^ (3 = 0.85^;?^ GeV"^ , 
0.7561°^? /? = 4.90^:iGeV-^ 



(3.30) 



Both solutions fit the first two ratios in Eq. p.28| ), but do not do as well for the 
third. Minimizing the has biased against the third ratio as a result of its large 
experimental error. For this ratio the g = 0.265 and g = 0.76 solutions give values 
which are 4 and 13 times too small respectively. For the first solution it is possible to 
improve the fit to the third ratio with reasonably sized counterterms. For instance, 
simply taking = 2 gives B{D* Ds7c'^)/13{D* Ds^f) = 0.036. As we will see 
below, a large g solution with < 1 is only possible if g increases to ~ 0.9 and P 
increases to ~ 6.0GeV~^ (c.f. Fig. |3l^) . 
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The experimental limit r{D*^) < 0.13 MeV ^ translates into an upper bound on 
the value of g. Since B{D*^ D^l) is small, this bound is almost j3 independent 
and to a good approximation is 



c/ < 0.52 V VI + 3.01X-1 x = r(D*+)i'"^'V(0.13MeV). (3.31) 



For the situation in row five of Table ^]1] this excludes the hatched region in Fig. p75 
The limit on T{D*^) therefore eliminates the g ^ 0.76 solution at the two sigma level. 
Since this limit has not been confirmed by other groups it would be useful to have 
further experimental evidence that could exclude this solution. 

The central values in Eq. ( |3.30 ) have uncertainty associated with the parameter 



m^. Taking = 1.4 ± 0.1 GeV gives 0.25 < g < 0.28 and 0.79 GeV"^ < /3 < 
0.93 GeV-i for the first solution, and 0.72 <g < 0.80 and 4.6 GeV-^ < /5 < 5.3 GeV^^ 
for the second solution (in both cases the changes very little). There is also 
ambiguity in the solution in Eq. ( p.30|) due to the choice of scale /i (ie., the value 
of the counterterms ai, a^, ki and K5). Increasing /i to 1.3 GeV gives solutions 
{g = 0.28, P = 0.91 GeV^\ = 1-4) and {g = 0.78, P = 5.0 GeV~\ = 4.1), while 
decreasing /i to 0.7 GeV gives solutions {g = 0.25, /3 = 0.83GeV~^,x^ = 3.7) and 
{g = 0.72, (3 = 4.7 GeV~^, = 3.1). Note that the of the second solution remains 
large, while the of the first solution is reduced significantly by an increased scale. 

Another method of testing the effect of the unknown counterterms ai, a^, ki 
and K5 is to take their values at /i = 1 GeV to be randomly distributed within some 
reasonable range of values. We take —1 < ki, K5 < 1 and —2 < ai, < 2, with the 
motivation that the counterterms change the tree level value of Z° and by less 
than 30%, and give corrections that are not much bigger than those from the one-loop 
graphs. Near each of the two solutions 5000 values of g and (3 were then generated 
by minimizing the x^- This gives the distributions in Fig. |3.6| . The solution with 
g = 0.265 and /3 = 0.85 GeV^^ has fairly small uncertainty from the counterterms. 
The g = 0.76, f3 = 4.9GeV~^ solution has much larger uncertainty because the 



corresponding contour lines in Fig. 3.5 are almost parallel. For this solution the upper 
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Figure 3.6: Effect of tlie order counterterms ks, ai, and a^) on tlie solutions 
in Eq. ( p.30| ). The counterterms are taken to be randomly distributed with —1 < 
Ki,K5 < 1, —2 < «i,«5 < 2. For each set of counterterms g and (3 were determined 
at the new minimal x^- 5000 sets were generated near each of the two solutions. 



bounds are determined by the limits of a few MeV |57] on the D* widths. From this 
analysis we estimate the theoretical uncertainty of the solutions in Eq. ( |3.3U| ) to be 
roughly 



9 
9 



0.2651°:°^ (3 = 0.85+°:? GeV-^ 
0.76^°:? p = 4.9^^:? GeV"i , 



(3.32) 



at this order in chiral perturbation theory. The errors on g and /3 are positively 
correlated since the values of g and (3 are constrained in one direction by the small 
error on the D*^ rate ratio in Eq. ( p.28| ). 

From Eq. ( |3.31|) and Fig. |3.6| , we see that if the error in 



Dsn' 



(3.33) 



can be decreased by a factor of two, in conjunction with a limit of T{D*^) < 0.6 MeV 
then this could provide strong evidence that the g = solution is excluded. On the 
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other hand if the central values of the second and third ratios in Eq. ( |3.28| ) decrease, 
then a width measurement or stronger limit on r(Z}*+) will be needed to distinguish 
the two solutions. 

Using the extracted values of g and (3 gives the widths shown in Table |3.2| . The 
couplings were extracted at one-loop and order rriq ~ l/rric, so the predictions for the 
D* widths are made at this order. The experimental uncertainty in the D* widths 
is estimated by setting g and j3 to the extremal values in Eq. (|3.30|) , which gives the 
range shown in the second and fourth rows of the Table. The uncertainty from the 
unknown counterterms in the third and fifth rows is estimated in the same way using 
the uncertainties from Eq. ( p.32| ). Note that for the g = 0.265 solution the D* width 
is small due to a delicate cancellation in /is resulting from setting Z^^ xl/rric = l/rric- 
Keeping Z^^/nic to order rrig gives a D* width of 0.28 keV with a range of 0.1—0.4 keV 
for both the experimental and the counterterm uncertainties. 

Making use of HQS allows us to predict the width of the B* mesons from their 
dominant mode B* B'-f. Eq. (|3.23| ) gives the rate for B* B'-f with Q' = —1/3 
and rric mt,. Since the couplings f3i and /?2 are unknown these rates can not be 
determined at order 1 / m^fc, but we can include the order niq corrections. The B meson 
masses are taken from and we use mf, = 4.8 GeV ||107|| . We set 6 = 0.047 GeV 
and A = k ■ V = 0, but since the contribution Q' /nif, in Eq. ( p. 231) is numerically 
important it is kept in our estimate. For comparison the widths obtained with the 
g = 0.76 and (3 = 4.9 GeV~^ solution are also shown. 

As a final comment, we note that heavy meson chiral perturbation theory can also 
be used to examine excited D^*^ mesons, such as the p-wave states, Dq, Dl, Di, and 



D2 ||109| , |110| , |111| , |90| . To do so, explicit fields for these particles may be added to the 
Lagrangian giving a new effective theory. For interactions without external excited 
mesons (such as the ones considered here) these new particles can then contribute as 
virtual particles. However, since we have not included these heavier particles they 
are assumed to be 'integrated out', whereby such contributions are absorbed into the 
definitions of our couplings. 
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Widths (keV) 


D 


*o 


D*+ 






B* 


g = 0.265, 


1 


B 


26 
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9 = 0.76, 
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Table 3.2: Predicted widths in keV for the D* and B* mesons. The experimental and 
counterterm ranges are determined by the extremal values of g and jS in Eqs. ( p.30[ ) 
and ( |3.32| ). For g = 0.265 the D* width is small due to a delicate cancellation in fi^, 
as explained in the text. The uncertainty in the B* widths is large due to unknown 
^/fnc,b corrections. 



3.4 Summary 



For the D*^, D*^, and D*, the decays D* Dn and D* —>■ D'-f are well described by 
heavy meson chiral perturbation theory. Using the recent measurement of B{D*^ 
D'^'-f) the ratios of the Dj and Dtt^ branching fractions were used to extract the 
couplings g and (3. Here g and (3 are the D*D7i and D*D'~f couplings since order rriq 
and l/rriQ corrections have been absorbed into their definitions. Two solutions were 
found 



g = 0.265 /3 = 0.851;? GeV-^ 

g = O.retfati /3 = 4.91 j 1^7° GeV-^ (3.34) 

The first error here is the one sigma error associated with a minimized fit to the 
three experimental branching fraction ratios (see Fig. |3.5| ). The second error is our 
estimate of the uncertainty in the extraction due to four unknown counterterms ai, 
as, Ki and K5 that arise at order rriq (see Fig. p.6|) . 
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It is possible that the uncertainty from these counterterms can be reduced by 
determining them from other processes. For these corrections to contribute at low 
enough order in the chiral expansion we need processes with outgoing photons or 
pseudo-Goldstone bosons, such as semileptonic D decays to r], or vr. Here there are 
also SU(3) corrections to the left handed current which involve an unknown parameter 
riQ |8^. Information on Ki and k,[ can be determined from the pole part of the 



Ds Kiui form factor p2[. In a similar manner Dg rjiu^ can constrain Ki 
and K5, and a comparison of the form factors for K^iui and Dg rjiui 

gives information on k,[ and K5. These investigations were beyond the scope of this 
study. In principle, information about the constants ai, and 0:5 could be obtained 
from a measurement oi B ^ ^ii^e- The CLEO experimental bound on B iue 
{i = e, /i)| |112| | is roughly two orders of magnitude above the theoretical prediction, 
but due to the helicity suppression for B iug the branching ratio for B — > 'yiue 



may be up to an order of magnitude bigger |rT3| , |1 14 ] . 

Another possible approach would be to use large Nc scaling for the counterterms 
in 6Cg and 6C13. Terms that have two chiral traces are suppressed by a power of Nc 
compared to those with only one trace. In the large Nc limit the counterterms ki and 
«! would dominate, and K5 and 0:5 could be neglected, thus reducing the theoretical 
uncertainty. 

The smaller solution for g in Eq. ( p. 341 ) is fairly insensitive to the addition of 
the one-loop corrections (see Table p.l[ ). However, corrections at order rriq ~ l/rric, 
including the heavy meson mass splittings, were important in determining the solution 
with larger g. The limit r{D*'^) < 0.13 MeV [0] gives an upper bound on the coupling 
g (see Eq. (|3.31| ) and Fig. |3l5| ), and eliminates the g = 0.76, [3 = 4.9 GeV~^ solution. 
Experimental confirmation of this limit is therefore desirable. Note that the largest 
experimental uncertainty in our extraction comes from the measurement of B{D* 
DgTc^), and dominates the theoretical uncertainty due to decay via single photon 
exchange. A better measurement of B{D* Ds7r'^)/B{D* ^ Dgj) along with a 
limit r{D*^) < 0.6 MeV could provide further evidence that the g = 0.76 solution is 
excluded. However, if the central values of the second and third ratios in Eq. ( p. 28 ) 
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decrease then a width measurement or stronger hmit on r{D*^) will be needed to 
distinguish the two solutions. An improved measurement of B{D* —>■ Dg-n-^) may 
also give valuable information on the unknown couplings ki, K5, ai, and a^. 

The extraction of g has important consequences for other physical quantities [2- 
11]. For example^, for the B Tciui form factors with E^^ < 2 m^r, analyticity bounds 
combined with chiral perturbation theory give g Jb ^ MeV ||1 1 5|| . The solution 
g = 0.265 gives fs ^ 190 MeV for the B decay constant. However, for g = 0.76 
we have fs ^ 66 MeV, which is roughly a factor of three smaller than lattice QCD 
values, fs ^ 160 - 205 |l6], |ll|, [lT9i |l20|, |l2l|, [T22| . 



^ Glenn Boyd and Ben Grinstein, private communication. 
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Chapter 4 Vub from Exclusive 
Semileptonic B and D Decays 

The next generation of B decay experiments will test the flavor sector of the standard 
model at high precision. The basic approach is to determine the elements of the CKM 
matrix using different methods and then check for the consistency of these results. At 
the present time CP non-conservation has only been observed in kaon decay arising 
from — mixing. Many extensions of the minimal standard model (e.g., models 
with several Higgs doublets or low energy supersymmetry) have new particles with 
weak scale masses that contribute to flavor changing neutral current processes like 
— mixing, — B^ mixing, B —>■ K*^, etc., at a level comparable to the 
standard model. 

At the present time, the magnitude of the h ^ u CKM matrix element is deter- 
mined by comparing experimental results on the inclusive electron spectrum in the 
endpoint region with phenomenological models ||123|| , or by comparing experimental 
results on B ^ p i U£ and B n i ue with phenomenological models and lattice QCD 



results [ p.24|| . These two approaches yield remarkably consistent determinations of 



\Vub\, but have large uncertainties. 



In this chapter we discuss the proposal to determine \Vub\ [125, |126|| using a com- 



bination of heavy quark symmetry [20, ^ and SU{3) flavor symmetry. The basic 
idea is to compare D K* i vi with the Cabibbo suppressed decay D ^ piui. Using 
heavy quark symmetry the SU{3) violations in the form factors that occur in these 
decays are related to those that occur in a comparison of 5 — > K*i£ (or B K* ui i?£) 
with B —>■ piue. Therefore, experimental data on B —>■ K*£i in conjunction with 
data on D —>■ piue and D — > K*ii'£ can be used to determine \Vub\- This proposal 
is complementary to other approaches for determining \Vub\, since it relies on the 
standard model correctly describing the rare flavor changing neutral current process 
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B K*il 

In this chapter we compute corrections to these form factor relations which vio- 
late both chiral and heavy quark symmetry, and are non-analytic in the symmetry 
breaking parameters. We also reconsider the influence of long distance effects on the 
extraction of the B K* form factors from B —* K*i i. 

We denote the form factors relevant for semileptonic transitions between a pseu- 
doscalar meson P^'^\ containing a heavy quark Q, and a member of the lowest lying 
multiplet of vector mesons, V, by g^^^^\ /(^^^) and a^^~*^\ where 

{V{p\e)\q^,Q\H{p)) = ^g(^-^)e^,,^e*''{p + p')\p~pr, 
{V{p\e)\q^,^,Q\H{p)) = fi^-*y)e; + af-''\e*.p){p + p'), 

+ai^-''\e*-p)ip-p'),, (4.1) 

and e°i23 _ _5q-^23 = 1. We view the form factors g, f and a± as functions of the 
dimensionless variable y = v ■ v', where p = itlh v, p' = my v', and = {p — p'Y = 
m?^ + m\r — 2mH rriv y. (Although we are using the variable v ■ v', we are not treating 
the quarks in V as heavy.) The experimental values for the D K* I form factors 



assuming nearest pole dominance for the q dependences are ||127 | 

fi.D-.K^)^y^ a.9±0.1)GeV 



1 + 0.63 (y-1) ' 
_ (0.18 ±0.03) GeV"^ 
1 + 0.63 (?/- 1) 



= _ (0.49 ± 0.04) GeV- 

^ l + 0.96(i/-l) ^ ^ 



The shapes of these form factors are beginning to be probed experimentally ||127|| . The 
form factor a_ is not measured because its contribution to the D K* I decay am- 
phtude is suppressed by the lepton mass. The minimal value of y is unity (correspond- 
ing to the zero recoil point) and the maximum value of y is (m|) + m|^,)/ {2m£)mK*) — 
1.3 (corresponding to = 0). Note that f{y) changes by less than 20% over the whole 
kinematic range 1 < y < 1.3. In the following analysis we will extrapolate the mea- 
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sured form factors to the larger region 1 < y < 1.5. The full kinematic region for 
B pdvi'is I < y < 3.5. 

The differential decay rate for semileptonic B decay (neglecting the lepton mass, 
and not summing over the lepton type €) is 



dy 

Here S^^^^\y) is the function 

S^^-''\y) 



IT/ 12 

■™ m [yj 



48 7r3 



ruB 



(4.3) 



+ 4Re 
+ 4 



/(^-^)(y) (2 + y' - 6yr + 3r') 



(4.4) 



{y) m]jr\y'-lf 



+ 



^7(^^^)(y) m%r\l + r^-2yr){y^-l) 

/(^^^)(l/) \2 + y'- Qyr + ?>r') [1 + ^(^^^^(y)] 



with r = my /ran. The function ^^-^^^^ depends on the ratios of form factors 
^{H^v) I j(^H->v) g{H^v) I j(H^v) _ S^^^f'^y) can be estimated using combinations 
of SU (3) flavor symmetry and heavy quark symmetry. SU (3) symmetry implies that 
the B^ p+ form factors are equal to the B —>■ K* form factors and the B^ —>■ p° 
form factors are equal to 1/^/2 times the B — > K* form factors. Heavy quark sym- 
metry implies the relations ||125|| 
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(4.5) 



The second relation is obtained using a^^^^ "* = — a^^'^ \ valid in the large rric 
limit. 
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Using Eq. ( [4 .51 ) and SU (3) to get ^ £ U£ form factors (in the region 1 < 
y < 1.5) from those for D K*lv£ given in Eq. ( [4.2[ ) yields S^^^^\y) plotted in 
Fig. of Ref. ||126|| . The numerical values in Eq. (|4.2| ) differ slightly from those used 
in Ref. [ 126 |. This makes only a small difference in S^^^^\ but changes S^^^p^ more 
significantly. In the region I < y < 1.5, \S^^^^\y)\ defined in Eq. ( |4.4| ) is less than 
0.06, indicating that a^^''^ and g^^^P^ make a small contribution to the differential 
rate in this region. 

This prediction for S^^^p^ can be used to determine \Vub\ from a measurement of 
the B pii^e semileptonic decay rate in the region 1 < y < 1.5. This method is 
model independent, but cannot be expected to yield a very accurate value of \Vub\- 
Typical SU (3) violations are at the 10 — 20% level and one expects similar violations 
of heavy quark symmetry. 

Ref. ||126|| proposed a method for getting a value of S^^^p^ (y) with small theoretical 
uncertainty. They noted that the "Grinstein-type" pTl double ratio 



(y) 



(4.6) 



is unity in the limit of SU (3) symmetry or in the limit of heavy quark symmetry. 
Corrections to the prediction R{y) = 1 are suppressed by ms/rricfi {mu4 ^ 
instead of ms/AqcD or Aqcd/^^c,;;- Since Riy) is very close to unity, the relation 



S^^-p){y) = S^^-^*\y) 



fi^-P)(y) 



/ mB-rrip 
\mB — rriK* 



(4.7) 



together with measurements of |/*-^~*''^*''|, \f^^^'^^\, and S^^^^*^ will determine S^^^p"^ 
with small theoretical uncertainty. The last term on the right-hand-side makes 
Eq. ( ^4.7| ) equivalent to Eq. ([4.6| ) in the y —>■ 1 limit. The ratio of the (2 + — 
6yr + 3r^) [1 + 5^^^^\y)] terms makes only a small and almost y-independent con- 
tribution to in the range 1 < y < 1.5. Therefore, corrections to 
Eq. ( [4.7|) are at most a few percent larger than those to R{y) = 1. 

|y(D-^i4'*)| j^g^g oiYQg^dy been determined. \f^^~'^^\ may be obtainable in the fu- 
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ture, for example from experiments at B factories, where improvements in particle 
identification help reduce the background from the Cabibbo allowed decay. The mea- 



surement B{D p'^iue)/B{D K*'^iue) = 0.047 ± 0.013 pS] already suggests 
that |/(^~*/')//(-^~*^*)| is close to unity. Assuming SU{3) symmetry for the form fac- 
tors, but keeping the explicit my-dependence in S^^^^\y) and in the limits of the 
y integration, the measured form factors in Eq. ( [4.2|) imply B{D C,V()/B{D 
K*^Iv() = 0.044.Q 5'(^^'f^*) is obtainable from experimental data on _B ^ K* u^vi 
or B ^ K*i i. While the former process is very clean theoretically, it is very difficult 
experimentally. A more realistic goal is to use B K*ii, since CDF expects to 
observe 400 — 1100 events in the Tevatron run II (if the branching ratio is in the stan- 



dard model range) ||129|| . There are some uncertainties associated with long distance 
nonperturbative strong interaction physics in this extraction of S^^~^^*\y). To use 
the kinematic region 1 < y < 1.5, the form factor ratio f^^^p) j f^^'^^') in Eq. ( |4.7| ) 
must be extrapolated to a greater region than what can be probed experimentally. 
For this ratio, the uncertainty related to this extrapolation is likely to be small. 

The main purpose of this study is to examine the deviation of R from unity 
using chiral perturbation theory. We find that it is at the few percent level. The 
uncertainty from long distance physics in the extraction of S^^^^*^ is also reviewed. 
On average, in the region 1 < y < 1.5, this is probably less than a 10% effect on 
the B K*ii decay rate. Consequently, a determination of \Vub\ from experimental 
data on D ^ K*ii'£, D —>■ p iug, B ^ K*ii and B —>■ piue with an uncertainty 
from theory of about 10% is feasible. 



4.1 Chiral perturbation theory for the form factor 
ratio 

The leading deviation of R from unity can be calculated using a combination of 
heavy hadron chiral perturbation theory for the mesons containing a heavy quark 

^This prediction would be \Vcd/Vcs\'^ /2 ~ 0.026 with nip — ttik* ■ Phase space enhances D —t p 
compared to D ^ K* to yield the quoted prediction. 
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(section |2.2.1|) and for the lowest lying vector mesons (section p.2.2| ). We adopt the 
notations and conventions of Refs. |]79| , The weak current transforms as {3^, Iji), 
and at the zero recoil kinematic point there are two operators that are relevant for 
piQ) Y transition matrix elements (where P*^*-* = B, P^^^ = D, and V is one of 
the lowest lying vector mesons p,uj,K*,(f)). Demanding that the Zweig suppressed 
Dg ujIvi process vanishes relates the two operators, yielding ||8^ 



qal,{l - 75) Q = PT^mll.iX - 15)HPC, 



t 1 

ba\ ' 



(4.8) 



where Neb is given in Eq. ( |2.27| ). Here repeated SU{3) indices are summed and the 
trace is over Lorentz indices. contains the ground state heavy meson doublet, N 
is the nonet vector meson matrix ||58[, and /5 is a constant. The leading contribution to 
i?(l) — 1 arises from the Feynman diagrams in Fig. [4.1| . Diagrams with a virtual kaon 
cancel in the double ratio R. Neglecting the vector meson widths]^ these diagrams 
yield 



R{1)-1 



9 92 



12 7r2/2 



G(m^,A(^))-G'(m^,A(''))-G'(m^,A(^)) + G'(m^,A(")) , (4.9) 



where A*^^) = m^* — m^, A*^'^^ = m^* — ttid, and for m > A, 



G{m,A) 



7rm 



2A 



arctan 



y/m'^ — A^ 
A 



A Inm - 



(4.10) 



Here g2 is the pun coupling, g is the DD*tx coupling, and / ^ 131 MeV is the pion 



decay constant. In the nonrelativistic constituent quark model g = g2 = 1 fT^, while 
in the chiral quark model |5^ g = g2 = 0.75. Experimental data on t ^ uj tt Ur in 



the region of low 00 tt invariant mass gives g2 — 0.6 ||130|| . In chapter 3 we saw that 
the measured branching ratios for D* decays give g = 0.27 01 g = 0.76. 

For small A, Eq. ( |4.9| ) for i?(l) — 1 has a non-analytic y/fn^ dependence on the 
light quark masses. This cannot arise from corrections to the current in Eq. (|4.8|) or 



^The only significant width is that of the p meson. Since it occurs in the loop graph involving an 
77, neglecting the p width amounts to treating Tp/2m.fi ^ 1, which is a reasonable approximation. 
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p(Q) 




Figure 4.1: Feynman diagram that gives the leading contribution to -R(l) — 1. The 
dashed hne is a vr or an 77. The black square indicates insertion of the weak current. 

to the chiral Lagrangian, and must come from 1-loop diagrams involving the pseudo- 
Goldstone bosons vr, rj. Using the measured values of the pion and eta masses, 
Eqs. ( [4.9|) and (|4.1CI|) imply -R(l) = 1 — 0.035 (7 (72- There may be significant cor- 
rections from analytic terms of order ms/rric ~ 1/10 or from higher orders in chiral 
perturbation theory. However, the smallness of our result lends support to the expec- 
tation that -R(l) — 1 is very close to zero. There is no reason to expect any different 
conclusion over the kinematic range 1 < y < 1.5. 



4.2 Long distance effects and extracting S^^~ 



The decay rate for B — * K* u^Ui could determine S^^~^^*^ free of theoretical un- 
certainties. However, experimental study of this decay is very challenging. A more 
practical approach to extracting this quantity is to use B — > K*ii. The differential 
decay rate is 



dT{B K* 
dy 



a Y 2 



24 7r3 

x5(^-^*)(y)[l + %)]. 



\C9iy)\' 



\a '2 



[1 + A(l/)] 

(4.11) 



This and Eq. (|4.7| ) allow us to rewrite Eq. ( [4.3|) as 



dy 



\Vt*sVtb\^ a2 \c,(y)\^ + \C^o\^ 1 + A{y) 1 + div) m 



K* 



rriB - nip 
rriB - rriK* 



f(D-^'){y) 



dT{B K* 
dy 



(4.12) 
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which can be directly used to extract \Vub\- Unitarity of the CKM matrix imphes that 
|V^*\4fe| ~ l^^s^cfel with less than a 3% uncertainty. The fine structure constant, a = 
1/129, is evaluated at the ly-boson mass. d{y) parameterizes long distance effects, 
and will be discussed below. A{y) takes into account the contribution of the magnetic 
moment operator, O7 = (e/167r^) m^, (s^^ a^j, F'^'' (a factor of —4:GFV^*Vtb/V2 has 
been extracted out in the definition of operator coefficients). Ref. ||1 26|| (see also 



Ref. ||131|| ) found using heavy quark symmetry that A(y) —0.14 — 0.08(y — 1) in 
the region 1 < y < 1.5. Corrections to this are expected to be small since there are 
no l/rric corrections to A(l). Ciq is the Wilson coefficient of the operator Oio = 
(e^/lGvr^) (sl7^j6l)(£7'^75£). Cgi^y) takes into account the contribution of the four- 
quark operators, Oi — Oq, and the operator O9 = (e^/lGvr^) (sL'j^bL) (I'y^i). In 
perturbation theory using the next-to-leading logarithmic approximation 

C,{y) = C, + h{z,y)i3Ci + C2 + 3Cs + C^ + 3C, + Ce)-^h{0,y){C3 + 3C,) 
-^Kl, y) {AC, + AC, + 3^5 + C^e) + ^ (SCg + C4 + 3^5 + C^) , (4.13) 

where z = nic/mb- Here 



^/ N 81 ,8,4 2 ,^ , r- -\ In ^ + Vl^ - in ; x<l 

y 2 arctan(l/-\/a; — 1) ; a; > 1 , 

(4.14) 

where x = Au^rnl/ [m^g + m\, — 2m b ttik* y)- Using rrit = 175 GeV, rub = 4.8 GeV, 
rric = 1.4 GeV, as{mw) = 0.12, and as{mb) = 0.22, the numerical values of the 
Wilson coefficients are Ci = -0.26, C2 = 1.11, C, = 0.01, C, = -0.03, C5 = 0.008, 
Cg = -0.03, Cj = -0.32, C9 = 4.26, and Cw = -4.62. Of these, Cg and Cw are 
sensitive to mj (quadratically for rrit ^ rriw)- 

In Eq. (|4.13|) the second term on the right-hand-side, proportional to h{z, y) comes 



from charm quark loops. Since the kinematic region we are interested in is close to 
= Ami, ^ perturbative calculation of the cc loop cannot be trusted. Threshold 
effects which spoil local duality are important. It is these long distance effects that 
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give rise to the major theoretical uncertainty in the extraction of \Vub\ from the 
B £ differential decay rate using Eq. ( [4.12[ ).p| The influence of this long distance 

physics on the differential decay rate is parameterized by d{y) in Eq. ( [4.11|) , where 
setting d{y) = gives the perturbative result. 

For the part of the c c loop where the charm quarks are not far off-shell, a model for 
h{z, y) which sums over 1 cc resonances is more appropriate than the perturbative 
calculation. Consequently, we model the part of h{z, y) with explicit g^-dependence 
in Eq. ( ^4.14| ) with a sum over resonances ||134| , |135| , |136| , |137|| calculated using factor- 



ization 

The resonances V-^") have masses 3.097 GeV, 3.686 GeV, 3.770 GeV, 4.040 GeV, 
4.160 GeV, and 4.415 GeV, respectively, and their widths r^{„) and leptonic branching 
ratios B{-ip^"^ il) are known ||^. The factor k = 2.3 takes into account the 



deviation of the factorization model [|138|] parameter 02 from its perturbative value. 
Denoting the value of Cg{y) in this model by Cg{y), its influence on the differential 
decay rate is given by d{y) defined as 

my)\' + \Cio\' = i\C,{y)\' + \Cio\') [1 + d{y)] . (4.16) 

d{y) is plotted in Fig. [4.2| (solid curve). The physical interpretation of the 1 res- 
onances above 4 GeV is not completely clear. It might be more appropriate to treat 
them as DD resonances than as cc states. It is possible that for these resonances 
factorization as modeled by Eq. ( ^.15| ) with n = 2.3 is not a good approximation. 



Including only the first three 1 resonances in Eq. (|4.15| ), yields d{y) plotted with 
the dashed curve in Fig. [4.2| . 

This estimate of d{y) based on factorization and resonance saturation differs from 



that in Ref. [|126|] in two respects. Firstly, the phase of k is viewed as fixed because 



■^The four-quark operators involving light u, d, and s quarks also have uncertainty from long 
distance physics. However, this is expected to have a very small effect on the B K*£i rate. 
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Figure 4.2: d{y) defined in Eq. (^]T6|). The solid curve takes into account all six 1 
cc resonances according to Eq. ( ^.15| ), whereas the dashed curve is obtained including 
only the three lightest ones. 



recent data has determined the sign of the ratio of factorization model parameters, 
02/01, and the phase of ai is expected to be near its perturbative value ||139|| . Secondly, 
since the resonance saturation model only represents the cc loop for charm quarks 
that are not far off-shell, we have only used it for the part of h{z, y) in Eq. ( [4.14|) 
with explicit dependence, retaining the perturbative expression for the first two 
terms, — (8/9)ln2; + 8/27. The In^; term has dependence on m;,, which is clearly 
short distance in origin. This reduces somewhat the magnitude of d{y) and makes it 
more symmetric about zero (compare Fig. |4.2| with Fig. 6 of Ref. ||126|| ). It would be 
interesting to have a more physical separation between the long and short distance 
parts of the amplitude. 

Whether it is reasonable to use factorization for the resonances above 4 GeV can 
be tested experimentally, since these states cause a very distinctive pattern in dF/dy. 
In Fig. [4.3| the shape of dF/dy is plotted in the region 1 < y < 1.5 using the resonance 
saturation model for d{y) (solid curve). Experimental support for this shape would 
provide evidence that this model correctly describes the long distance physics param- 
eterized by d{y). Although d{y) gets as large as ±0.2, since it oscillates, its influence 
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y 

Figure 4.3: dT{B K*ei)/dy in units of [1^9(1)12 + |Ciop] m^* 
X I V^*Vt6p/(3847r^) as given in Eq. ( |4.11| ). The solid curve takes into account all 
six 1 cc resonances, the dashed curve includes only the three lightest ones, and the 
dotted curve is the perturbative result (i.e., d{y) = 0). 

on the B K*ii decay rate in the region 1 < ?/ < 1.5 is about —8% compared to 
the perturbative result (which is plotted with the dotted curve in Fig. |47^ ). Even 
if our estimates of this long distance physics based on factorization and resonance 
saturation has a 100% uncertainty (a prospect that we do not consider particularly 
unlikely), it will only cause about a 4% uncertainty in this determination of \Vub\- In- 
cluding only the first three 1 resonances in the sum in Eq. ( [4. 151) yields the dashed 
curve in Fig. In this case d{y) causes a —13% change in the B K*ii decay 
rate in the region 1 < y < 1.5. 



4.3 Nearer term prospects 

Without information on the y spectrum for the B decay rates in Eq. ( [4.12| ), it is still 
possible to determine \Vub\ by comparing the branching ratios for B —>■ piug and 
B — > K*ii in the region 1 < y < 1.5. Integrating Eq. ( [4.12|) over 1 < ?/ < 1.5 we can 
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write 



\V.,,X' 8^2 



y<1.5 



Here the barred quantities, Cg, A, and d denote the averages of \Cq{ii)\'^, A{y), and 
d(y) weighted with S^^^^*\y). Using the shape for S^^'^^*^ predicted from heavy 
quark symmetry, we find Cg = 4.58, A = —0.16, and d = —0.08. Note that the 
y-dependence of Cg is small and Cg is close to Cg. In Eq. (|4.171) the y-dependence 
of the ratio f^^'^^\y)/f''^^^*\y) has been neglected. If the shape of these form 
factors can be approximated with a pole form, then the pole masses of 2.56 GeV 
for fi^-*^*) and 2.45 GeV for f(^^P^ (corresponding to D** and to D**) imply that 
|/(i?-p)(^)//(s-^*)(^)|2 varies by less than 1% over the range 1 <y < 1.5. SU{3) 
symmetry and the measured D K* form factors imply that S^^'^p^ contributes 
only about 23% of the D pii^e decay rate. Using this prediction for 5^^~*''\ 
and assuming that f^^~*p^ and f(^^^*) have the same ?/-dependence, yields B{D 
pO iu,)/B{D ^ K*'> iu,) = 0.044 |/(^-^^)(l)//(^^^*)(l)|2. 

In the region = {pi+PiY < ""^j/^ (corresponding roughly to y > 2), one cannot 
use the double ratio and Eq. ([4.12|) . Moreover, the O7 contribution to the B K*ii 
rate is large and proportional to 1/g^, so the (leading order) heavy quark symmetry 
relations between the tensor and (axial-)vector form factor^ introduce a significant 
uncertainty. For < Tn^jj^., one can do better using SU{'i) flavor symmetry alone 
to predict dr(i? 7Tiui)/dq'^ from a measurement of dr(i? Kii)/dq'^. Since this 
region is far from g^^x' ^* P^le contribution |]72|, ^ |141| , [73( 1 is unlikely 
to upset the SU{3) relations. The O7 contribution to dr(i? — Kdtj/dq^ is at the 

10 — 15% level, fairly independent of q^. In the region (1 — 2)GeV^ < q^ < rn^jj^i 
^It was argued in Ref. |14C| that heavy quark symmetry can be used even at small q^. 
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neglecting m'j^ .j^/m'^^, 



(4.18) 



A similar relation also holds for integrated rates. 

A measurement of the B —>■ K*ii decay rate is unlikely before the Tevatron run 
II. Without this measurement, one has to rely on predicting the B ^ p form factors 
from D ^ p using heavy quark symmetry, or from D K* using both chiral and 
heavy quark symmetries. As discussed following Eq. ( |4.7] ), recent experimental data 
|128| | suggests that the SU{3) relation between fi^-^^*) and f^^^f^ is not violated 



by more than 15%. Furthermore, a prediction for B K*'y can be made using 



heavy quark symmetry and extrapolating the D K*iiy form factors | 142 |. The 
surprising agreement of this prediction with data may indicate that heavy quark 
symmetry violation in the form factor relations is smaller than anticipated. Heavy 
quark symmetry and the measured D K* form factors in Eq. ( [4.2|) imply that the 
B^ — > p'^iui branching ratio in the region 1 < y < 1.5 is 5.9 The measured 



decay rate ^ p+iu^) = (2.5 ± 0.4;°;^ ± 0.5) x 10"^ m together with \V, 



'ub\ 



0.003 imply that about 20% of B^ p'^ii^e decays are in the range 1 < y < 1.5. 

Despite the presence of long distance effects associated with the cc resonance 
region, the B —>■ K*ii rate can be used in Eq. ( [4.12|) to determine \Vub\ with a 
theoretical uncertainty that is about 10%. Experimental verification of the distinctive 
y-dependence of the differential rate associated with the 1 resonances above 4 GeV 
(see Fig. [4.3|) would reduce the theoretical uncertainty from long distance effects. A 
precise value of \ Vub\ may be available from other processes, e.g., the hadronic invariant 
mass spectrum in inclusive B — * X^ii^e decay [143, H^] or from lattice QCD results 



on exclusive form factors ||145|| before the B K*ii decay rate is measured. In that 



case, Eq. ( [4.12| ) gives an accurate standard model prediction for the B K*ii decay 



rate in the region 1 < y < 1.5. Comparison with data may signal new physics or 
provide stringent constraints on extensions of the standard model. 
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Chapter 5 NN Scattering 



5.1 Introduction 

Effective field theory is a useful tool for studying nuclear interactions. To describe 
low energy processes involving nucleons and pions in a model independent way, all 
possible operators consistent with the symmetries of QCD are included in an effective 
Lagrangian. A further advantage of effective field theory over potential models is 
that theoretical errors can be estimated in a systematic way. Contributions to an 
observable are organized by a power counting in Q/A, where Q is a momentum scale 
which characterizes the process under consideration, and A is the range of validity of 
the effective theory. A disadvantage of the effective field theory method is that the 
expansion parameter may not be very small, so that the description is not precise at 
low orders. 

In an effective field theory, ultraviolet divergences must be regulated and a renor- 
malization scheme defined. The ultraviolet divergences give a constraint on the power 
counting, because when a divergent loop graph occurs a contact operator that can 
absorb the divergence must be included at the same or lower order in Q. This is 



familiar from pion chiral perturbation theory as discussed in section |2]T|. The choice 
of regulator cannot affect physical results, but may make implementing a renormal- 
ization scheme easier. The renormalization scheme and power counting are also tied 
together. In a natural scheme, the renormalized coefficients of the operators in the 
Lagrangian are normal in size based on dimensional analysis with A. Once a power 
counting is established one can translate between different renormalization schemes at 



a given order in Q without changing the physical predictions. Recall from section 
that counting powers of Q/A in the nuclear effective theory is a subtle issue because 
of the large S-wave scattering length, a. In Refs. |^ Kaplan, Savage, and Wise 
(KSW) devised a power counting to take this into account. 
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Two different calculational techniques for the effective theory of nucleons are used 
in the hterature. In one approach, the power counting is apphed to regulated N- 



nucleon potentials and the Schroedinger equation is solved |2^, Eq, 146, 147, 148 



1491 , pO| , 115011 . Solving the Schroedinger equation is equivalent to including all ladder 



graphs with the potential as the two-particle irreducible kernel (see, for e.g., ||61|| ). 
The second approach, advocated by KSW, is like ordinary chiral perturbation theory 
in that the power counting is applied directly to the Feynman graphs which con- 
tribute to the amplitude. As discussed in section p.3.1| , a non-relativistic propagator 
is used which includes the kinetic energy term to regulate the infrared divergence at 
zero kinetic energy. In the Feynman diagram approach, dimensional regularization is 
the most convenient regulator, and analytic results are readily obtained. In the po- 
tential method, the Schroedinger equation is usually solved numerically. In practice, 
divergences are regulated and renormalized couplings are defined using a finite cutoff 



scheme. In Ref. |p.51|| , it has been explicitly shown that without pions the potential 
method can deal with large scattering lengths, and gives an expansion in Q/A. 

An important aspect of the KSW analysis is the use of a novel renormalization 
scheme, power divergence subtraction (PDS). In PDS, loop integrals in Feynman 
graphs are regulated using dimensional regularization, and poles in both d = 3 and 
d = 4 are subtracted. The subtraction of c? = 3 poles gives a power law dependence on 
the renormalization point, fi^, to the coefficients of four-nucleon operators. Choosing 
f^R ~ graphs with an arbitrary number of Co'*'*(/i/j) {s =^ Sq, or ^^i) vertices scale 
as 1/Q and must be summed to all orders. This is precisely the set of graphs that 
sums corrections that scale like (Qa)". Higher order contributions form a series in 
Q/A. In Ref. |Q, it is emphasized that it is possible to phrase the power counting in 
a scheme independent manner. The choice of scheme is simply to give natural sized 
coefficients which make the power counting manifest. PDS is one example of such a 



scheme. In Ref. [[152|| , it is shown how the KSW power counting can be implemented 
by solving the Schroedinger equation in a finite cutoff scheme. 

Pions can be added to the effective field theory by identifying them as the pseudo- 
Goldstone bosons of the spontaneously broken chiral symmetry of QCD. All operators 
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with the correct transformation properties are added to the effective Lagrangian. This 
includes operators with insertions of the hght quark mass matrix and derivatives, 
whose coefficients are needed to cancel ultraviolet divergences from loop graphs. In 
dimensional regularization, these ultraviolet divergences are of the form p^'^m^/e. 
For instance, for nucleons in the ^Si channel, the two loop graph with three pions 
and a two loop graph with two pions and one Cq have ultraviolet divergences of the 
form p^/e. This pole must be cancelled by a counterterm involving a four-nucleon 
operator with 2 derivatives. Because divergences of the form p^"/e must be cancelled 
by local counterterms, pion exchange can only be calculated in a model independent 
way if higher derivative contact interactions are included at the same order that 
these divergence occur P5|, |153|] . In Weinberg's power counting, pion exchange 
is included in the leading order potential. Therefore, graphs with arbitrary numbers 
of pions are leading order, while the counterterms necessary to cancel the ultraviolet 
divergences in these graphs are subleading. However, the potential method can still 
be used. As higher order derivative operators are added to the potential the accuracy 
is systematically improved, because the onset of the model dependence of the pion 
summation appears at higher order in Q/K. For example, the cutoff dependence of 
the two pion graph with one Co will be cancelled by cutoff dependence in C2. At a 
given order, the left over cutoff dependence in this method is a measure of the size of 
higher order corrections. 

Different estimates of the range, A^r, of an effective theory of nucleons with per- 
turbative pions exist in the literature. Some authors [154, 155 , |^ argue that A^ is 
as small as m^r, so that including perturbative pions is superfluous. One estimate 
of the range is given by KSW who conclude that A^ ~ 300 MeV. They point out 
that in PDS the renormalization group equation for the coefficient Co(/i_R) is modified 
by the inclusion of pions in such a way that for jiR > 300 MeV, Cq^hr) scales like 
/i^j instead of Since the power counting is no longer manifest above this scale. 



KSW conclude that the effective theory breaks down at this point. In Ref. ||155|| dif- 
ferent renormalized couplings are obtained. Here a breakdown of the power counting 
for C2{fiR) at fiji ~ m-TT is observed. A crucial question is whether a breakdown in 
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the running of the couphng constants is a physical effect or simply an artifact of 
the renormalization scheme. It is dangerous to draw conclusions based on the large 
momentum behavior of the coupling constants because the beta functions of the cou- 
plings are scheme dependent^. In section 5^, a momentum subtraction scheme is 
introduced where the power law dependence of the coupling constants persists even 
in the presence of pions, and for all values of /x/j > 1/a. This scheme is called the OS 
scheme, since in a relativistic theory it might be called an off-shell momentum sub- 
traction scheme. In Ref. ||6^, a similar scheme is applied to the spin singlet channel 
in the theory without pions, where it is shown to give results identical to the PDS 
scheme. The OS scheme is a natural scheme that works with arbitrary partial waves 
and with pions. Thus, the range of validity of the effective theory is not limited by 
the large fi^ behavior of the couplings. PDS is still a useful scheme in which to calcu- 
late observables. If one splits Cq{^r) into a non-perturbative and perturbative part, 
Co(/ii?) = Coif^Rj+Co^if^n)^ then CQ^ifiR) ~ l/fiR for all fiR > 1/a. We will see that 
this split is also necessary if we wish to avoid having the location of the pole in the 
amplitude shifted by chiral rriT^/A corrections. Once this split has been performed, it 
is straightforward to establish relations between the OS and PDS schemes order by 
order in perturbation theory, and any prediction for an observable will be identical 
in the two schemes up to the order in Q/A^^ to which it is calculated. Since in both 
schemes there is no scale where the power counting breaks down, it is possible that 

> 300 MeV. The importance of looking at results in several schemes is that it 
allows us to disentangle which results are physical and which are scheme dependent. 

Physically, one expects the effective theory to be valid up to a threshold where 
new degrees of freedom can be created on-shell. For elastic nucleon scattering, 
the relevant physical threshold is production of A resonances which occurs at p = 
^Mn{Ma - Mjv) = 525 MeV (the S-wave channels couple only to the A A inter- 
mediate state so p = ^/2Mn{Ma - Mn) = 740 MeV |5|]). Above this scale, the 
A must be included as an explicit degree of freedom. Below this scale, the A can 



^This is in contrast with dimensionless coupling constants like g in QCD. In that case the first 
two coefBcients of the beta function are scheme independent, so conclusions based on the behavior 
of the running coupling constant at small coupling (e.g., asymptotic freedom) are physical. 
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be integrated out leaving an effective theory of pions and nucleons. Rlio exchange 
becomes relevant at a scale, p ~ nip = 770 MeV. There is also a A^*(1440)A^ inter- 
mediate state with a threshold of p = 685 MeV. One might expect A^^ to be of order 
these thresholds. However, there is an intermediate scale of 300 MeV associated with 
short distance contributions from potential pion exchange^. Using dimensional anal- 
ysis, a graph with the exchange of n + 1 potential pions is suppressed by p/300 MeV 
relative to a graph with n potential pions. Comparison of the size of individual graphs 
is scheme dependent (for example the size of graphs differ in MS and in MS). The 
300 MeV scale applies only to a subset of graphs, and may change once all graphs at 
a given order in Q are included in the estimate. Therefore, 300 MeV can be taken as 
an order of magnitude estimate for the range of the theory, but the actual range may 
be enhanced or suppressed by an additional numerical factor. 

This then motivates the important question: How does one determine the range 
of the effective field theory? This is obviously a question of great practical impor- 
tance. Theoretical arguments can only give an approximate estimate for the range. A 
good example comes from SU(3) chiral perturbation theory. In this strong coupling 
theory, it is natural to expect that the range of the theory is the chiral symmetry 
breaking scale A^ ~ 2^/2U = 1200 MeV [IS^, |l58l. However, the convergence of 
the momentum expansion will depend on the particular process under consideration. 
For instance, in tt — tt scattering the range of the expansion is set by the threshold 
for p production, nip = 770 MeV. In this chapter, the range of the two nucleon effec- 
tive theory will be estimated using nucleon-nucleon scattering data. Our results are 



consistent with A.^ ~ 500 MeV. As we will explain in section ^.9| , the error analysis 
is applied to 6 rather than to pcot6 as in Ref. [|155|| . This range does not depend on 
the value of the renormalization point chosen, and is found in both the OS and PDS 
schemes. However, only next-to-leading order calculations have been used so it is hard 
to estimate the error in this value. When higher order corrections are computed, it 
should be possible to obtain a reasonably accurate estimate of the range of the two 

^ The phrase "potential pion exchange" wiU be used for a perturbative pion with energy inde- 
pendent propagator. This is sometimes called static pion exchange. 
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nucleon effective field tlieory witli perturbative pions. Tliis 500 MeV estimate is based 
solely on the phase shift data. The accuracy of predictions for deuteron observables 
11591 indicates Q/A ~ 1/3, which for Q ~ is A ~ 400 MeV. 

In section we review the power counting method of KSW 0, ^, and the PDS 
scheme. The importance of being able to count factors of the large nucleon mass in 
a non-relativistic effective field theory is discussed. We review the OS scheme, which 
is compatible with the KSW power counting. We describe the procedure for defining 
the renormalized couplings using local counterterms for each of these schemes. 



In section ^.3|, we discuss the theory with only nucleons, where A ~ m.,^. Local 
counterterms for both the PDS and OS schemes are computed. These counterterms 
are used to obtain the beta functions for the four-nucleon operators, and we explain 
why the beta functions for the most relevant operators in this theory are one-loop 



exact. In section f)A\ we explain how treating part of Cq perturbatively allows us to 
reproduce the effective range expansion with an amplitude that has its pole in the 



physical location at every order in perturbation theory. In section |5.5| we show how 
calculations with a cutoff regulator reproduce the dimensional regularization results. 

The theory with nucleons and pions is analyzed in section f).6\ In the "^Si channel, 
there are corrections to the PDS beta functions at all orders in Q. As examples, 
we compute the PDS beta functions for Cq to order Q, and for C2 to 

order In this channel, even in the limit 0, there are logarithmic divergences 

(poles of the form p^/e in dimensional regularization). In the OS scheme, the ^5*1 beta 
functions can be calculated exactly. We compute the exact beta functions for Co{fiR), 
C2{^J'r), and Ca^{^r) in the OS scheme in the ^Sq and ^5*1 channels. In section pTT] , 
the counterterms for the coupling constant D2{^ir) are derived in the OS and PDS 
schemes. 



In section p78| , we discuss why it is important to have /ir independent amplitudes 
order by order in the expansion. In the OS scheme amplitudes are /i^ independent, 
while in PDS jiR independent amplitudes can be obtained by treating part of Cq{iir) 
perturbatively. If this is not done then the sensitivity to is larger than one might 
expect [ [L54j| , for reasons we explain. Fits to the data are presented for different values 
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of fiji and the coupling constants in both OS and PDS are shown to evolve according 
to the renormalization group equations. 



In section |5.9| , an error analysis similar to a method due to Lepage is used to 
investigate the range of the effective field theory with perturbative pions at next-to- 
leading order. Weighted fits are performed for the scattering data in both the ^5*0 and 
^Si channels. Our results rule out ~ m^, and are consistent with ~ 500 MeV. 



5.2 Power counting and renormalization schemes 

In this section, the KSW power counting and compatible renormalization schemes are 
discussed. The theory containing only nucleon fields is considered first. The renormal- 
ized couplings are then defined in terms of local counterterms, and the KSW power 
counting for coefficients of four- nucleon operators is reviewed. Next, we consider the 
theory including pions. We review the power counting for potential pions, and ex- 
plain the origin of the 300 MeV scale associated with potential pion exchange. The 
PDS renormalization scheme is then discussed and we introduce the OS momentum 
subtraction scheme, which is also compatible with the power counting. 



Recall from section |2.3.1| that the Lagrangian in the two nucleon sector is given 

by: 

oo 

Cmn = Art ^,g^ + ^V(2M) + . . . ] iV - ^ c£ o£ . (5.1) 

s m=0 

The appearing in Eq. (|5.1| ) are bare parameters. To renormalize the theory, the 
bare coupling is separated into a renormalized coupling and counterterms as follows: 

oo 

/^bare /^finite ruv/^ /^finite / , \ \ ^ r"/^ /, \ /'K o^ 

n=0 

Note that we divide the counterterms into two classes. The first, which have the 
superscript uv, contain all genuine ultraviolet divergences. These include 1/e poles, 
if dimensional regularization is used, or powers and logarithms of the cutoff if a hard 
cutoff is used. We will also include some finite constants (e.g., the — 7 + ln(47r) that is 
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subtracted in MS) if this proves to be convenient for keeping expressions compact. By 
construction, these counterterms are fiR independent, but will depend on Cf^'**^. The 
renormalized coupling is denoted C2m(yUi?). The remaining counterterms, S"'C2m{fJ'R), 
contain no ultraviolet divergences and will be referred to as the finite counterterms. 
The choice of the finite counterterms differentiates between the schemes studied here. 
An infinite number of finite counterterms are needed because an infinite number of 
loop graphs are included at leading order. The renormalization is carried out order 
by order in the loop expansion. The superscript n indicates that 6"'C2m is included 
at tree level for a graph with n loops. When higher loop graphs are considered, the 
S^C2m counterterm takes the place of n loops . For example, at three loops we have 
three loop diagrams with renormalized couplings at the vertices, two loop diagrams 
with a 6^C counterterm, one loop diagrams with either one or two 5^C"s, and a 
tree level diagram with 6^C. Examples are given in section |5.3| . 

For the nucleon theory, the kinematic part of the power counting is very simple 



27| , p2| . Q is identified with a typical external momentum characterizing the process 
under consideration. For instance, in elastic nucleon-nucleon scattering Q ^ p, where 
p is the center of mass momentum^ Each nucleon propagator gives a Q^^, each 
spatial derivative a Q, each time derivative a Q^, and each loop integration a Q^. 

In the theory with only nucleons, the only graphs relevant to 2 — > 2 scattering 
are bubble chains. Consider a graph Q with L loops in the non-relativistic limit. In 
dimensional regularization, each loop will give a factor Mp/An, and there are L + 1 
vertices, each giving a factor — zCf^'^'^p^'". If the operator appears rim times in 
the graph (L + 1 = J2m ^rn) the result is: 

°° , /f /^finite 



^ = S n ( 4^'" ) P'^ ^here j = 5^2mn^ + L. (5.3) 

m=0 m=0 



For the scattering N{q + p) + N{q — p) ^ N{q + p') + N{q — p') it is useful to define 
p = \J MEtot — (p + ie, where Etot is the total incoming energy, and M is the nucleon mass. To 
simplify the notation we will work in the center of mass frame, q — 0, where p^ ~ p^ = p' ^ ~ AdE, 
and E is the center of mass energy. For external particles, one can always translate between E and 
p using the equations of motion. 



72 

If one matches onto the effective range expansion in MS one finds 



61 



(5.4) 



Note that all graphs Q are proportional to 1/M in agreement with the discussion in 
section p.3.1 . The large S-wave scattering lengths enhance the importance of some 
graphs compared to the p power counting. This affects the power counting for S-wave 
couplings, and through the mixing, couplings with L and/or L' = 2 and S = 1. For 
other channels we have the usual chiral power counting of p's. The power counting 
for insertions of four-nucleon operators is P2 



2m 



q{s, m) 



ig(s,m) 



where 



m - 1 for L = L' = 

m for 5 = 1 and (L, L' = 0,2), or (L, L' = 2, 0) 

m + 1 for S = 1 and L, L' = 2,2 

2m for all other S, L, and L' 



With the coefficients scaling as in Eq. (|5.5|), the graph Q scales as 



where i = g'(s, m) + L 



(5.5) 



(5.6) 



Note that the power of Q is less than or equal to the power of p, i < j. A useful 
mnemonic for this power counting is 1/a ~ Q, however, the power counting is still 
valid for Qa ^ 1. 

This Q power counting will be manifest in any renormalization scheme^ in which 
the C2m{fi'R) scale with /i^j ~ Q in such a way that Eq. (|5.5| ) is true. At leading order 
the counterterms 5"'C2m(/^R) will have the same Q scaling as the coefficient C2m(/^i?)- 
These schemes may differ by contributions in C2m(At_R) that scale with a larger power 
of fin/ A, since this will not change the power counting at low momentum. 

^Although the power counting can be implemented in different scliemes, the PDS scheme intro- 
duced in Ref. ||31| was very useful for initially working out the power counting. 
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Let us now discuss the theory with pions. To add pions, we identify them as 
the three pseudo-Goldstone bosons which arise from the breaking of chiral symmetry, 
SU{2)l X SU{2)ji SU{2)y. With the pions included in this way, we are doing an 
expansion in vHt^I and "pj K.,^. Note that in this theory, no matter how small p is 
made the expansion parameter will never be smaller than 111^^/ h.^^. This theory still 
includes the four-nucleon operators in Eq. (|2.3CI| ), but the short distance physics pa- 
rameterized by the coefficients is different because the pion is no longer integrated 
out. In the pion theory, the short distance coefficients should be independent 
of the scale m^r. All the m^r dependence is now contained explicitly in powers of the 
light quark mass matrix in the Lagrangian. 

Pions will be encoded in the representation, E = = exp (2in//), where 



n = „ ^ , (5.7) 




and / = 130 MeV is the pion decay constant. Under SU{2)l x SU(2)ji the fields 
transform as S ^ LLR\ ^ — ^ L^t/^ = U^K^, and N ^ UN. The chiral covariant 
derivative is L"^ = S'^ + K^^'^^t+^t^M^). with pions we have the following Lagrangian 
with terms involving 0, 1 and 2 nucleons: 

£^ = :^Tr(9^S9^St) + £^TrKE + m,St) 
o 4 

+ '-^N^am^^-C^^m + N^{^Do + ^)N (5.8) 

- J^^iioii - Di'^uTrim^m^Pt^NyiN^Pt^N) + .... 

s,Tn 

Here = \{^mqC, +^^"^g^^), f^q = diag(m„,m(i) is the quark mass matrix, = 
w{mu + mfi) where w is a constant, and = 1-25 is the nucleon axial- vector coupling. 
The ellipsis in Eq. (|5.8|) denote terms with more derivatives and more powers of m^. 



With pions there are additional complications to the power counting p5| , |153 | 



which are similar to those encountered in Non-Relativistic QED and QCD |T9l 22, 44 



The complications arise because there are two relevant energy scales for the pions. 
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~ Q'^/M for potential pions, and E^^ ^ Q for radiation pions. When the energy 
integral in loops is performed via contour integration, the graphs with potential pions 
come from terms in which one keeps the residue of a nucleon propagator pole. In these 
loops, the energy of the loop momentum is ~ Q'^/M and the energy dependent pieces 
of the pion propagator are suppressed by an additional Q'^ jNP' . Nucleon propagators 
give a Q^^ and the loop integrals give . There are also radiative pion graphs, in 
which the residue of the pion pole is kept. The power counting for these graphs is 
discussed in Chapter 6. In general graphs with radiation pions are higher order than 
those with potential pions. The combined propagator and vertices for a single pion 
exchange give so the pions can be treated perturbatively 0]. 

In general, pion exchange gives both long and short distance contributions. The 
short distance contributions from potential pions are important since they may limit 
the range of the effective field theory. A single potential pion exchange gives 



■ 9a Q-^^ap q- cr^s ^ 



2/2 g 2 + ^2 



^2 



'2/2 



q-aa(3 q- cr-ys m^^ q ■ a^p q ■ a^s 



q- 



g 2 (g 2 _|_ ^2 -) 



n ■ T2 



(5.9) 



where the spin indices connect to nucleon fields NlNpN^^Ns which belong to external 
lines or propagators. The first term dominates for (f ^ m2, and can be isolated 
by taking the limit m^r -^0. To study the dominant short distance contribution of 
pion exchange the m,r dependence can be neglected. Graphs with radiation pions are 
suppressed by powers 

of (m^/M) 1/2. In the non-relativistic limit, with only potential 
pions, the only loop diagrams are ladders. Consider an arbitrary graph Q with 
four point vertices, Cf^^*'', and k potential pions. For L loops, this graph has a total 



of L + 1 = k + Ylm ^rn vcrticcs, and with m^^ 







^ - O (if) ^ 



2 , fc °° 

i 



finite '\nm 
2m ) 



A-K (-iMg\ 



2 N k 



m=0 

oo 



M V 8 7r/2 



m=0 



Air 



(5.10) 

oo 

where j = 2 m rim + L . 



m=0 
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In the ^5*0 channel, the relation in Eq. ( |5.10| ) becomes an equality. The graph ^ ~ 



where i is given in Eq. ( |5.6| ). The power counting of the S^"'C2m counterterms is 
determined by the need to cancel ultraviolet divergences, and will not spoil the scaling 
for the renormalized coefficients, since i < j. For graphs with only potential pions 
{rim = 0), it appears that our expansion is in p/(300MeV) since 

T^ = |^-(300MeV)-^ (5.11) 
Ann ottj 

Comparing the size of potential pion graphs therefore predicts a range of 300 MeV, but 
the size of these graphs may change depending on the renormalization scheme (i.e.. 



the finite subtractions). As mentioned in section ^]T] the value of the loop graphs 
and contact interactions separately do not have unique values. It is not known a 
priori how the contact interactions will affect the range of the effective theory. The 
scale 300 MeV is therefore an approximate estimate for the range of the effective field 
theory with perturbative pions. A further discussion of this issue will be taken up in 
section IV. 

Next, consider the power counting for coefficients that multiply operators with 
powers of rriq. If we are interested in momenta of order m^, then one counts rriq ~ 

~ Q^. Therefore, any interaction term that has an operator with a total of 2m 
powers of p and m^^ will scale as Q^C^.m) ^j-^g^g q(^g^ ytl) is given in Eq. (|5.5| ). For exam- 
ple, ^"''m^ ~ . It is important to understand that in the KSW power counting 
D2 should be treated perturbatively even though the structure of the operator it 
multiplies is similar to that of the leading four nucleon operator with no derivatives. 
Graphs with radiation pions will also give contributions with powers of m^. 



5.2.1 The PDS scheme 

PDS is one scheme in which the KSW power counting is manifest. In PDS, we first let 
d = A and take the 5^^C2m counterterms to subtract 1/e poles as in MS. We use the 
notation yU/j for the renormalization point, and /i for the dimensional regularization 
parameter. In PDS, like in the MS scheme, one takes /i = ^r. In a momentum 
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subtraction scheme this is not necessary. The next step in PDS is to take d = 3 
and define the finite counterterms, S'^C2m{fJ'R), to subtract the l/{d — 3) poles in the 
amphtude. Graphs which contribute to S"'C2m{fJ'R) are those whose vertices have a 
total of 2m derivatives. When calculating the S"'C2m{f^R) we can take m,r = since, 
for instance, counterterms proportional to renormalize coefficients like D2{^r)- 
After making these subtractions everything is continued back to four dimensions. It 
is this second set of finite subtractions that gives the right power law dependence on 
^R. To define the coefficients that multiply operators with powers of m^, a similar 
procedure is followed except we count the powers of at the vertices. In PDS with 
just nucleons, all the graphs that affect the running of C2m{.l^R) are order Q^C^,™)^ 
except for those with intermediate states of different orbital angular momentum. For 
example, the beta function for C^"^^^ has contributions ~ Q {q{^Si,A) = 1), as well 
as contributions ~ from graphs with two C2^^^ vertices. When pions are 
included there are additional graphs that are sub-leading in the power counting and 
affect the running of the couplings. In fact, in section we will show that there will 

/3c \ 

be corrections to the PDS beta function for Cq ^ {(J'r) at all orders in Q. 
5.2.2 The OS scheme 

Another renormalization scheme that can be used to reproduce the power counting 
in Eq. ( |5.5| ) is a momentum subtraction scheme. A simple physical definition for 
the renormalized couplings can be made by relating the couplings to the amplitude 
evaluated at the unphysical momentum p = ifij^. This scheme will be called the OS 
scheme, since in a relativistic field theory this would be referred to as an off-shell 
momentum subtraction scheme. We start by dividing up the full amplitude as 

oo 

zA' = zJ2A^2n. + --- ■ (5.12) 

m=0 

Here contains the Feynman diagrams that will be used to define the coupling 
^2m(/^i?) (or equivalently the counterterms S"'C2m)- The ellipsis in Eq. ( p.l2|) denotes 



pieces that vanish as which are not needed to define C2m{l^R)- ^2m defined 
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to contain the remaining graphs that scale as Q'^^^'"^\ where q{s,m) is defined in 
Eq. ( [5 .51 ). The definition for the renormahzed couphng is then 

^■^U ^ = = -tC^i^if,^) . (5.13) 

rriTr = 

As we will see, this ensures that 6*2™ (/^h) scales in the desired way. In general, there 
may be divergent graphs scaling as Q * and p^™" {i < 2m) whose 1/e poles need to be 
absorbed by a S^^C2m counterterm. For example, consider the graph with two pions 
and one Cq shown in row four of Fig. This graph has a p"^/ e pole which is cancelled 
by a counterterm 5^^6*2. The finite part of this graph is used in Eq. (|5.13|) to define 
C'iifJ'R) because the graph is order Q. The key point is that since q{s,m) < 2m, 
an ultraviolet divergence that appears in a graph of a given order can always be 
absorbed into a coefficient that appeared at the same or lower order in the power 
counting. Therefore, we will define 5"^C2m in MS to subtract all four dimensional 1/e 
poles so that these subtractions are independent of the renormalization point. The 
finite counterterms are then fixed by the renormalization condition in Eq. (|5.13|) . 

In the OS scheme, the couplings CoifiR) and C2(yU_R) are defined by the renor- 
malization condition in Fig. |5]1|. This condition is to be imposed order by order in 
the loop expansion so that graphs with n loops determine 5"Co(/i_R). The = 
part of pion graphs contribute to 6*2™, (a^k) for m, > 1 in which case the condition 
rriTr = in Eq. ( p.l3| ) is important. In the theory with pions, we also need to define 
couplings multiplying powers of rrig, like D2 in Eq. (|5.8| ). To define these couplings 
we will not include all the terms in the amplitude proportional to m^. In particular, 
pion exchange graphs give long distance non-analytic contributions which will not be 
used to define the running of the short distance coupling D2{fiR). The idea that long 
distance physics must be excluded from the short distance coefficients is discussed in 
Ref. | ]155| . A detailed discussion of how we define D2{iir) in the OS scheme will be 
left to section V. 

Note that in the OS scheme there is another approach for calculating an amplitude 
in terms of renormalized couplings. One can calculate all loop graphs in in terms 
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Figure 5.1: Renormalization conditions for Cq^hr) and C2{fJ,R) in the OS scheme. 
iA^"^) is the four point function with Cq{^r) and 5"'Co(/iR) vertices, evaluated be- 
tween incoming and outgoing ^5*0 or ^5*1 states. The amphtude A^^^ contains graphs 
with one C2 or one potential pion dressed with Co bubbles. 

of the finite (or MS) parameters and then demand that the renormalization condition 
in Eq. (|5.13|) is satisfied. This gives expressions for the renormalized couplings in terms 
of the constants Cf^^**^. The amplitude can then be written in terms of renormalized 
couplings by inverting these equations. This simplifies higher order calculations. 

In the OS scheme, when an amplitude is written in terms of renormalized couplings 
it will be explicitly hr independent at each order in Q. The hr dependence in PDS 
with pions is cancelled by higher order terms. It is possible to obtain ^r independent 
amplitudes in PDS if part of Cq{^r) is treated perturbatively|^. Consequences of 
this /iij dependence will be discussed in section |5]^. In section ^]3| we will see that 
for the theory with just nucleons the OS scheme gives very similar definitions for the 



renormalized couplings to those in PDS. In section |5.6| , we investigate the running 
couplings in both schemes in the theory with pions. 
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2>0< X 

3XX 2>o< X 



Figure 5.2: One and two loop counterterms for Cq and C2. The solid lines are nucleon 
propagators, and symmetry factors are shown explicitly. The generalization to higher 
loops is straightforward. 

5.3 Theory with pions integrated out 

In this section, we compute the renormalized couplings in the non-relativistic nucleon 
effective theory without pions. We expect A ~ 171.,^. This theory will be examined 
in both PDS and the OS scheme. The renormalization program is implemented by 
explicitly calculating the local counterterms. In Ref. |Q, it is shown that the PDS 



and OS schemes give the same renormalized coupling constants in the ^Sq channel. 
Here we also consider the spin-triplet channel and higher derivative operators. Di- 
vergences in loop integrals are regulated using dimensional regularization. For the 
OS scheme, the same renormalization program can be carried out using a momentum 



cutoff regulator as shown in section p.5| . Following Ref. we will multiply each 
loop integral by {ix/2)^^~'^\ and define d = A — 2e. Since there are no logarithmic 
divergences in the nucleon theory, 5^^C2m = in dimensional regularization. 

In both PDS and OS scheme, it is straightforward to derive the finite counterterms, 
S"'C2m{fJ'R)- The tree level graphs with Co(/^/?) and C2(/ii?) satisfy the renormalization 
condition in Eq. (|5.13|) . Therefore, in both PDS and OS, S^Cq = 5°C2 = 0. At one 
and two loops we have the graphs in Fig. f).2[ In d dimensions, the two graphs in the 
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first row give 

(-.Co)^(^)r(3^) + ^ S'Co , (5.14) 

determining 6^Cq. In PDS, we define the counterterm to cancel the d = 3 pole in 
Eq. ( ^.14 ) and then continue back to four dimensions. In the OS scheme, we take 



d = A and demand that the contribution to the amplitude in Eq. ( |5.14| ) satisfies the 



condition in Fig. |5.1| . The counterterms calculated in each scheme are the same (with 
fi = fiji in PDS). In both schemes the counterterms determined from the graphs in 
Fig. p.2| are 



S'C2{fiR) = 2(^) C.MCoifiR) , = -3(^^y C2{fiR)Co{mf ■ 

Note that it is essential that loop graphs also have vertices with insertions of the 
counterterms. For instance, the contribution to the amplitude from all the graphs in 

is 



the second row of Fig. 



(5.16) 



If the one-loop graph with a S^Cq counterterm had been left out then the answer 
would have been proportional to (p^ + /i^) which is not correct. Since the loops in 
the nucleon theory factorize, the renormalized n-loop graph gives {ip + fin)^. Loop 
graphs will not always factorize once pions are included. 

It is straightforward to extend this calculation to n loops and to include higher 
derivatives. In both the OS and PDS schemes, this gives the following counterterms 
(s = %,3^i, n>l): 



^^0 



n+l 
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n+l f 



-ly 



(n + l)cf")M 



\n—l 



X 



47r 



(5.17) 



(■'•Si), 



47r 



<5"cf-^-^)w = (-1)-^ {^^y &:'^\,nr C?^-'^^\, 



<^"cr^^ w = (-1)"+' 



47r 



2 



Note that with fiR ^ Q, the counterterms have the same Q scahng as their corre- 
sponding couphng constant. In the PDS scheme, there arc also subleading terms that 
come from the mixing of angular momentum states. In PDS 



5"Cf^H//«) = (-1) 



47r 



(5.18) 



where the last term is suppressed by Q^. In the OS scheme 
5-cf'^\^n) = (-1)"+^ 



C^.^'^'ifiRr , (5.19) 



which is the same as the -S'o channel. In the OS scheme, graphs with two C2 
couplings and any number of Cg '^^^'s contribute to the beta function for Cg since 
they are order Q^. One might also ask about channels where the large scattering 
length does not effect the power counting. In this case C2m(A*i?) ^ Q^-> ^^"^ recover 
the usual chiral power counting. In our OS scheme, the counterterms (5"C2^(//ij) in 
these channels are either zero or a constant independent of /ir. 
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From Eq. ( |5.2|) one can derive the beta functions using 

P2m = /iiJT^ C2m{m) = I^RT; (5"C'2m (/ifl) • (5.20) 

The first few beta functions are 



'Si.'Di: (5.21) 

4''-' = (^)C■^''(..)^ 

/3f'-'°"=(^g^)cf"(^H)Cf'-"'-'(^„), 



in agreement with Refs. ^T], For S = states the beta functions are one loop 



exact in the sense that the contribution in Eq. ( ^.21|) comes from the one-loop graphs 



with the higher order graphs giving contributions which cancel. The reason for this 
cancellation is that the only loop corrections are in the bubble chain, and they form 
a geometric series. The sum of bubble graphs is just the chain of irreducible one loop 
bubbles for the full (point-like) propagator. An analogy would be QED, if the only 
possible graphs were the two point photon graphs with electron loops. In this case 
the beta function would also be one-loop exact because the graphs that are not IPI 
do not contribute. In general, the beta functions of higher order couplings may have 
contributions beyond one- loop in cases where angular momentum mixing is present. 

Expressions for the running coupling constants can be derived by summing the 
counterterms in Eq. ( |5.2| ) or by solving renormalization group equations. For s = ^So 



or '^5'i this gives 
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/^finite 



(47r) 



^finite 



finite \ 2 



1 _ MjlR 

Cjnitc 47r 



2 ' 



(5.22) 



where Cq"'**' and Cf'^'*'^ are constants which can be determined by specifying boundary 
conditions. Since the theory should be good for arbitrarily small momenta, one 
possibility is to demand that the amplitude reproduces the effective range expansion, 
pcot (5) = — 1/a + \rQp^ + 0{p'^). In Refs. |^ Q this matching was done at /x/j = 
giving Cq = Cf"'*'' = etc. We could equally well have chosen a different 

matching point (such as iir = 1/ a), and obtained the same results. For ~ Q, the 
running couplings in Eq. ( 5.22| ) have the scaling in Eq. ( |5.5|) . Written in terms of 
renormalized couplings the amplitude in the ^5*0 or ^Si channels is 



47r 



+ 



47r C2{^lR) 



MCo(pr) 



MCoM 



+ 0{Q) 



(5.23) 



and satisfies Eq. ( 5.13| ). The amplitude A is independent. It is interesting to note 
that we can choose a renormalization point where all loop corrections vanish giving 



m=0 



m=l 



An 



1 



An 
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H 2 , 

Ml/a + ip M\l/a + ipJ 2^ 



(5.24) 



The amplitude exactly reproduces the effective range expansion by construction. 
From Eq. ( |5.24| ) the range of the effective field theory can be estimated as A ~ 
2/ro ~ as expected. 



5.4 Reproducing the pole in the amphtude 

It is possible to choose the boundary condition for Cq{^r) to change the location 
of the pole that appears at each order in the expansion. For instance, consider the 
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following expansion of the amplitude in the theory without pions: 



A = — 



1 



M [-l/a + fp'^ + ...-ip 



-An 



+ 



An 
M 
+ A 



M [1/a + A + ip (l/a + A + ip) 



1 



-1/a - A + A + fp2 _^ _ 



(5.25) 



where A < 1/a. The series with A = and with A 7^ will both reproduce effective 
range theory, but differ in the location of the pole that appears at each order in the 
perturbative expansion. In the ^Si channel, the pole of the physical amplitude is at 



— = 7 



^/MEd = 45.7 MeV, 



(5.26) 



where Ed is the binding energy of the deuteron. For comparison, 1/a = 36.3 MeV 
in this channel. For A = 0, the pole that appears at each order in the perturbative 
expansion will be off by 30%. For some calculations, such as processes involving the 
deuteron [0, |160| , |161| , |162| , |163|| , a better behaved perturbation series is obtained by 
choosing 1/a + A = 7. The pole in the amplitude occurs at pcot 6{p) = ip so 



1 2 

A = 7 - 1/a = -ro7 . 



(5.27) 



Therefore, although the second term in Eq. ( ^.251) has a double pole, the residue of 
this unphysical double pole is zero. 

If we want to reproduce the expansion in Eq. ( |5.25| ) in the theory without pions 
then part of Cq^hr) must be treated perturbatively, Cq^iir) = Cq^^iir) + CQ^fin), 
where Cq^^hr) ~ l/Q and Cg(/i/j) ~ . Choosing the pole to be at p = ^7 gives 
^finite _ /^-j^l^^Mr^y Qg^gg amplitude becomes 



A' = — 



An 



1 ^An CliiiR 



1 ^An C2ifiR) p^ 



MlJ + ^p M(Co"^(/x^))M7 + *pP M(Co"^(/x^))M7 + ^p) 



(5.28) 



where the first term is order 1/Q, and the second and third terms are order Q^. This 



85 

is simply a reorganization of the perturbative series. The RGE's are 

l^n^C-'im) = ^C-'imf, (5.29) 
t'R^C'oim) = 2^C-'{^^n)Cl{^^n) + 0{Q). 

These can be derived by substituting Cq{hr) = Cq^{^r) + CI{hr) into the renormal- 
ization group equation for Cq{^b)- They can also be derived using the counterterm 
method described above. If we demand that the observed scattering length and ef- 
fective range are reproduced at this order then we find 

47r CI{^r) _1 47r_C^(/iR)_ ^ ro 

In order for the power counting of Cq{^r) to be consistent we must treat 7 — 1/a ~ Q^. 
From Eq. ( |5.27| ) we see that the first relation in Eq. (|5.3CI|) could have been derived 



by demanding that the residue of the double pole in the amplitude vanishes 



= 0. (5.31) 

jp=7 



At higher orders there will be order parts to Cq whose values are fixed by conditions 
analogous to Eq. ( |5.31| ). These conditions ensure that the position of the pole is not 
shifted by perturbative corrections. This is analogous to what we do in computing 
perturbative corrections to the electron mass in QED. 



5.5 Loop integrals with a momentum cutoff regu- 
lator 



Although the analysis in section |5.3| used dimensional regularization to regulate di- 
vergent loop integrals, the results for the coefficients C2m(/^i?) in our momentum 
subtraction scheme are independent of this choice. As an exercise we will derive the 
counterterms for Cq{hr) and C2{ijI'r) using a momentum cutoff regulator, A. This 
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will give us the chance to see what type of complications can arise using a different 
regulator. Note that this is not the same as using a finite cutoff scheme. There 
the momentum cutoff plays a double role as both a regulator and as part of the 
subtraction scheme. 



The graph in the first row first column of Fig. p.2| gives 



A + ^-ptanh-i(| 
2 A 3A2 



(5.32) 



An ultraviolet counterterm cancels the linear divergence, 



4tt \ n J 



(5.33) 



and the same finite counterterm, S^CoifiR) in Eq. ( ^.17|) is used to satisfy the condi- 



tion in Fig. |5.1| . The renormalized graph is then the same as calculated in dimensional 
regularization in section III. Note that contributions of order have been neglected 
in defining ColfiR) as required by our renormalization condition. An added complica- 
tion with a cutoff is that graphs with only Cq's give a contribution to the amplitude 
proportional to p^. However, as A — > oo, ptanh~^(p/A) 0, so these terms can be 
completely neglected. This will remain true even for higher loops since the countert- 
erms will always cancel dangerous powers of A that appear in the numerator. At n 
loops we find an ultraviolet counterterm of the form 



-M 
An 



\n+l 



(5.34) 



while the finite counterterms are given by Eq. ( |5 . 1 7| ) . 

The graph in the third row first column of Fig. gives 



2zCoC2- 



2 Jo {2n)W-p^ 



3 

' ' (5.35) 



Note that there are different contributions from this graph when the vertices are in 
the order C0C0C2 or C0C2C0. At order p^, this graph gives a correction to the 
counterterm S^'^'^Co, i.e., S^'^'^Cq ^ (J^'^^Cq + S^*'''^Co, where 

^i*,uv^^ = _^2Co^'^'**^C'f"^*"— . (5.36) 
47r Sir 

Unhke the contribution to 6^'^"Co in Eq. ( |5.34| ), 5^*'^^ Co is to be treated perturba- 
tively, so that at it only appears once in any graph. The justification of this fact is 
that this contribution to the counterterm appeared at order Q'^ (a purely formal trick 
to recover this counting is to take A ~ /i/j ~ Q). The counterterm 5^'"^C2 is fixed by 
considering the order terms in Fig. |5]^, row 3. From Eq. ( |5.35D (the tanh~^ piece 



can again be thrown away) we have 



1 M n n / 2A\ 



51."- C2 = — 2 C°>^''" . (5.37) 

An V ^ / 

The calculation for higher loops is similar and there are again corrections 5^*'^'^Co to 

/-/If \ " / 9A\""^ A3 
<5"*'-Co = (-—) n{n + l) (Co^"*^)" Cf"*"^ 



Att J ^ ■ ^ ^ ^ ^ ^ V y 37r ' 



The finite counterterms are the same as in Eq. (|5.17|) . Thus the running couplings 



and amplitudes with a cutoff are the same as found using dimensional regularization. 



5.6 Theory with nucleons and pions 

In this section, we study the renormalization of contact interactions in the effective 
field theory with pions. In the ^Si channel, graphs with two or more consecutive 
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Figure 5.3: Zero, one, and two-loop graphs with Cq and 5"Co vertices and potential 
pion exchange. The dashed lines denote potential pion propagators. 

potential pions do not factorize and give poles of the form p^™/e where d = A — 2e. 
We explicitly compute these poles for two loop pion graphs. There are also m'^/e 
poles in both the ^5*0 and ^5*1 channels at order [^, Because of these 1/e 
poles, pions cannot be summed to all orders in a model independent way. The finite 
counterterms in PDS and OS are different in this theory. Throughout this section we 
will take = 0, since we are only interested in the couplings C2m{fJ'R)- The D2{fiR) 
counterterms will be considered in section V. We compute the PDS counterterms and 
beta functions for Cq^hr) and C2(/i/j) to order Q. In PDS, Cq{hr) no longer obeys 
the Q scaling for > 300 MeV|^. This can be fixed by treating part of the coupling 
Co(/ii?) perturbatively as discussed in Section VI. The exact expressions for Cq^iir), 
C2{fJ'R), and Ci^^R) are given in the OS scheme and exhibit the correct Q scaling for 
all fiR > 1/a. Therefore, it is no longer apparent that the power counting breaks down 
at 300 MeV. The 300 MeV scale does appear in the short distance contribution to the 
amplitude from pion exchange, however, it can only be taken as an estimate for the 
range of the effective field theory once pion and contact interactions are both included. 
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Figure 5.4: The basic order Q graphs in the ^Si channel whose loop integrals do not 
factorize even for m.^ = 0. 

In section VI, we will discuss how experimental data suggests that At, > 300 MeV. 

To determine how the pions contribute to the beta functions for C2m(/ii?), we use 
the rules in section II. Some of the pion graphs that will be needed are shown in 



Fig. 5.3 



In both PDS and OS, the first step is to subtract 1/e poles with d = A — 2e. For 
two nucleons in the ^5*0 channel the spinor indices in Eq. ( ^.9|) are dotted into 6a5 
Therefore the m,r = piece of pion exchange reduces to a contact interaction and 
gives no 1/e poles. In the ^5*1 channel, graphs with two or more consecutive pions do 
not factorize and may have 1/e poles. Order Q graphs with two consecutive potential 
pions are shown in the first row of Fig. |^]^, and labeled a), 6), and c). We find 

4f4)V^). (5.39) 

1 14 , ^ , / VTU^ 

27 + 4 In (2) + 2 In ' ^ 



2 \2p) \ A'K 
.finite [9a V f-wM-^ 



2P 

-3i ((7fi°it°^2 f 9a\^ f—ipM^^^ 



2fJ \ 4n 



e 3 V — p2 _ is 

--37-61n2 + - + 31n^ ""^^ 



e 6 V —p^ — ie 



Graphs b) and c) have been written with Cq°'*'^ vertices to emphasize that the uv 
counterterm which cancels their divergent part is independent of /i/j. The divergence 
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in b) is cancelled by a tree level graph with the counterterm 



finite 



6C, 



— -7 + ln (7r) + 2-21n(2) 



(5.40) 



where the superscript 2 indicates that the counterterm comes in at two loops. The 
extra factor 2 — 2 In (2) is included because this leads to simpler analytic expressions. 
Expanding the Cq bubble graph (second row, first column of Fig. |5.3|) in e gives 



pM 



{Co? 



1 + e 



2-7-21n(2) + ln 



TT/i 



-p^ 



(5.41) 



When graphs with 1/e poles are dressed with Cq bubbles, the factors of [2 — 7 — 
21n2 + In (vr)] that appear are cancelled by similar factors from the counterterms. 
In fact, 5'^'^^C2 is the only uv counterterm we need for two potential pion exchange 
with = 0. The 1/e pole in c) is nonanalytic since it is proportional to p^. When 
graph c) is added to graphs d) and e) the poles cancel. These cancellations continue 
to occur when more Cq bubbles are added to h) and c). After including graphs with 
52.^-C2 we find 
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(5.42) 



- + 21n 
6 
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Note that for /i ~ p there are no large numerical factors from these graphs. 

In the ^Si channel, potential pion graphs without contact interactions also have 
p^ /e poles. The two loop graph with three potential pions (fourth row, second column 
in Fig. p.3|) is equal to 



(5.43) 
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- + . . . 

Le J 



In the Q power counting, this graph is order and will not be considered here. 
Because of these 1/e poles it is not possible to sum pion ladder graphs to all orders. 
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Now that the ultraviolet divergences have been removed from graphs b) and c), the 
finite subtractions can be performed. 



5.6.1 PDS 

For PDS in the ^Sq channel, we can compute the effect of potential pions on the 
C2m{f^R) counterterms to all orders in Q (neglecting relativistic corrections). For 
Co{^r), the relevant zero, one, and two loop graphs are shown in Fig. p.3[ The 
Co{fiji) and C2(/U/j) counterterms are 



°)(/i^) = (-1) 



V 47r 



5"C^r°^(/ii?) = (-l)"+'(ri + l) 



Air 
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2f 

Co(M + |^J C.iiin)- (5.44) 



The PDS counterterms in the ^Si channel will only be computed to order Q since 
the loop graphs with consecutive pions do not factorize. For this case it is essential 
to use the counterterms to carry out the PDS renormalization program. To define 
Co(/iR) at order Q, we set up the finite subtractions as in Fig. but leave out 
all graphs with more than two potential pions since they are 0{Q^) (we also neglect 
relativistic corrections that are order Q but come with an additional 1/M^). Note 
that in (i = 3 only the overall divergence (oc l/{d — 3)" for n loops) is needed since 
loops with counterterms will cancel the sub-divergences. Evaluating the graphs in 
Fig. |5.4j with = 3 and then continuing back to c? = 4 gives 
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(5.45) 



Using these values we find 



V Air 
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(5.46) 
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1, 



l)(n + 4)Co(M«r-^(|i 



2i 



for n > 2 



Note that for graphs with two consecutive potential pions, the fiR dependence does 
not come in the hnear combination fin+ip. For instance, adding the PDS counterterm 
to graph a) in Fig. ^ gives the hnear combination 3ip/2 + dfiR. 

In PDS, hke in MS, the renormahzed couphng C2{^ir) will depend on hi[ix\/ iiq) 
in such a way that the ln(/i|j) dependence in the amplitudes in Eq. (|5.39| ) is cancelled. 
Here /^o is an arbitrary scale expected to be of order A^r. At order Q we find 
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Note that the part of 5"C2(yU_R) proportional to In /ig) has a coefficient that sums 
up to at this order. From Eqs. (|CT| ), and ( |CT[ ) we find 
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(5.48) 
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Coim) + 1|] c^2(/iii) - 12 j Coim) [1 + m^^Coif^R) 



Note that in the ^5*0 channel all contributions to the beta functions beyond one- 



loop cancel, leaving them one-loop exact. In Ref. @], the last two terms in 
are absent, but should be included in the complete order Q calculation. Dimensional 
analysis implies that the ^5*1 beta functions can have corrections at all higher orders in 
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Q, since there is nothing to prevent the dimensionless factor {fijig\M)/{87if'^) ~ Q 



from appearing. In Ref. ||155|| , expressions for the beta functions are derived by 
demanding that dA/d^R = 0, but these are not the PDS beta functions. Since in 
all renormalization schemes dA/dfiR = 0, this condition is not sufficient to fix the 
renormalization scheme uniquely. 

When the beta function is not exactly known, the large /iij behavior is ambiguous. 
For example, the PDS beta function for Cq{hr) is [Q 
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47r 



2P 



0{Q) . (5.49) 



Two solutions which satisfy this equation to order are 
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(5.50) 



where fig = /i/j/Atvat and we have chosen Co(0) = Ana/M. The first solution is 
obtained by computing the counterterms 6^Co{fiR) to order and summing them. 
This solution falls as 1/fiR for all fiR > 1/a, and is numerically close to the qa ^ 
solution. The second solution is obtained by truncating and solving Eq. (|5.49| ). This 
solution approaches a constant as fiR ^ oo. The two solutions both solve the beta 
function to order but have very different large fiR behavior. Since the beta function 
for Cq ^^\fiR) can have corrections at any order in Q its large fiR behavior is unknown 
and not meaningful. This issue can be dealt with by expanding the C2m coupling 
constants in a series in Q and solving the beta function for each order separately. 



5.6.2 OS 

In the OS scheme, there is no such ambiguity since at a given order in Q the running 
of all the coupling constants that enter at that order are known exactly. The coupling 
CQ\fiR) has contributions only from the nucleon graphs discussed in section II and 
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therefore has the same beta function. For C2\iJ'r), the order (5° graphs in A2 include 
the nucleon graphs from section II, as well as the graphs with one potential pion and 
any number of Cq vertices. At tree level we add a finite counterterm to cancel the 
= part of the tree level pion interaction at p = ifiji 



9\ 



2/^ /^l 



(5.51) 



This counterterm is order Q"^ like C^i^y^n) itself. Since all the graphs in A2 factorize 
the higher loop counterterms are the same as in the theory without pions, so 5"C2 
for > 1 are given in Eq. ( |5.17| ). The exact beta function is then 
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Note that the finite ln(/i^/(— — ze)) terms in Eq. ( ^.39|) are order Q and in the OS 
scheme do not affect the running of C2{fiR), but rather Ci{^R). In terms of the finite 



constants Cq^^^*^ and Cf°'*'^ we have solutions 
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Although it may seem that the piece of C2\hr) that goes as will spoil the 

power counting for low momentum, in fact, the part dominates entirely until 

~ 1/a, since Cq^^^^ ~ a, Cg"'**^ ~ a?- Written in terms of renormalized couplings 
the = part of the next-to-leading order OS amplitude is 
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(5.54) 



which is order Q as desired. 

One might still ask if the problem with the 300 MeV scale will reappear in higher 
order coefficients. To check that this is not the case we compute the running of 
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the coupling Ci{^{i) in the OS scheme. The easiest way to compute this running 
couphng constant is to compute the order Q amphtude in terms of the finite couphngs, 
^2rn°i ^"^^ then demand that the amphtude satisfies the renormahzation condition 



in Eq. (|5.13| ). The graphs we need to compute include those with 



ii) 
in) 
iv) 



one C4 and any number of Cq's , 

two C2's and any number of Cq's , 

one C2, one potential pion and any number of Cq's , 

two potential pions and any number of Cq's . 



(5.55) 



Computing these graphs in terms of the finite couplings and then demanding that 
the amplitudes satisfy the renormahzation condition gives the OS couplings 
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(5.56) 
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, , M ^finite 



where here /i is an unknown scale expected to be of order A^. Again the pion contri- 
butions do not spoil the fiR scaling behavior, since they are suppressed by factors of 
the large scattering length. Note that at order Q the PDS coupling Ca^iir) 0] is the 
^ limit of Eq. i^^. 

In this section, expressions for the renormalized couplings Cq{iir), C2{fiR), and 
Ci^jjiR) were derived in the PDS and OS schemes working to order Q. For the ^Si 
channel, we have shown that Co{^r) has corrections at all orders in Q in PDS. Unlike 
PDS, the OS couplings C2m{fJ'R) can be computed exactly because they only have 
contributions at one order in Q. The OS couplings exhibit the correct fiR scaling for 
all fiR > 1/a. 
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5.7 The coupling D2{ij.r) 

In this section, the OS and PDS counterterms for D2{^r) are computed. To define 
-D2(/^ij) in the OS scheme, we take 



(5.57) 



where A^{D2) contains terms in the amphtude that are analytic in and propor- 
tional to m^. Only terms that are analytic in are kept because it is unnatural 
to put long-distance nonanalytic contributions that come from pion exchange into 



the definition of the short distance coupling ||155|] . For example, one potential pion 
exchange gives a m^/p^ ln(l + 4p^/m^) term. Including this in A^{D2) would give 
D2{f^B) both a branch cut at /i/j = 171^^/2 as well as explicit dependence on the scale 
m^. In the OS scheme, D2{^r) will be calculated as follows. First m1/e poles are 
subtracted. The finite counterterms are then determined by including graphs with 
a single D2{^ir) or potential pion and any number of Cq{^r) vertices in A^{D2)- 
Contributions from these graphs that are non-analytic in are dropped. 

There is a ^ra^/e pole in the 0{Q^) graph in the third row and third column of 
, so we have a counterterm 



Fig. 5.3 31 
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Note that when this counterterm is dressed with Co bubbles the extra factors of 
2 — In 2 from Eq. ( |5.41| ) will cancel without the need for an additional finite term in 
^2,uv^^_ After subtracting this counterterm the value of the two- loop graph is 
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For PDS we set fi = fiR and then find finite counterterms 
6'D2{fiR) = 2(^)Co(^«)Z^2(/i«), 



(5.60) 
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Here /xq is an unknown scale expected to be of order A^. 

In the OS scheme, the 5^D2{^ir) counterterm is the same as in PDS. In dimensional 
regularization logarithms of the form ^(/i^/m^) will appear in loop graphs. To make 
the /i^ dependent part analytic in we write 
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and then subtract the Inlfi^ / fij^) term with the counterterms. This will give -D2(/^i?,) 
a fiR dependence which cancels the ln(yU^/m2) in the amplitude. In the OS scheme, 
the in Eq. (|5.59|) gets subtracted along with the logarithm. We find 
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for n > 2. Summing the counterterms the solutions for D2{^r) are then 
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which can be written as 
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The parameter fi must be determined by fitting to data. With m^r ~ Q ~ P-r, 
D2{fJ'R)TTT''i ~ Q° in both OS and PDS, implying that D2{iir) should be treated 
perturbatively. 



5.8 Schemes and amplitudes 



In this section, the amplitudes in the ^5*0 and ^Si channels are presented to order 
both in PDS |31, § and OS. Fits to the ^Sq and "^Si phase shift data are done 



in both schemes for different values of ^r. As pointed out in section |5.4| , one has 
the freedom to split Cq{hr) into perturbative and nonperturbative pieces: Cq{hr) = 
Co^(/XH)+Co(/i_R), where Cq^^hr) ~ Q^^ and CI^Hr) ~ Q^. This division is necessary 
in PDS in order to obtain ^r independent amplitudes at each order. Furthermore, 
Cq^Ihr) ~ 1/ fiR so the coefficients scale in a manner consistent with the power 
counting for all fiR > 1/a. For convenience we will drop the superscript np in what 
follows. Some issues that arise in matching the pion theory onto the effective range 
expansion are also discussed. 

First, we give the nucleon-nucleon scattering amplitudes in the PDS and OS 
schemes. In PDS, the amplitudes were calculated to order in Refs. [^, ^. At this 
order, amplitudes in the ^5*0 and ^5*1 channels have the same functional form. 



(5.65) 



In both OS and PDS we have 
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The remaining part of the order PDS amphtude is 



1 2 



Coim? 2 2/2 W) 



+ C'oM 



(5.67) 



Note that since we have made a different finite subtraction than KSW the second 
term has a prefactor of 1/2, rather than a 1 as in Ref. @]. In the OS scheme 
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(5.68) 



The OS and PDS couphngs that appear at this order are related by a change of 
variables. Couplings on the left are in PDS while those on the right are in the OS 
scheme. 
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(5.69) 



In the PDS scheme, there are order Q° contributions to f3o (c.f., Eq. (|5.48| )). If the 
order contributions are separated from the order 1/Q pieces, the beta function for 
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(5.70) 



This equation has the solution 
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where K is a. constant which must be order for Cq{^r) ~ . (Recall, from 
Eq. ( |5.30| ) that ^^ = 7 — l/a< 1/ain the pure nucleon theory.) Including CI{iir) 
makes the PDS amplitudes explicitly /i/j independent. In performing fits to the data 
the constant K and coupling D2{^ir) cannot be determined independently. In what 
follows we will drop K when fitting and simply remember that the values of D2{^ir} 
extracted from the fits may differ from the renormalized coupling in the Lagrangian. 
In PDS, if CI{^r) is omitted from our expressions then D2{iir) does not follow the 
renormalization group equation, as we will see below. 

In the OS scheme, the constant g\/{2p) in Eqs.( ^.70l) and ( ^.711) is not present, 
so CI{^r) /Cq{iir)'^ is /iij independent. The OS scheme amplitudes A^~^'^ and A^^'^ are 
therefore hr independent without Cq^hr) as can be seen by examining Eqs. ( |5.53| ), 
and ( |5.6ij| ). In OS the constant K will also be absorbed into D2{^ir) for fitting. 

Using the Nijmegen phase shifts between 7 and 100 MeV, we fit the coeffi- 
cients Cq{^r), C2(yU_R) and D2{^ir). We took = 137 MeV. Clearly we would 
like to bias the fit towards the low momentum points since that is where the the- 
oretical error is smallest. This can be accomplished by assigning a percent error, 
^ p/(300MeV), to the data and then minimizing the function. In Tables O and 

we show the values^ of Co{^r), C2(/i_R) and D2{fiR) extracted from the fits for 
^R = 70, 100, 137, 160, 280 MeV. These values exhibit the fiR dependence predicted 
by the RGE's to ~ 1% in the ^Sq channel and ~ 4% in the ^Si channel. Since the 
amplitudes are explicitly fiR independent this deviation is a measure of the accuracy 
of the fitting routine. In Fig. |5.5| the results of the fits are shown. The results of the 
fits shown in the figure are identical in both schemes. Higher order corrections will 
give contributions to 6 of the form p^/A^. The error in 5 at p = 300 MeV is consistent 
with > 500 MeV. 

For processes involving the deuteron it is convenient to fix Cq{^r) using the 
deuteron binding energy, Cq^uIt^) = — 5.708 fm^. Fitting with this constraint we 

^ The coefficients extracted from our fits differ from those in Ref . Q because we have emphasized 
the low energy data as opposed to doing a global fit. It is interesting to note that using our PDS 
value C2(/ii? = f37MeV) = 11.5 fm'*, the prediction for the RMS charge radius of the deuteron 
becomes f.966fm which is within 1% of the experimental result. 
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Fit to % Fit to ^Si 







C2M 


D2M 




C2M 




70 


-6.48 


10.11 


-0.532 


-22.73 


171. 


-70.41 


100 


-4.71 


5.36 


1.763 


-9.93 


32.7 


-4.157 


137 


-3.53 


3.01 


2.000 


-5.88 


11.5 


1.500 


160 


-3.05 


2.25 


1.869 


-4.69 


7.32 


1.897 


280 


-1.79 


0.772 


1.105 


-2.19 


1.57 


1.004 



Table 5.1: ^5*0 and '^Si couplings in the PDS scheme. Co(/i_R) (in fm^), C2{hr) (in 
fm^), and D2{^r) (in fm^) are fit to the Nijmegen data at different values of (Xr. 

Fit to % Fit to 



/iii(MeV) 




C2{liR) 


D2{^^R) 


Cq^Hr) 


C2{^^R) 


D2{^^R) 


70 


-6.50 


9.75 


-6.047 


-24.1 


121. 


-170.1 


100 


-4.73 


5.33 


-1.143 


-10.0 


27.3 


-20.18 


137 


-3.54 


3.00 


0.378 


-5.92 


10.5 


-4.124 


160 


-3.06 


2.25 


0.658 


-4.74 


6.89 


-1.671 


280 


-1.80 


0.779 


0.692 


-2.23 


1.61 


0.2985 



Table 5.2: ^Sq and ^Si couplings in the OS scheme. Cq^jr) (in fm^), C2{hr) (in fm''), 
and D2{^r) (in fm"^) are fit to the Nijmegen data at different values of hr. 

find 6*2(^77) = 10.80 fm^ and D2{mj^) = 1.075 fm^ in the PDS scheme. The fit to the 
phase shift data with these values is as good as that in Fig. |5.5| . 

In PDS, it is necessary to break Cq{^r) into perturbative and non-perturbative 
parts to obtain amplitudes that are ^r independent order- by-order. If CI{iir) is 
omitted then the values of D2{^ir) determined from the fit will not follow the RGE. 
To see this we define the iir independent quantity 

and choose the constant c so that R = 1 for = 137 MeV. For other values of ^r the 
deviation of R from 1 gives the discrepancy between the values predicted by the RGE 
and those extracted from the fit. For ^r = 70, 280 MeV we find R = —0.53, 7.25 in the 
^5*0 channel and R = —0.52, 11.4 in the ^Si channel. These large deviations disappear 
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Figure 5.5: Fit to the phase shift data emphasizing the low momentum region. The 
sohd hne is the Nijmegen fit to the data 0, the long dashed line is the order l/Q 
result, and the short dashed line is the order result. 



if Co(yUi?,) is included. Without Co(/ii?), the PDS amplitude is fiR independent to 
the order that one is working. However, this residual fiR dependence gives larger 
corrections than expected ||154|| . The reason for this (as explained below) is that the 
residual /i^ dependence makes the tuning that was setup to give the large scattering 
length fiR dependent. 

Integrating out the pion gives low-energy theorems for the coefficients Vi in the 
effective range expansion [|], 



/ r\ ^0 2 4 fi S 

pcot[d) = -- + —p +V2P + V3P + V4P + . . . . 



(5.73) 



Performing a matching calculation between the two theories gives expressions for a, 
vq and the Vi in terms of the parameters in the pion theory. Since the theory with 
pions is an expansion in Q these predictions take the form of Taylor series in Q/A^ 
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- = 7 + E^' 







(5.74) 
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where Bn^ ~ Q*. At this time only the first coefficient in each series is known since 
pcot6 has only been calculated^ to order Q^. The notation 



7 



+ fJ^R 



(5.75) 



will be used to denote the location of the perturbative pole in the amplitudes. In 
PDS 



(2) 



(0) 



MP 



M ) Coim 



2 2 Knil) 



27 



-K. 



8nP 



+ 



2f 
ml 



(5.76) 



Note that if Cq{hr) had been neglected then B^^^ would not be /i/j independent. With 
fiR = rriT, Eq. ( |5.76| ) agrees with Ref. if their definition of D2{^r) is adopted. In 
the OS scheme we have 
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The value of the remaining Bn'^ determined at this order are the same in both schemes 
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(5.78) 



For n = 2,3,4, Eq. ( p.78[ ) gives the low-energy theorems derived in Ref. 0. 

The Vi in Eq. (|5.73| ) can be predicted in terms of one parameter, Cq^{iir), which 
is fixed in Ref. [0] by the condition An / [MCq"^ {^r)] + fiR = 1/a. Corrections to these 
predictions are expected to be 30 — 50% due to higher order Q/A terms. The 
extracted from the phase shift data |P, |165|| disagree with the low-energy theorems 



calculation of the low energy theorems at order in the ^Sq channel was presented at the 
INT Workshop on Nuclear Physics with Effective Field Theory, Seattle, WA, 25-26 Feb 1999 (S. 



Fleming, T. Mehen, and I.W. Stewart) |164]. This result will not be discussed here. 
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Figure 5.6: The effective field theory and Nijmegen Partial Wave analysis values 
of pcot5 are compared. The solid lines use pcot{S^^^ + 5^^^), the dashed lines use 
Eq. (|5.73| ) with the Vi from Ref . , and the dotted lines use the values of f , from the 
low-energy theorems. 

by factors of order 5. In Fig. |5]^, we see that the agreement of pcot(5'^'^^ + 6^^^) (solid 
lines )[] with the Nijmegen partial wave analysis is comparable to that of the effective 
range expansion with the Vi from the fits in Refs. |165|| (dashed lines). Note that 
our fit is more accurate at low momentum than the global fit in Ref. 0]. However, 
keeping only the first five terms from the low-energy theorems (dotted lines) gives 
larger disagreement at 70 MeV. This is not surprising since the pion introduces a 
cut ai p = i 171^^/2, so the radius of convergence of the series expansion of pcot (6) in 
Eq ( [5. 731 ) is ~ 70 MeV. At p = 70 MeV, one expects large corrections from the next 
term in the series. However, the fit values of Vi give good agreement with the data 
even at 70 MeV. It is possible that uncertainty from higher order terms in the Taylor 
series has been absorbed into V2, fs, and in the process of performing the fits. For 
this reason, the uncertainty in the values of Vi that were found from fitting to the 
data may be considerable. 

7 Note that when expanded in Q, pcotS ^ ip + 47r/[MA(-i)] - 47rA(o)/[A/(A(-i))2] + 0{Q^), 
which differs from p cot( (5^°^ + (5*^^^) by terms of o rder Q"^. The latter expression is used since the 



parameters in Table (5T) were fit using Eq. ( 5.80 ) 
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To get an idea of the error in f2, we will specialize to the ^5*1 channel. The 
Nijmegen phase shift analysis lists two data points for p < 70 MeV: p = 21.67 MeV 
where = 147.747±0.010 °, andp = 48.45 MeV, where 6^^^'^ = 118.178±0.021 °. 



Using a = 5.420 ± 0.001 fm and tq = 1.753 ± 0.002 fm ||165|| in the effective range 
expansion and fitting to the lowest momentum data point, we find V2 = —0.50 ± 
0.52 fm^, where the error in a, tq, and pcot6 have been added in quadrature. This 
differs by one sigma from both the value predicted by the low-energy theorem, f = 
— 0.95fm^, and the value from the fit, = 0.04 fm^. Since the range of the pure 
nucleon theory is 70 MeV, there will also be a ~ 0.1 fm^ error in this extraction from 
V3 and higher coefficients. This error was estimated by comparing the theoretical 
expression for p cot 6 with the first three terms in its series expansion. If we instead use 
the higher momentum point we find V2 = 0.03 ± 0.04 fm'^ with ~ 0.5 fm^ theoretical 
uncertainty. The uncertainty in these values of V2 is too large to make a definitive 
test of the low-energy theorems. 

Recall that the unnaturally large scattering length a is a fine tuning that was 
accounted for by demanding that in Eq. ( |5.75| ), Co(/i/j) is close to its ultraviolet fixed 



point, and 7 ~ 1/a. Examining the expression for 1/a in Eq. ( ^.74|) it may seem that 
this could be destroyed by chiral corrections. If D2{fiR) ~ Cq^hrY then the first term 
gives B^^ ~ 205 MeV. In fact from Table ^.3| , we see that the fit gives B^^ < 1/a. 
The reason for this small value is that since A^'^'> oc (A^^^^)^ the amplitude has a 
double pole. Since this pole is spurious (occurring from the perturbative expansion) 
the residue of the double pole must be small in order to fit the data. This leads to a 



good fit condition]^ which will be approximately satisfied 



0. (5.79) 



-?.p=7 



The condition in Eq. ( [5.79[ ) implies b'^'' ~ 47r7^/M. In fact this gives the right order 
of magnitude for the values of Bq determined from the fits in Table |5.3|. Similar 



good fit conditions occur at higher order keeping the coefficients Bq^ small. The 
division of Cq{^r) into nonperturbative and perturbative pieces is arbitrary, allowing 
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% Fit Fit 

/ifl(MeV) 7 l/g r-o 7 l/g rp 

70 -10.18 2.05 -8.124 0.01468 48.39 -15.82 32.57 0.01101 

137 -10.16 2.04 -8.121 0.01480 48.96 -16.76 32.19 0.01098 

280 -10.23 2.12 -8.105 0.01484 46.39 -12.64 33.76 0.01111 

Table 5.3: Values of 7, , l/o, and tq (in MeV) obtained from our fits. Three values 
of fiR are shown to emphasize that the value of the extracted parameters depends 
weakly on fiR. 

us to set up the theory so that the Q expansion for 1/a is well behaved. By imposing 
conditions like Eq. ( |5.79|) , chiral corrections to the location of the pole are absorbed 
into perturbative pieces of Coifin) order by order. Thus, we choose to spoil the short 
distance nature of Cq by giving it dependence ||166|| to keep the pole in the right 
place. 

In Table |5.3| we see that when B^^ is added to 7, values of 1/a are obtained which 
are close to the physical values, l/a(}So) = —8.32 MeV and l/a{^Si) = 36.4 MeV. 
It is encouraging that the value of 7 found from fits in the ^Si channel are close to 
the physical pole in the amplitude which corresponds to the deuteron, 7 = 45.7 MeV. 
Values for tq can also be predicted from the fits using Eq. (|5.77|) . Experimentally, 



roi^So) = 0.0139 Me V"^ and ro{^Si) = 0.00888 MeV"\ so the values in Table p| 
agree to the expected accuracy. 



'■TT 



5.9 Determining the range 

Here we will examine the phase shift data to see what it tells us about the range of 



the effective field theory with perturbative pions. In Ref. [|155|| , a Lepage analysis 
is performed on the observable p cot 6{p) in the ^Sq channel. Near 350 MeV the 
experimental ^Sq phase shift passes through zero. Therefore, the error Ipcot^^^^"^ — 
p cot 5^^""" I is greatly exaggerated since pcot6{p) -^00. To avoid this problem we 
will use the ^5*0 and ^Si phase shifts as our observables, since A6 = |(5NPWA _ ^EFT| 
remains finite for all p. The next-to-leading order amplitudes given in section V will 
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Figure 5.7: Error analysis for the phase shifts in the ^So and ^Si channels. A6 is the 
difference between the the effective field theory prediction and the Nijmegen partial 
wave analysisp|. The long and short dashed lines use the 0(<5°) and 0{Q) theoretical 
phase shifts respectively. 



be used. The phase shifts have an expansion of the form 5 = S'^^^ + 5^^^ + 0{Q'^/'^), 
where H 



-^In 
2 



271 



5(1) 



pM 

47r 1 + i^^(-i) 

ZTT 



(5.80) 



Recall that a momentum expansion of 6 would result in terms with only odd powers 
of p. However, the expansion for 6 in Eq. ( p.80| ) is not simply a momentum expansion, 
so the next-to-leading order calculation can have errors which scale as p^/A^. For 
example, once pions are included we can have a term tan~^(2p/m^) which is odd 
in p, order Q^, and scales as for large momenta. 

In Fig. |5.7| , we plot A6 versus p using log-log axes. Note that the sharp dips in 
Fig. |5.7| are just locations where the theory happens to agree with the data exactly. 
The Nijmegen data[0 is available up to p = 405 MeV. In a theory with just a 
momentum expansion the errors will appear as straight lines on the log-log plot as 
pointed out by Lepage pO[. In the pion theory the expansion is in both m^/A^ and 
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p/ At^i so this is no longer true. For p > we expect the errors to be of the formQ 
. (l,^, ...)£,..., (5.81) 



An Al J An \ An J Al 

TV n 1^ ^ n 

The fact that there is always a p/An error arises from the fact that, as seen in 
Eq. ( p.74| ), ro is reproduced in the effective field theory as an expansion in m^/A^. 
For p/An ^ vtinj An the slope of the lines on the plot should indicate the lowest power 
of p that has not been included. At low momentum the error in A^^"-^ is dominated 



by the pm^/A""*"^ term and the lines should be parallel. From Fig. 5/7 we see that the 
error is smallest at low momentum and increases as the momentum increases, which 
is how the theoretical error is expected to behave. 

It is clear that even for p ~ 400 MeV the next-to-leading order calculations are 
reducing the error in the phase shift. Because two new parameters are added at 
next-to-leading order it is always possible to force exact agreement at some value 
of p. However, if one were to force the data to agree too well at high momentum 
then this would destroy the agreement at low momentum. Since the improvement 
of the fit in Fig. |5.7] at high momentum does not come at the expense of the fit 
at low momentum this is evidence that the error is being reduced in a systematic 
way. At high momentum one expects that the error is ~ p^ / A^. From Fig. ^.7| , 
A(5 ~ 0.26 radians for p = 400 MeV, implying A^^ ~ 800 MeV. This is only a rough 
estimate for the range because we cannot yet exclude the possibility that the next-to- 
next-to- leading order phase shift has an anomalously small coefficient. Even though 



the lines in Fig. are not straight they should still cross at approximately the range 
of the theory since at this point higher order corrections do not improve the agreement 
with the data. This error analysis is consistent with the possibility that the range is 
> 500 MeV. 

Further information on the range of the effective field theory can be obtained by 



^ At momenta 1/a ^ p ^ nin we could have AS^^'' ^ B^^'' /p ^ rn^/Anp. However, as explained 

(2) 

in section VI, _Bq < 1/a so this term is very small. 
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examining electromagnetic processes involving the deuteron |]32| , |160| , |161| , |162| , |163 



such as the deuteron charge radius, electromagnetic form factors, deuteron polariz- 
ability, and deuteron Compton scattering. For these observables errors are typically 
~ 30 — 40% at leading order and ~ 10% at next-to-leading order. This is what one 
would expect if the expansion parameter m^/A^ ~ 1/3, implying ~ 410 MeV. 
This is consistent with our previous estimate for A^^. If the range is this large one 
should expect that the error in deuteron properties will be at the few percent level 
once next-to-next-to-leading order calculations are performedPI. 



5.10 Summary 

In this chapter the structure of the effective field theories of nucleons with and with- 
out pions is studied. We discuss a momentum subtraction scheme, the OS scheme, 
which obeys the KSW power counting. The method of local counterterms is used 
to obtain the renormalization group equations for the coupling constants in these 
theories. Using local counterterms defines the OS and PDS renormalization schemes 
unambiguously. Two-loop graphs with potential pions in the ^Si channel are com- 
puted and shown to have p'^/e poles. The presence of 1/e poles implies that the only 
model independent piece of pion exchange is the part that can be treated perturba- 
tively. We obtain the renormalized couplings Co(/i_R), C2(/i_R,) and Ci{iiR) at order Q 
in the OS and PDS schemes. 

We have emphasized why it is important to have independent amplitudes order 
by order in Q. Such amplitudes are obtained automatically in the OS scheme. In 
PDS hr independent amplitudes may be obtained by treating part of Cq^iir) pertur- 
batively. It is also necessary to treat part of Co(/U_r) perturbatively if we wish to keep 
the pole in the amplitude in a fixed location order by order in the chiral expansion. 
Another result concerns the large behavior of the couplings in this theory. In the 
OS scheme the couphng constants obey the KSW power counting for all fiR> 1/a. In 



^ Recently calculations of the deuteron quadrapole moment 1 167 and the ^5*0 phase shift |16J , 164| 
have been carried out at this order. 
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PDS the breakdown in the power counting for Cq{iir) is avoided if Cq{iir) is spht into 
non-perturbative and perturbative parts. Therefore, the breakdown of the scahng in 
the couphng constants is artificial. 

Next-to-leading order calculations of nucleon-nucleon phase shift data 0| provide 
fits to data at large momenta which are far more accurate than one would expect 
if the theory broke down completely at 300 MeV. Of course, this does not mean 
that nucleon effective theory can be applied at arbitrarily high energies. The scale, 
Mg\j [Sirf'^) ~ 300 MeV, is associated with short distance contributions from pion 
exchange and provides an order of magnitude estimate for the range. In the S-wave 
channel, A production and p exchange become relevant at ~ 700 MeV, which sets 
an upper limit on the range of the expansion. To get a better understanding of 
the range of the nucleon effective theory with perturbative pions one must examine 
experimental data. An error analysis of the S wave phase shifts with next-to-leading- 
order calculations seems to be consistent with a range of 500 MeV. Though next-to- 
next-leading order corrections need to be compared with data and other processes 
investigated, we remain cautiously optimistic that the range could be as large as 
500 MeV. 



Ill 



Chapter 6 Radiation and Soft Pions 

In the KSW power counting for nucleon-nucleon interactions, pions are included per- 
turbatively. In evaluating diagrams with pions, three types of contributions can be 
identified: potential, radiation, and soft. These pion effects differ in size and there- 
fore each have a different power counting. In this chapter the power counting for 
radiation and soft pions are discussed. The distinction between pion contributions 
arises because there are several scales associated with two nucleon systems. In this 
respect the theory is similar to NRQCD and NRQED |9|, |2 . 



In NRQCD there are three mass scales associated with non-relativistic systems 
containing two heavy quarks: the heavy quark mass M, momenta ~ Mf , and ki- 
netic energy ~ Mf^, where v is the relative velocity. QCD effects at the scale M 
are integrated out and appear as local operators in NRQCD. The remaining low en- 
ergy contributions can be divided into potential, radiation (sometimes referred to as 
ultra-soft), and soft pieces ||, |T6|, |4|, ||, |170], [171], [ITg. Potential gluons have 



energy of order Mf ^ and momentum of order Mf , radiation gluons have energy and 
momentum of order Mf ^, and soft gluons have energy and momentum of order Mv. 
The power counting for radiation gluons requires the use of a multipole expansion 
at a quark-gluon vertex ||169| , |4^ . The v power counting of potential and radiation 



gluons can be implemented in the effective Lagrangian by introducing separate gluon 



fields and rescaling the coordinates and fields by powers of f 1 45, 46]. In Ref. | 170 
the separation of scales was achieved on a diagram by diagram basis using a thresh- 
old expansion. The potential, radiation, and soft regimes were shown to correctly 
reproduce the low energy behavior of relativistic diagrams in a scalar field theory. In 
Ref. ||171| , |172|| it was pointed out that these effects may be reproduced by an effective 



Lagrangian in which separate fields are also introduced for the soft regime. Note 
that soft contributions come from a larger energy scale than potential and radiation 
effects. The heavy quark system does not have enough energy to radiate a soft gluon. 
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so they only appear in loops. In Ref. ||170|| it was shown that soft contributions to 
scattering do not appear until graphs with two or more gluons are considered. 

In the nucleon theory there is another scale because the pions are massive. For the 
purpose of power counting it is still useful to classify pion contributions as potential, 
radiation, or soft. For a pion with energy go momentum q, a potential pion 
has go ~ q^/M where M is the nucleon mass, while a radiation or soft pion has 
% ^ Q ^ iTT'TT- In a non-relativistic theory, integrals over loop energy are performed 
via contour integration. Potential pions come from contributions from the residue 
of a nucleon pole and give the dominant contribution to pion exchange between two 
nucleons. For these pions, the energy dependent part of the pion propagator is treated 
as a perturbation because the loop energy, go ~ g^/M ^ g. The residues of pion 
propagator poles give radiation or soft pion contributions. The power counting for 
soft and radiation pions differs. For instance, the coupling of radiation pions to 
nucleons involves a spatial multipole expansion, while the coupling to soft pions does 
not. 



6.1 Radiation pions 

In chiral perturbation theory the expansion is in powers of momenta and the pion 
mass m^. For power counting potential pions it is convenient to take the nucleon 
momentum p = Mv ~ [H, 0, so i; = itLt^/M ~ 0.15. The situation is different 
for radiation pions. There is a new scale associated with the threshold for pion 
production, which occurs at energy E = m^r in the center of mass frame. This 
corresponds to a nucleon momentum p = Qr, where Qr = y/Mrn^ = 360 MeV. 
Because the radiation pion fields cannot appear as on-shell degrees of freedom below 
the threshold E = 171^^, one expects that the radiation pion can be integrated out for 
p Qr- (Potential pions should be included for p > 171.^/2.) Another way to see 
that radiation pions require p ^ Qr is to note that in order to simultaneously satisfy 
kQ = k'^ + and ko ^ k ^ Mv"^ requires v ~ \Jm.,^/M ~ 0.38. 

The full theory with pions has operators in the Lagrangian with powers of 
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which give all the dependence. If the radiation pions are integrated out, then the 
chiral expansion is no longer manifest because there will be dependence hidden 
in the coefficients of operators in the Lagrangian. One is still justified in considering 
the same Lagrangian, but predictive power is lost since it is no longer clear that 
chiral symmetry relates operators with a different number of pion fields. Also, the m^r 
dependence induced by the radiation pions may affect the power counting of operators. 
For example, as shown in Ref. ||173| , p.74| , |175|| , integrating out the pion in the one- 



nucleon sector induces a nucleon electric polarizability a^; oc I/itLt^. Alternatively, 
one can keep chiral symmetry manifest by working with coefficients in the full theory 
and including radiation pion graphs. This is the approach we will adopt. 

The presence of the scale Qr modifies the power counting of the theory with 
radiation pions. In the KSW power counting, one begins by taking external momenta 
p ~ ~ Q. The theory is organized as an expansion in powers of Q. To estimate the 
size of a graph, loop 3-momenta are taken to be of order Q. However, potential loops 
within graphs with radiation pions can actually be dominated by three momenta of 
order Q.,.. To see how this comes about, consider as an example the graph shown in 
Fig. |6.1| c. Let q be the momentum running through the pion propagator, and let k 
be the loop momentum running through a nucleon bubble inside the radiation pion 
loop. The poles from the pion propagator are 



(6.1: 



ql-q'^-ml + ie (go - ^q^ + ml + ie){qo + ^q"^ + ml - ie) ' 



so the radiation pion has Iq'oI > ^tt- This energy also goes into the nucleon bubbles. 
The k integrand is largest when the nucleons are close to their mass shell. But since 
the energy going into the loop is ~ m^r, this occurs when k'^/M ~ m^r, i.e., k ~ Qr- 
We will begin by considering the contribution of radiation pions to elastic nucleon 
scattering at the threshold, E = m^. At this energy, external and potential loop 
momenta are of the same size and power counting is easiest. Because p ^ Qr it is 
obvious that we want to count powers of Qr rather than Q 

The power counting rules at the scale Qr are as follows. A scheme with manifest 
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power counting will be used, so that Cq{^r) ~ l/(M/iR), C2(/i/j) ~ l/(MAyu|,), etc., 
where A is the range of the theory. We will take p ~ /i/j ~ Qr- A radiation loop 
has qo ^ q ^ m^^ so (i^g ~ Q\jM'^^ where q is the momentum running through the 
pion propagator. A radiation pion propagator gives a M'^ jQ\^ while the derivative 
associated with a pion-nucleon vertex gives Ql/M. A nucleon propagator gives a 
M/Ql. External energies and momenta are kept in the nucleon propagator since 
E ~ p^/M ~ Ql/M. Furthermore, it is appropriate to use a multipole expansion 
for radiation pion-nucleon vertices which is similar to the treatment of radiation 



gluons in NRQCD ||169|| . Therefore, radiation pions will not transfer three-momenta 
to a nucleon. This is usually equivalent to expanding in powers of a loop momentum 
divided by M before doing the loop integral. The multipole expansion is an expandion 
in f = ^Jm-j^jM. A potential loop will typically have running through it either an 
external or radiation loop energy ~ Q'l.jM . Therefore, in these loops the loop energy 
^0 ~ Q'^./^^i while the loop three momentum k ^ Qr, so d'^k ~ Q^jM. It is not 
inconsistent for k ^ Qr while q ~ Ql/M, since three momenta are not conserved at 
the nucleon-radiation pion vertices. At the scale Qr potential pion propagators may 
still be treated in the same way, i/ (k^ — k"^ — ml) = — z/ (A;^ + m^) + O^kl/k"^), which 
has an expansion in Ql/M'^. We will see through explicit examples that this power 
counting correctly estimates the size of radiation pion graphs. 

Note that only the potential loop measure gives an odd power of Q^, so without 
potential loops the power counting reduces to power counting in powers of rriT,. The 
power counting here therefore correctly reproduces the usual chiral power counting 
used in the one nucleon sector [^, ^]. 

Graphs with one radiation pion and additional higher order contact interactions 
or potential pions are suppressed by factors of Qr/A relative to graphs with a single 
radiation pion and Cq vertices. The Qr expansion is a chiral expansion about = 0, 
so there is a limit of QCD where it is justified. The scale A is unknown. One possible 
estimate is Ann = 811 f"^ / {M g\) = 300 MeV since a graph with m+1 potential pions 
is suppressed by p/300MeV relative to a graph with m potential pions. However, 
this order of magnitude estimate only takes into consideration a partial subset of the 
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Figure 6.1: Leading order radiation pion graphs for NN scattering. The sohd hnes are 
nucleons, the dashed hnes are pions and 6M, 6Z are the mass and field renormahzation 
counterterms. The fiUed dot denotes the Cq[iir) bubble chain. There is a further field 
renormahzation contribution that is calculated in text, but not shown. 

graphs of the theory. As argued in chapter 5, it is possible that the range is of order 
the scale of short range interactions that are integrated out, implying A ~ 500 MeV. 
In fact, the accuracy of NLO computations of nucleon-nucleon phase shifts is in 
agreement with this physically motivated estimate of the range. We will assume in 
this section that an expansion in Qr/A is valid. This hypothesis will be tested further 
by seeing how well the effective theory makes predictions at p ~ 300 MeV. For 
example, processes with external pions could be considered. If the Qr/A expansion 
is not convergent, then application of the theory is restricted to p < Qr- 

The radiation pion graphs that give the leading order contribution to nucleon- 
nucleon scattering are shown in Fig. p.l\ . The filled dot denotes the leading order 
interaction between nucleons, a Cq{iir) bubble sum. We illustrate the power counting 
with an example, the graph in Fig. |6.1| d. For the moment, replace the Co bubble sums 
with single Cq vertices. Each Co gives a factor of 1/MQr and each nucleon line gives 
a factor M/Q^. The derivatives from the pion couplings combine with the radiation 
pion propagator to give a factor of unity. The radiation loop gives Q^/M'^, while the 
nucleon bubble loop gives Q^/M. There is also a factor of 1//^ from pion exchange. 
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and two factors of l/47r from the radiation loop giving a 1/A^. (A^ ~ 1 GeV is the 
chiral symmetry breaking scale.) Combining all factors, we find that this graph scales 
like Ql/{M^K^. This graph is suppressed relative to the leading order amplitude, 
by a factor of Ql/{M^K\) = ml/K\. Note that Cq bubbles are summed on 
external nucleon lines as well as in the interior of the radiation loop, and each graph 
in the sum has the same size. It is straightforward to verify that all graphs in Fig. |6.1| 
scale the same way. For external bubble sums we can simply use the vertex iA'^~^'> 
where A'^~^^ is the leading order S-wave amplitude, 

^'-" = -1^, (6.2) 

and the pole 7 = Att / MCq{^r) + ~ 1/a. Graphs with two radiation pions are 
suppressed by at least Ql/{M'^h!^ = m^/A^ and will not be considered. 



The first graphs we consider are those in Fig. |6.1| a,b. These graphs have contri- 
butions from potential and radiation pions, and it may not be obvious that a clean 
separation occurs. Here the energy integrals will be evaluated without any approxi- 
mations, after which the graphs split into radiation and potential parts. The graph 
in Fig. |6TT| a gives: 

,a{-i)9a_ f i i 

2P] (2vr)'^| + ^„_(£^ + ,e f-go-^ + ^e%^-g^-< + ^e• 
(6.3) 

(Throughout this chapter we will include a factor of (/i/2)^~'^ in the loop measures.) 
Performing the go integral gives a term from the residue of the nucleon pole and a 
term from the pion pole. 



2/2 j (27r)" (g - p)^ - ME ^2 + _ [| _ M!]2 

.2 r Jn„ 7^2 1 



(6.4) 



(6.5) 
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where n = d — 1. Eq. ( |6.4|) is the potential pion contribution. Expanding in — 
^W-? = [^^2^? gi^^s the result in Ref. [0, |. The subleading terms in this 
expansion are suppressed byQ m^/M^. Eq. (|6.5| ) is the radiation pion contribution. 
With |g| < M, we may take {q — p)'^/M — > p'^/M in the last two propagators, which is 
the same approximation that is made by performing the multipole expansion. Finally, 
we use the equations of motion to set E —p^/M = 0. It is important to note that we 
have not neglected E relative to For n = 3 — 2e, Eq. (|6.5|) becomes 



4/27 (27r)"(g 2 + ^2)3/2 (47r/)2Le 3 V 712 



, (6.6) 



where /i^ = ^'^'kc ""^ . Note that this integral is finite in three dimensions {n = 2) 



The next graph we consider is shown in Fig. |6.1[ b. We have chosen to route loop 
momenta so that q runs through the pion and ±k and ±{k + q) run through the 
nucleon lines. The momentum k is potential, while q can be potential or radiation. 
Doing the ko contour integral and combining the two terms gives: 

1)12^ f ^ f d^q q^ 1 1 



(27r)^go2_^-.2_^2+,, E-^^ + te E-^^ 

M M 



E 

A 



[E-^W-SI-^o + ^e][E-i^^-^ + qo + ^e] 

(6.7) 



Doing the go integral gives two terms, but the radiation and potential contributions 
are still mixed. Combining these gives 

■ ,^i-i).2 9A_ f d^kd^q g2 ^ ^ 

^ ^2/27 (27r)2n ^/^^^T^ E-f^ E-^^ 



1 



p _ (fc+g)^ _ fe^ _ 7^2 I ™2 



(6.8) 



^ Note that rn^/q^ ^ nriT^/AI, but wc have kept the term in the potential pion propagator 



in Eq. (6.4). We could consider expanding in m^/g using the asymptotic expansion techniques 
discussed in section B, but for p ~ TOtt these terms would have to be resummed. Unlike the soft and 
radiation contributions, there is no issue of double counting for potential pions, so for simplicity we 
will simply keep the in the propagator. 
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which can be spht into potential and radiation parts 
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(6.9) 
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The first term in Eq. ( p. 91 ) is the two- loop potential pion graph evaluated in Ref. [Q. 
The factors of {2k ■ q + q'^)/{2M) appearing in the denominators can be dropped 
because the loop integral is dominated hy k,q <ti M and therefore {2k-q+q'^) / (2M) ^ 



\/q'^^ m^. For the second term, which is the radiation pion contribution, this is 



equivalent to the multipole expansion. Momenta k ~ \fMrn^ and q ~ dominate 



the integrals in the second term. In Ref. ||1 70|| , the potential and radiation parts of 
the graph in Fig. |6.1| b were evaluated in the limit m^r = 0, and shown to correctly 
make up the corresponding part of the fully relativistic calculation. The calculation 
in Ref. [|17CI|| agrees with Eq. ( |6.9| ) for m^r = 0. Note that the radiation part would 
not agree if we assumed k ~ m,r and used static nucleon propagators in the radiation 
loop0 For = 3 — 2e the radiation part of Eq. (|6.9| ) is 
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where 



h(x) 



3 r(-|) 
2 r(|) 



F2 {-- 



1 1 



1 5. 



7^ 



^ '4'4^'^2'4^'''' 



(6.11) 



^ Furthermore, if static nucleon propagators are used one obtains a linear divergence requiring a 
non-analytic counterterm oc TOtt ]17(|. 
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For n = 2 the loop integral in Eq. ( |6.1CI|) is finite, except for the Ii term which 
has a jP'jin — 2) pole. In PDS this pole would effect the running of C^i^^i-R)^ but 
as we will see, contributions proportional to Ji will cancel between graphs. Since 
~ \l{MQr) the resuhs in Eqs. (|6]|, |6T0D are order Ql/{M^kl) as expected. 
At one- loop the 1/e pole in Eq. ( |6.6| ) is cancelled by a counterterm 



finite 9a 




- 7b + In (vr) 



(6.12) 



For higher loops the 1/e poles in Eq. ( |6.6|) , Eq. (|6.10| ), and 6^^'^"'D2 dressed with Cq 
bubbles cancel. Note that the 0(e) piece of the bubbles give a finite contribution. 
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(6.13) 



The result of combining Figs. |6.1|a,b and 5^^' D2 dressed by Cq bubbles is: 



a)+h) = 3i[A^-'^] 
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(6.14) 



Next we consider the graphs in Fig. |6.1| c,d,e. The loop integrals in these graphs 
vanish if the pion pole is not taken so there is no potential pion contribution. As 
pointed out in Ref. Q], emission of the radiation pion in these graphs changes the 
spin/isospin of the nucleon pair. Therefore, if the external nucleons are in a spin- 
triplet (singlet) state, then the coefficients appearing in the internal bubble sum are 

^°\f^R) (Cq The notation Co (Cq) will be used for vertices outside (inside) 

the radiation pion loop. We begin with Fig. |6.1| c. The contribution from the graph 
with m nucleon bubbles in the internal bubble sum is 
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where we used the multipole expansion and then the equations of motion to ehminate 
E and p from the first two propagators. All nucleon propagators have a go pole above 
the real axis, while the pion propagator has one pole above and one below. Therefore, 
the go contour is closed below. The d'^k integrals are also easily performed giving 



P 



-C^{^iR)MT{l-n/2) 



(47r)"/2 



■2 



(27r)" (g 2 + m2 )3/2 



(6.16) 



Note that the size of the loop momenta k in the nucleon bubbles is ~ y/Mrn^ even 



for p < y/Mm^. The fiR inside the brackets comes from inclusion of the PDS or 
OS S"'Cq{^r) counterterm graphs for the internal bubble sum. The integral will be 
dominated by g ~ m,r so the graph will scale as 

1 ml /v^M^x- 



Since ~ y/Mm^, all graphs in the sum are of order Qj?/(M^A 



For Figs. |6Tl|c,d,e the sum over bubbles should be done before the radiation loop 
integral. The reason is that an arbitrary term in the bubble sum has a much different 
dependence on the energy flowing through it than the sum itself. This can be seen 



in the g dependence in Eq. ( |6.16 ). If we integrate over g then terms in the sum may 



diverge whereas the integral of the complete sum is finite. In fact, for n = 3, Eq. ( |6.16|) 
has divergences of the form r(— 1 — m/4:)F{E'^ /ml) and r(— 1/2 — m/4)EF{E'^ /ml) 
where F is a hypergeometric function. These divergences are misleading because for 
momenta > 1/a we know that the correct form of the leading order four point function 
falls off as 1/p. Recall from section p.3.2| that for p ^ 1/ a the summation is 

infrared physics that we have built into our theory. For this reason the summation is 
performed before introducing counterterms to subtract divergences. (This approach 
is also taken in the analysis of three body interactions in Ref. |[1771 , |178|| ). Summing 



over m, Eq. ( |6.16|) becomes: 
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where 7' = 47r/MCQ(/iij) + fiR ~ 1/a. As expected the graph scales as Qf. In the 
hmit 77. ^ 3, /2 is finite and given by 
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I2 is manifestly fiR independent and is also finite as n ^ 2. 



Next we consider the graph in Fig. 6. Id. Integrals are done in the same manner 



as that of Fig. 3.1c. For n = 3 — 2e, Fig. 6. Id is 
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(6.20) 



Fig. |6.1| d is finite for n = 2. The 1/e pole in Eq. (|6.20|) is cancelled by a new tree 
level counterterm i5^'''^'^D2Tn1 where 5"^'^^D2 has the same form as Eq. (|6.12| ) except 
with a —12 instead of a —3. 



Evaluation of Fig. 3.I3 is also similar to Fig. 3.1c. For n = 3 — 2e we find: 
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This graph is finite for n = 2 except for the Ji term. A D2 counterterm cancels the 
divergence in the last line, 



5^-.i^D2 = 6 {C^ 
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(6.22) 



For two and higher loops the remaining 1 /e poles cancel between Eqs. ( p.2(]| , |6!2li|6.22|) 
and 5"^'^'^ dressed with Cq bubbles, so no new counterterms need to be intro- 
duced. The 0(e) piece of the bubbles again give a finite contribution. Combining 
Figs. p?T|c,d,e, and 5^^'^'^D2 and 5'^^'^^D2 dressed with Co bubbles gives 
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Fig. |6.1| f shows a two loop graph with a nucleon self energy on an internal line. It 
is important to also include graphs with the one-loop wavefunction and mass renor- 
malization counterterms, 6Z, 6M inserted on the internal nucleon line. We will use 
an on-shell renormalization scheme for defining these counterterms, which ensures 
that the mass, M, appearing in all expressions is the physical nucleon mass. The 
counterterms are: 
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The result from the graphs in Fig, 



. Olf is then 
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When Eq. ( |6.25| ) is added to Eqs. ( |6.14]J6.2^) the terms proportional to Ji cancel. 

To implement PDS we must consider the value of the graphs in Fig. |6.1[ f using 
Minimal Subtraction with n = 2. For = 2 + e we have 5M = 3(7^m^yu/(167r/^ e) 



and SZ = 0, which makes the sum of graphs in Fig. |6.1| f finite except for the Ji term. 
Finally, renormalization of the bare nucleon fields in the Lagrangian, 
Z = 1 + 6Z, induces a four- nucleon term 



6C 
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Since 6Z Q this term is treated perturbatively. A tree level counterterm 
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is introduced to cancel the 1/e pole. Dressing the operator in Eq. ( |6.26| ) with Co 
bubbles gives 
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Again, for n = 2 we have 5Z = so no new PDS counterterms were added. Note 
that if we had instead used bare nucleon fields then there would be no correction of 
the form in Eq. (|6.26|) . However, Eq. (|6.25|) would be modified because the last graph 
in Fig. |6.1[ f is no longer present. When this is combined with the contribution from 
the LSZ formula the sum of Eq. (|6.25|) and Eq. ( |6.28|) is reproduced. 

For PDS, the sum of graphs in Figs. |6.1| a-f are finite for n = 2 so no new finite 
subtractions were introduced. For n = 3, counterterms are introduced to renormalize 
the terms with in Eqs. (|6.14|J6.23|JO^ ) (in PDS = ^ir). In OS only terms 

analytic in are subtracted (including rri^ ln(/i^)). We find D2{ij.r) -D2(/Xi?) + 
AD2{fiR), with 
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Here k = 1/3 in PDS and k = in OS, and fiQ is an unknown scale. Note that the 
logarithm in Eq. ( |6.29| ) gives a contribution to the beta function for D2{f^R) of the 
form 



^"^'-Jfj^^^C^oW^ (6,30) 



This disagrees with the beta function of Ref. [^], because in that paper the beta 
function was calculated including only the one-loop graphs. 

Adding the contributions in Eqs. ( |6 . 1 4| , |6 . 23| , |6 . 2 5| , |08|) gives the total radiation 
pion contribution to the amplitude at order Q^: 



47r 



The first term here has the same dependence on the external momentum as an in- 
sertion of the D2 operator dressed by Co bubbles. Its dependence is cancelled 
by iiR dependence in AD2{fiR). Note that due to cancellations between graphs, this 
term is actually suppressed by a factor of '~f/Qr relative to what one expects from 
the power counting. The second term in Eq. (|6.31|) has a nontrivial dependence on 
E and is suppressed by an even smaller factor of These cancellations were 

not anticipated by the power counting. In the next chapter it will be shown that the 
dependence of Eq. (|6.31|) on (7 — 7') follows from the fact that the nucleon theory 



obeys Wigner's SU(4) symmetry in limit that 7,7' —>■ 0. Therefore terms at order 
will likely give the leading contribution of radiation pions to NN scattering. 

If we now consider momenta p ~ ^ Q Qr, we should fix fiR at the threshold, 
fiR = \/Mrn^, and expand in E/rriT, giving /2(-E/m^) = —3.9A+0{E/m^). Therefore, 
the dominant effect of the graphs that occur at order is indistinguishable from a 
shift in D2{fiR). Integrating out the radiation pions amounts to absorbing their effects 
into the effective D2 in the low energy theory. The result in Eq. (|6.31| ) is suppressed 
relative to the NLO contributions in Ref. by a factor of roughly 2(7— 7')/M ~ 1/10. 
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Since this is smaller than the expansion parameter, Q/A ~ 1/3, it can be neglected 
at NNLO. 

It is useful to consider more closely how the size of a calculation at p ~ Qr changes 
when we take p ~ m^r. At p ~ Qr the power counting gave contributions of size 
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where k is an integer. Taking p ~ m^r gives p/Qr = ^Jl<A~[rn^^ so terms with k < —\ 
will grow in size. It is the external bubble sums which are responsible for the factors of 
momentum in the denominator since A^"^^ ~ 1/p. The graphs in Fig. |6.1| have either 



zero, one, or two external bubble sums. Thus, there can be at most two factors of p in 
the denominator, and the biggest possible enhancement is by Mjm.,^. Therefore, the 
Q\ graphs are expected to have terms of order nij'^ ^ m^r, and mj"^ which agrees with 
what was found above. However, since the terms in Eqs. ( |6.14|) , (|6.23|) , and ( |6.25| ) 
that are proportional to l\ cancel, there is no term in the answer proportional to 

1/2 

vtltI . There is no obvious reason why this cancellation had to occur. It is interesting 
to note that the Ji terms scale as 1/-\/M, which is a higher power of M than the 
leading order amplitude. 

Recall, that in evaluating the graphs in Fig. |6.1| a mulitpole expansion was used, 
which in this case is an expansion in t> = a/^^/M. If the first correction in the 

1/2 

multipole series were to multiple a nonzero term of order 712-1/ then this would give 



an order m^r contribution. However, we have checked that for all the graphs in Fig. |6.1| 
the first correction in the multipole expansion gives a vanishing contribution. This 
occurs because these terms involve loop integrals with numerators that are odd in the 
loop momentum. 



At order graphs such as those in Fig. 3.2 will contribute to NN scattering. The 



graph in Fig. |6.2| a includes an insertion of the operator 

£ = iG2 [N^P.^'^N]^ [N^P^^aji^dj^^ - ^^dj^N] + h.c. . (6.33) 
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a) \ rx b) r\ c) 



G2 \ ^ ^ Cq 



Figure 6.2: Examples of order radiation pion graphs for NN scattering. 

(Note that because of the hermitian conjugate this operator is the same for s =^Sq 
and s =^5*1.) This graph will be dressed with Co bubbles inside and outside the 
radiation pion loop. The renormalization group equation for Gi gives G2{fiR) ~ 
l/(M/i|) ~ 1/(MQ2). When this is combined with the remaining factors of Qr, 
Fig. [6l^ a is of order Qf/{M'^A^) and is therefore suppressed by Qr/M relative to a 



graph in Fig. |0|. Power counting the graphs in Fig. |0[b,c gives Qj/{M^Af^]<fA 



2) 



yielding a factor of Qr/^NN relative to a graph in Fig. |6.1| . This provides an example 
of how graphs with potential pions seem to restrict the range of the effective field 
theory to Ann ~ 300 MeV. The 300 MeV scale applies only to a subset of graphs and 
may change once all graphs at this order are included. 

A complete NNLO calculation of the NN phase shifts at p ~ m^r requires a cal- 
culation of terms of order m^r. For graphs with one radiation pion at order Qf, an 
enhancement by 1/m^ could give an order contribution. Therefore, it would be 
interesting to learn if the same cancellation that occurred at order continues at 
higher orders. Although terms proportional to [A*^^^^]^ will appear at this order, 
they are factorizable. Therefore, the same cancellation that occured for the order 
graphs will occur in these terms. 



6.2 Soft pions 

In this section soft pion contributions will be discussed. We will see that there are 
graphs with non- vanishing soft contributions that should be included for p > m„. 
In soft loops, two scales appear: m^r and p = Mv. It will be shown below that it 
is necessary to take p ~ Qr when power counting graphs with soft loops in order 
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to avoid double counting. In other words v should have the same value as in the 
radiation pion calculation. A soft loop has energy and momentum go ~ ~ Qr? so 
d'^q ~ Q^. The mass of the soft pion is smaller than its momentum, and is treated 
perturbatively. Nucleon propagators in a soft loop are static (like in heavy quark 
effective theory, see Eq. ( p.21| )) since the loop energy is greater than the nucleon's 
kinetic energy ||171| , |172|| . Therefore, these propagators count as l/Qr- This power 



counting is identical to that proposed in Ref. @] except powers of Qr are counted 
rather than Q. 

Unlike potential pions, both soft and radiation pion pieces come from taking the 
pole in a pion propagator. Therefore, care must be taken not to double count when 
adding these contributions. This is accomplished by taking p ~ Q,. when evaluat- 
ing both soft and radiation pion graphs. This ensures that the soft and radiation 
modes have different momenta (~ Qr and ~ respectively). Integrals involving 
the scales Qr and can be separated using the method of asymptotic expansions 
and dimensional regularization [|170| , |179| , |180| , |181| , |182| , [L83|] . Consider splitting a 



loop integral into two regimes by introducing a momentum factorization scale L such 
that m,r < L < Qr- After the pion pole is taken in an energy integral over go, the 
remaining integral is of the form 

(Tq = I (radiation) + / (Tq (soft) , (6.34) 



L 



which is obviously independent of L. In Eq. ( |6.34| ) the power counting dictates that 



expansions in / Ql should be made so that each integral becomes a sum of integrals 
involving only one scale (m^ for radiation and Qr for soft). In dimensional regular- 
ization power divergences vanish, while logarithmic divergences show up as 1/e poles. 
Therefore, after expanding we can take L — > cxd in the radiation integral and L ^ 
in the soft integral. Taking the L oo and L limits may introduce ultraviolet 
divergences for the radiation integral and infrared divergences for the soft integral. 
When the radiation and soft contributions are added any superfluous 1/e poles will 
cancel. This will be illustrated with an explicit example below. The asymptotic ex- 



128 



a) b) c) d) 

► ■ ■ ► ► ■ ■ ► ► > ■ ► 



Figure 6.3: Examples of one-loop graphs which have soft pion contributions. Graphs 
a)-d) also have a radiation pion contribution, while in addition graph a) has a potential 
pion contribution. 

pansion procedure has been rigorously proven for Feynman graphs with large external 
Euclidean momenta and large masses ||184| , |185| , |186| , |187| , |188|| . It has also been shown 



to work for the non-relativistic threshold expansion of one and two-loop graphs ||170|| . 

Notice that it is crucial to expand the soft pion propagator in powers of m'^/Qf., 
because otherwise the radiation pion contribution may be double counted. As an 
example, consider the graph in Fig. 3J.a. Taking p ^ Qr implies Mv"^ ~ m.^. For 
the radiation pion contribution go ~ 9 ~ Mv"^ ~ rrij^, so we keep the in the 
denominator of Eq. ( |6.6|) . When computing the soft contribution, we assume go ~ 
g ~ Qr ^ f^n, and must expand the denominator in powers of rn^/Ql. The q 
integration is now scaleless so the soft contribution to Fig. |6.1| a vanishes in dimensional 
regularization. If we did not expand in nij^/Qr when evaluating the soft contribution, 
we would have double counted the radiation contribution. The same argument can 



be applied to all the diagrams in Fig. |6]1|. In each case the soft contribution vanishes. 

Examples of graphs which have a non-vanishing soft contribution are shown in 
Fig. |0| . These diagrams were calculated in Ref. [|189|| (although the S-wave channels 



were not analyzed there). In the KSW power counting they must be dressed on the 



outside with Co bubbles. (If Fig. |673| a or |673| b are dressed on the inside with Co bubbles 
then the soft contribution vanishes.) To see how these graphs obtain contributions 
from the soft and radiation regimes consider Fig. o^a. Unlike the massless case [170|, 



this graph has a radiation like contribution. In the ^5*0 channel the loop integral for 
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Fig. |0|a is 



-ig\\'^ f d^q i i q'^ 

W) J (2vr)^ I + go - (^f^ 1-50-^ + ^'^ q'o-q'-ml + ze 

^^^^^^ , (6.35) 

go - (g - ty - ml + ie 

where t = p' — p and ±p and ±p' are the incoming and outgoing nucleon momenta. 
Unhke the graphs in Fig. |6.1| , we are forced to route an external momentum, t, 
through a pion propagator. Taking a nucleon pole in Eq. ( |6.35p gives the potential 



pion contribution proportional to Mp. Taking the contribution from the pion poles 
gives soft and radiation contributions. Our power counting tells us that the leading 
order soft contribution will be ~ Q^, while the leading order radiation contribution 
will be ~ Qf./M^. In Eq. ( lOSi) the factors of E/2 - (g + p)V(2M) can be dropped. 
In the soft regime the factors of E/2 — {q + pY/{2M) are order Ql/M, and are 
dropped relative to go ~ Qr leaving static nucleon propagators. In the radiation 
regime E/2 — (g + p)^/(2M) after using the multipole expansion and equations 
of motion. This leaves 

2 \2f) J (27r)" g2 _ f)2 \ [^-2 + ^2]3/2 _ ^ ^2]3/2 



2V2/2y J (27r)« [g2 + m2]3/2 lg2_(^_f)2 + ,g g'2_(^_f)2_,g. 

(6.36) 

where n = d — 1. The singularity at g ^ = {q — tY is cancelled only in the sum of 
terms in the first line of Eq. ( |6.36| ). These terms can be calculated separately by 
introducing an ie in this denominator ||170|| , giving an average over ±ie as indicated^. 



The factor of {q~ty in the numerator can be removed by partial fractioning. For the 
soft contribution the scale of the loop momentum is set by the external momentum, 



The second line of Eq. ( |6.36 ) is more easily split into soft and radiation contributions. If we 



had used the integrand on the first line we would also have to consider {q — i)'^ ~ 
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Q'o ~ 1^ ~ 1^ ~ Qr- Expanding Eq. (|6.36|) in m^/q"^ gives 



9 a 



-2f 



cTq 

(2^ {2q-t 

1 



k1 



r(-i/2) 



1927r2V2/2 



t ± ze) ,^=0 

2 



r(-l/2 -m)r(m + 1) 



18 m: 



(6.37) 



where we have kept the first few terms in the expansion. The soft contribution 
starts at order Qf. as expected. The first 1/e pole in Eq. ( |6.37| ) is an ultraviolet 
divergence, while the second is an infrared divergence. For the radiation contribution 
Qo ^ lol ^ <^ Expanding in {2t ■ q}/t'^ gives 



2f 



9\ 



dJ'q 1 

(27r)" [g 2 + ^2]3/2 
2 



2q- 1 



r m=o =■= 



+ 



1927r2 V2/2 
- + In 



ml 



90 



72 

4 



- + In 



t2 



140^ + ... 



24 m: + 39' 



m: 



482 ml 
+ . . . 

3 ^ 



(6.38) 



The radiation contribution starts out as order Qj/M"^ as expected. Note that only 
powers of m^r = Q^/M appear. The 1/e poles in Eq. (|6.38| ) are ultraviolet divergences. 
When the soft and radiation contributions are added the infrared poles in Eq. ( |6.37| ) 
cancel a subset of the ultraviolet poles in Eq. (|6.38| ) . Adding Eq. (|6.37|) and Eq. ( |6.38|) 
we find 



1927r2 V2/2 



J.2 2 

t — ml 



^ + 181nr4U721n^^^ 



m^ 



(6.39) 
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where the remaining ultraviolet 1/e poles are cancelled by counterterms for C2 and 

If we are interested in making predictions for p ~ 777,7^, then powers of m^/t^ must 
be summed up. Summing the series in Eq. ( |6.39|) gives 



3a) 



1927r2 V2/ 



9 A 



90 m: 



1 



+ In 



58m: 



;i28m^ + IGmlt'^ - t^) / ^f^ + Ami - t 
t^Jf^ + Aml ^^^t^ + Ami + ^ 



(6.40) 



Since the coefficients in the series in Eqs. ( |6.37| , |6.38D diverge for p ' = Eq. ( |6.40| ) 
should be used when integrating over dt to obtain the ^5*0 partial wave amplitude 
(even for p ~ Qr). Eq. (|6.40|) agrees with the result of evaluating Eq. ( |6.36|) exactly. 
Although the asymptotic expansion is not necessary for evaluating Fig. p.3| a, it allows 
us to identify the radiation and soft contributions to this graph and verify that the 
power counting for each regime works. We also see that adding soft and radiation 
pions reproduces the correct answer without double counting. Recall that power 
counting with p ^ Qr was necessary to avoid double counting for graphs like those in 



Fig. 3.1 



For p ^ Qr the diagrams in Fig. |6]^ are order Q^/(/ A^), and are larger than the 
order graphs with a single radiation pion in Fig. |6.1| . Decreasing p to p ~ m.^ the 
graphs in Fig. |6.3| give contributions of the form 



P {ATlff 



F(j^/m: 



2 



(6.41) 



where F is a function. For p ~ m^r the graphs in Fig. |6.3| are order ml which is 



smaller than the graphs in Fig. |6]^ which include order m^^, mvr, and m'^'^ terms. 
It is interesting to note that the relative importance of these graphs changes with p. 



3/2 



The graphs in Fig. |6.3| dressed by Cq bubbles give a contribution that is the same size 
as a four nucleon operator with 6 derivatives, C^{^r)p^ ^ and are N^LO in the KSW 
power counting. Note that dressing these graphs with Co bubbles does not result in 
enhanced terms for p ~ m^r, unlike the graphs in the previous section. 
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It would be nice to see the expansion used in evaluating the radiation contribution 
to Fig. |6.3| a implemented at the level of the Lagrangian. It is not clear to us how to 
do this at the present time. In the radiation regime, the pion whose pole is taken 
can be thought of as a radiation pion. However, the other propagator gives factors 
of q- tjt^^ etc., which look more like insertions of non-local operators than the 
propagator of a field. Also, since in general p 7^ p ', the couplings for this second 
propagator change the nucleon momenta and therefore do not involve a multipole 
expansion. Finally, the result in Eq. ( |6.4CI| ) does not have an expansion in E/m.,^. 
So unlike the radiation pion contribution computed in section A, this contribution 
cannot be integrated out for p < y/MmT^. For these reasons, the use of the term 
radiation for this contribution differs somewhat from the usage in section A. 

To summarize, a power counting in factors of Qr = \fMm^ has been introduced 
which is appropriate for graphs with radiation pions. The order Q\ radiation con- 
tributions to NN scattering were computed and found to be suppressed by inverse 
powers of the scattering length. Soft pion contributions also have a power counting 
in Qt- For p ^ Qr they are ~ Q^, but are higher order than the radiation contri- 
butions for p ~ itLt^. Higher order corrections are suppressed by factors of Qr/-^, 
and whether or not this expansion is convergent is an open question. If the range of 
the two-nucleon effective field theory with perturbative pions is really 300 MeV, then 
contributions from radiation pions induce an incalculable error of order m^/A^ to the 
NN scattering amplitude in this theory. The validity of the Qr/A expansion can be 
tested by examining processes at p ^ 300 MeV such as those with external pions. 
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Chapter 7 Wigner's SU(4) Symmetry 
from Effective Field Theory 

In this chapter it is shown that in the hmit where the A^A^ ^5*0 and ^5*1 scattering 
lengths, a^^^°^ and a^^^^\ go to infinity, the leading terms in the effective field theory 
for strong A^A^ interactions are invariant under Wigner's SU(4) spin-isospin symmetry. 
This explains why the leading effects of radiation pions on the S-wave A^A^ scattering 
amplitudes, calculated in the chapter 6, vanish as a^^^'^^ and a^^^'^^ go to infinity. 
The implications of Wigner symmetry for A^A^ — > A^A^axion and 'yd^np are also 
considered. 



Wigner's SU(4) spin-isospin transformations are ||190| , |191 



SN = iaf,^a^'T'' N , N={ ^ ). (7.1) 



n 



In Eq. ( [7.1|) , = (1, <?), = (1, r), and a^^, are infinitesimal group parameters (the 
notation here is that greek indices run over {0, 1, 2, 3}, while roman indices run over 
{1,2,3}). The a matrices act on the spin degrees of freedom, and the r matrices 
act on the isospin degrees of freedom. (Actually the transformations in Eq. (|7.1| ) 
correspond to the group SU(4) xU(l). The additional U(l) is baryon number and 
corresponds to the aoo term.) 

Consider first the effective field theory for nucleon strong interactions with the 
pion degrees of freedom integrated out. The Lagrange density is composed of nucleon 
fields and has the form C = Ci + C2 + ■ ■ ■, where denotes the n-body terms. We 
have 



£1 = N^\idt + ^'^/{2M) 



N + ... 
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s 

Cs = (iVtiV)^ + . . . , 

A = Eo {N^Nf + ..., (7.2) 



where the elhpses denote higher derivative terms and Pj^^^ was defined in Eq. ( p.31|) 



for s =^5*0, ^5*1. Without derivatives, higher body contact interactions vanish because 
of Fermi statistics. Fermi statistics also imphes that there is only one four body term, 
which is invariant under Wigner symmetry. Furthermore, there is only one term in 
£3 which is also invariant^ To count operators it is useful to classify them by their 
transformation properties under SU(4). The three nucleon and anti-nucleon fields 
must be combined in an antisymmetric way, so the three A^'s (A^^'s) combine to a 4 
(4) of SU(4). Combining the 4 and 4 gives 1 © 15. However, the 15 does not contain 
a singlet under the spin and isospin SU(2) subgroups so there is only one three body 
operator with no derivatives. A group theory argument can also be used to show 
that there are only two two-body operators. Combining two A^'s in an antisymmetric 
manner gives a 6, while combining two A^"'''s gives a 6. Then 6x6 = 1© 15© 20, 
and only the 1 and 20 contain singlets under the spin and isospin subgroups. The 
Lagrange density £2 can also be written in a different operator basis: 



2 



C^iN^Nf + C^{N^aNf + . . . , (7.3) 



where C^'^"^ = Cl - ?>C'^ and C^'^'^ = Cl + . In this basis it is the term 
that breaks the SU(4) symmetry (as well as some of the higher derivative terms). 
The Lagrangian in Eq. |7]^ will be SU(4) symmetric at leading order if Cq = in the 
two-body sector and if higher derivative operators are suppressed in the three and 
four body sectors. 

In the MS scheme, neglecting two derivative operators the A^A^ scattering ampli- 
^P.F. Bedaque, H.W. Hammer, and U. van Kolck, private communication and Ref. p9^ 



tude is 



where 
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■^^'^ = T7^' (7-4) 



as was discussed section |2.3.2|. The scattering lengths are very large: a^^^°^ 



-23.714 ± 0.013 fm and = 5.425 ± 0.001 fm H, and even have the opposite 



sign. Therefore, the value of Cq^ and Cq^ "^^^ are very different. In this scheme the 
Lagrangian in Eq. ( |7.3| ) does not have Wigner symmetry. Nonetheless, for p ^ 1/ a^'^^ 
the amplitudes become, A^^'^ = Aixi/ [Mp). The equality of the ^5*0 and ^Si ampli- 
tudes is consistent with expectations based on Wigner symmetry. The p-dependence 
is consistent with expectations based on scale invarianceQ, since the cross section 
= 47r/p2. 



Recall from section 2.3.2 that in the MS scheme each term in the bubble sum in 
Fig. |2.1| gives a factor of ap, so for p > 1/a the series diverges. For these momenta 
minimal subtraction is not natural from a power counting perspective. In natural 
schemes like PDS and OS, each term in the bubble sum is of the same order as the 
sum. It is in these "natural" schemes that the fixed point structure of the theory and 
Wigner spin-isospin symmetry are manifest in the Lagrangian. In PDS or OS the 
coefficients depend on the renormalization point and for a^''-' ^ oo we have 

Cl:\m) = _\/^(.) - , (7.6) 

which is the same in both channels. In this limit 



^The scale transformations appropriate for the non-relativistic theory are x —> Ax, t X^t, and 
N ^ A-3/2iV. 
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and 

C, = --—{N^Nf + .... (7.8) 
Mur 

The first term in Eq. (|7.8|) is invariant under the Wigner spin-isospin transformations 
in Eq. ( [7. 1| ) . The ellipses in Eq. ( |7.8|) denote terms with derivatives which will not be 
invariant under Wigner symmetry even in the limit a*^*) oo. However, these terms 
are corrections to the leading order Lagrange density and their effects are suppressed 
by powers oi p/A (where A is a scale determined by the pion mass and Aqcd)- In the 
region l/a(^) < p < A Wi gner spin-isospin symmetry is a useful approximation and 
deviations from this symmetry are suppressed by 

C^ifiR) oc (l/a('^«) - l/a('^^)) , (7.9) 

/I C N 

and by powers of p/A. The measured effective ranges are Tq =2.73 ± 0.03 fm and 
= 1.749±0.008 fm H]. A rough estimate of the scale is 1/A ~ [rj'^''^ -rf,'^' V2 = 
0.49 fm, or A ~ 400 MeV. In PDS or OS, the limit a^'^^ ^ cxd is clearly a fixed point 
of Cq'^'' (/i/j) since nnd/dfiR [ij,iiC^^\jjii)] = 0. Also, scale invariance is manifest since 
f^R Z^r/A under scale transformations. 

In the KSW power counting potential pion exchange is order and is treated 
perturbatively. Thus, pion exchange is higher order than the iterated Cq bubbles and 
the theory still has Wigner symmetry at leading order. 

Wigner symmetry is useful in the two-body sector even though a^^^'^^ and a'^^'^^^ 
are very different. This is because for <^ p <^ A corrections to the symmetry 

limit go as (l/a^^'^''^ — rather than {a^^^°'' — a^^^'^^). This is similar to the 



heavy quark spin-flavor symmetry of QCD |gO|, 21], which occurs in the mq oo 
limit. Heavy quark symmetry is a useful approximation for charm and bottom quarks 
even though mfj/mc — 3. 

As an application of the symmetry consider A^A^ NN axion, which is relevant 
for astrophysical bounds on the axion coupling [|193|| . The axion is essentially massless. 
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Figure 7.1: Graphs contributing to NN — > NN axion at leading order. The sohd 
hues denote nucleons and the dashed hues are axions. 

If the axion has momentum k, and the initial nucleons have momenta p and —p then 
the final state nucleons have momenta q — k/2 and —q — k/2. Energy conservation 
implies that p^/M = q^/M + k'^/{AM) + k where p^\p\,q^\q\, and A; = |^|. In the 
kinematic region we consider q,p k, and the axion momentum can be neglected in 
comparison with the nucleon momenta. Again, this is just a multipole expansion (see 
chapter 6). In this limit the terms in the Lagrange density which couple the axion to 
nucleons take the form 

(7.10) 



where is the axion field and 5'o,5'i are the axion-nucleon isosinglet and isovector 
coupling constants. Associated with spin-isospin symmetry are the conserved charges 



(7.11) 



and the axion terms in the action are proportional to these charges 



Q- 



J3 



x=0 



(7.12) 
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The charge Q^^ is the total spin of the nucleons which is conserved even without taking 
the a^^^ — > oo hmit, however Q^^ is only conserved in the a^^^ — > oo limit (and also in 
the limit 

ai'So) _ a('Si))_ Since conserved charges are time independent, only a zero 
energy axion couples in Eq. (|7.12|) , and these terms will not contribute to the scat- 
tering amplitude. We conclude that NN(^Sq) NN{^Si)X^ vanishes in the limit 
a(') oo and that A^A^(^S'i) NN{^Si)X^ vanishes for all scattering length^. Cal- 
culation of the Feynman diagrams in Fig. shows that the leading order —>■ ^Si 
scattering amplitude does indeed vanish, and the NN(^Sq) NN{^Si)X^ amplitude 
is 



^ ^k^lJ- 1 ir 1 1 r 1 .7..^ 



where e'is the polarization of the final ^S*! A^A^ state. This is proportional to (l/a'-^'^"-' — 
l/a*^^"^^)) and is consistent with our expectations based on the Wigner symmetry. The 
fact that the graphs in Fig. |7[T| vanish as a^^^ 00 lends some support to statements 
in the literature ||194|| which claim that one-pion exchange is sufficiently accurate to 
describe the matrix element for A^A^ —>■ NN axion (at least for momenta ^ 1 /a such 
as in neutron stars). This process has contributions from different partial waves, 
and in all but the S-wave a single perturbative pion exchange is the leading order 
contribution in the KSW power counting. For the ^Sq — > '^5'i transition the graphs in 
Fig. [7.1| are small, so the first sizeable S-wave contribution occurs at NLO (the same 
order as the other partial waves) . It involves one-pion exchange and insertions of 
dressed by Cq bubbles. Besides tree level pion exchange these NLO contributions to 
A^A^ —>■ NN axion have not been considered in the axion literature. Note that to 
properly treat the S-wave contribution in neutron stars Pauli blocking effects would 
have to be incorporated in the calculation of the loop graphs. 

Coupling of photons to nucleons occurs by gauging the strong effective field theory 
and by adding terms involving the field strengths E and B. In the kinematic regime 

^ NN{^So) NN{^So)X° vanishes due to angular momentum conservation since the axion is 
emitted in a P-wave. 
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where the photon's momentum is small compared to the nucleons' momentum the 
part of the action involving the field strengths is 



Si^i = dtB^ 



KoQ'^ + K^Q^')+ ... , (7.14) 

x=0 ' 



where kq and ki are the isosinglet and isovector nucleon magnetic moments in nuclear 
magnetons, and the ellipses denote subdominant terms. The term proportional to ki 
in Eq. ( [7.14| ) gives the lowest order contribution to the amplitude for '-yd npi^So). 
The form of the coupling above implies that like the axion case, this amplitude is 
proportional to (l/a'-^'^"-' — l/a'^^'^^^). 

As our last example, we discuss the corrections to A^A^ scattering due to radiation 
pions discussed in chapter 6. As pointed out in chapter 6, one should perform a 
multipole expansion on the coupling of radiation pions to nucleons. The first term in 
the multipole expansion is: 



Q'' , (7.15) 

5=0 



where qa — 1-25 is the axial coupling and / ~ 131 MeV is the pion decay constant. 
Radiation pions also couple to a conserved charge of the Wigner symmetry in the 
large scattering length limit. (A multipole expansion is not performed on the cou- 
pling to potential pions so they do not couple to a conserved charge.) This implies 
that only a radiation pion with = will couple, which is incompatible with the 
condition ~ ^Jk"^ + m^, so in the symmetry limit radiation pions do not contribute 
to the scattering matrix element. In Eq. ( |6.31| ) we saw by explicit computation that 



graphs with one radiation pion and any number of Cq*'''s give a contribution that is 
suppressed by at least one power of l/a('^i)-l/a('^»). This suppression was the result 
of cancellations between the Feynman diagrams shown in Fig. |6.1| . Wigner symme- 
try guarantees that the leading contribution of graphs with an arbitrary number of 
radiation pions are suppressed by inverse powers of the scattering lengths. 

It has also been shown that Wigner symmetry is obtained in the large number 
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of colors limit of QCD ||195| , |196|| . The implications of Wigner symmetry in nuclear 
physics were studied in Ref. ^9^, |198|, |T99|, |00|, ^ So far the applications 

in this chapter have been specific to the two-nucleon sector, however Wigner symmetry 



m 



is observed in some nuclei with many nucleons. Recent progress |P04 , P05| , 177 , 178 
the three body sector suggests that the (A^^A^)'^ contact interaction is not subleading 
compared with the effects of the first two body term in Eq. (|7.8|) . (This conclusion 
is not uncontroversial, see Ref. ||206|| .) If the higher body operators with derivatives 
can be treated as perturbations, then the above discussion shows that approximate 
Wigner symmetry in nuclear physics is a consequence of the large A^A^ scattering 
lengths. 
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Chapter 8 Predictions for the ^Si Di 
Mixing Parameter, ei 



This chapter briefly discusses resuhs for the A^A^ — > A^A^ ^5*1 — ^ Di mixing param- 
eter, ei, at next-to-leading order (NLO) in the theory with pions. A number of 
observables have been computed at NLO in the KSW power counting. These include, 
nucleon-nucleon phase shifts [0, |, p07| , p.68| , |6[], coulomb corrections to proton-proton 
scattering |p08| , p09|| , proton-proton fusion ||210| , |2 1 1|| , electromagnetic form factors for 
the deuteron |[, deuteron polarizabilities ||160| , np ^ d'y [ pi2| , |2T3| , Compton 



deuteron scattering [|162|, |214|| , and parity violating effects [|163| , |161|] . Typically errors 
of order 30%-40% are found at LO and of order 10% at NLO. As mentioned in chapter 
5, this is consistent with an expansion parameter Q/A ~ 1/3. Since the expansion 
parameter is fairly large, calculations at next-to-next-to-leading order (NNLO) are 
useful. In the theory without pions, calculations at this order can be carried out in 
a straightforward manner p. With pions, the loop graphs become more difficult, 
but even two and three loop graphs can be evaluated analytically^. To truly test the 
convergence and range of the theory with pions, the above observables need to be 
calculated to NNLO. The phase shift in the ^5*0 channel has been calculated at this 
order ||168|| (and independently in Ref. ||164|| ), while work in the ^5*1, ^Di Si, and 



^Di channels is near completion | 215 |. In this chapter a presentation of the prediction 
for ei at order Q is givenQ. 

At order 1/Q there is no contribution to ei, which is consistent with the fact 
that this angle is much smaller than the ^Si phase shift. At order the graphs 
for ei include single potential pion exchange and pion exchange dressed on one side 



by Cq bubbles as shown in Fig. |8.1| M. This prediction does not involve any free 



^The basic reason that loops in three dimensions are simpler, is that the integrals can be done 
in position space where the Bessel functions redu ce to exponentials. 

^The material in this chapter appears in Ref. [ 216 with a more detailed discussion 
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Figure 8.1: The two order Q° diagrams that contribute to ei [Q. The sohd hues are 
nucleons and the dashed hues are potential pions. 

parameters. The order Q calculation includes a two derivative operator which causes 
transitions from to ^Di with coefficient C^^\ This operator is normalized so 
that on-shell in the center of mass frame the Feynman rule is 

^s,\^D, = -^ci'^^v\ (8.1) 

where p is the center of mass momentum. In this section the PDS renormalization 
scheme will be used. It was shown in chapter 5, that in the PDS scheme Cf^\^R) ~ 
so this operator enters at order Q. The Feynman diagrams that can contribute 
to ei at order Q include: 

i) one €2^^^ and any number of Cq '^^^'s 

a) one C2 one potential pion and any number of Cq "^^""'s 

in) one D2 one potential pion and any number of Cq ^^^^s (8-2) 

iv) two potential pions and any number of Cq s 

v) radiation pion corrections . 

The graphs for i) through iv) are shown in Fig. |8.2| . It can be shown that there are 
no radiation pion contributions to the mixing parameter at this order. 

The value of C2 D2 Cq and Cq^ are fixed from the ^Si phase shift 
calculation at order Q^. The only free parameter is Cf^\iXR), which is varied to give 
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Figure 8.2: Order Q diagrams for ei. The filled circle is defined in Fig. ^]T|, and the 
diamonds in c) denote insertions of the ^5*1 operators with C2 or D2 coefficients. 



a reasonable fit. The result of the ei calculation at order Q in the theory with pions is 
given by the dot-dashed line in Fig. |8.3| . The order result in the theory with pions 
^ is shown by the dotted line. The stars in Fig. |8.3| are data from Virginia Tech 
0. The open circles are the Nijmegen single energy fit to the data [§] whose quoted 
errors are invisible on the scale shown. The solid line is the Nijmegen multi-energy 
partial wave analysis 0. The value of C''2^\^r) corresponding to the fit in Fig. ^]3| 
is not given since we have not specified what constants are subtracted along with a 
p^/e pole that appears in one of the graphs with two potential pions. 

For comparison results have been shown for the theory without pions |]^. The 
long dashed line is the LO result and involves a one parameter fit. The smaller 
dashed line is the NLO result and involves fitting two free parameters. With fewer 
free parameters, the theory with pions does better than the theory without pions 
for p > 60MeV. In fact the theory without pions breaks down around mTr/2, which 
is where we expect it to since this is where the pion cut begins (see chapter 5). It 
has been noted in the literature | |217| that many observables do not test the power 
counting for perturbative pions. As can be seen from Fig. the mixing parameter 
provides an example in which perturbative pions clearly give improved agreement 
with the data. 
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Figure 8.3: Predictions for the ^Si Di mixing parameter ei in degrees. The sohd 
hne is the multi-energy Nijmegen partial wave analysis The long and short dashed 
lines are the LO and NLO predictions in the theory without pions ||^. The dotted 
line is the LO prediction in the theory with pions from Ref . . The dash-dotted line 
is the NLO prediction in the theory with pions. The stars are data from Virginia 
Tech [0 and the open circles are Nijmegen single energy data whose quoted errors 
are invisible on the scale shown. 



The dot-dashed line in Fig. |8.3| improves over the order Q*^ result for p < 140 MeV. 
For larger values of p, the order Q prediction grows, while the Nijmegen partial 



wave analysis has ei (p) < 3° for p < 300 MeV. In Ref. ||148|| the mixing angle was 
calculated using Weinberg's power counting and solving with a potential. In this 
approach graphs with potential pions are summed up. A direct comparison with this 
calculation is difficult since graphs with A's were included, and more parameters were 
varied in the fit. When the potential method is used there is cutoff dependence in 
the result which cannot be cancelled by cutoff dependence in the coefficients. This 
is because contributions that are formally higher order are included in the answer. 
Since the results are numerical these can not be thrown away. For a cutoff of order 



nip, the prediction in Ref. [148] also grows with p. In the potential approach the 
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cutoff dependence gives a measure of the uncertainty due to higher order corrections. 
Varying the cutoff from 0.6 mp to 1.3 mp the prediction varies by ~ 1.0° at p = 
150 MeV. Therefore, at 150 MeV the theory with perturbative pions seems to be 
doing no worse than a calculation where the pions are treated non-perturbatively. 
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Chapter 9 Conclusion 

In this thesis effective field theory techniques are used to describe the interaction of 
heavy particles at low momentum in a model independent way. The focus has been 
on interactions with pions, described using chiral perturbation theory techniques. 
Chapter 2 gives an overview of the formalism for describing interactions of one or two 
heavy particles. 

In chapter 3, the decays D* — > Dn and D* — > D7 are studied using heavy meson 
chiral perturbation theory. With the recent measurement of B{D*'^ D'^l)-, the 
and Dl branching fractions can be used to extract the D*D7r and D*D^ 
couplings g and The D* —>■ Dj decays receive important corrections at order y^m^ 
and, from the heavy quark magnetic moment, at order l/rric- Here all the decay 
rates are computed to one- loop, to first order in nig and l/rric, including the effect of 
heavy meson mass splittings, and the counterterms at order m^. A fit to the exper- 
imental data gives two possible solutions, g = 0.27 ^02- 02) P — 0.85^'^ ^"J GeV^"*^ or 
g = 0.761: fsti, 13 = 4.901 J 1^7° GeV-^ The first errors are experimental, while the 
second are estimates of the uncertainty induced by the counterterms. (The experi- 
mental hmit Tr)*+ < 0.13 MeV excludes the g — 0.76 solution.) Predictions for the 
D* and B* widths are given. 

In chapter 4, the prospects for determining \Vu}\ from exclusive B semileptonic 
decay are discussed. The double ratio of form factors 

is calculated using chiral perturbation theory. Its deviation from unity due to contri- 
butions that are non-analytic in the symmetry breaking parameters is very small. It is 
concluded that combining experimental data obtainable from B ^ p^i^e, B ^ K*£i 
and D pii^e can lead to a model independent determination of \Vub\ with an 



147 

uncertainty from theory of about 10%. 

In chapter 5, an effective field theory description of nucleon nucleon interactions 
is investigated. A momentum subtraction scheme (OS) is introduced which obeys the 
power counting of Kaplan, Savage, and Wise (KSW). The KSW power counting was 
developed for systems with large scattering lengths, a. Unlike the power divergence 
subtraction scheme (PDS), coupling constants in this scheme obey the KSW scaling 
for all > 1/a. This chapter explains in detail how the renormalization in the OS 
and PDS schemes is implemented using local counterterms. The main complication 
is the need to include an infinite number of counterterms since the leading order 
result includes an infinite number of loop graphs. Fits to the NN scattering data are 
performed in the ^5*0 and ^Si channels. An error analysis indicates that the range of 
the theory with perturbative pions is consistent with 500 MeV, so it can be concluded 
that there is no obstruction to using perturbative pions for momenta p > m^. Some 
comments are made on the low-energy theorems derived by Cohen and Hansen^]. 

In chapter 6, radiative pion interactions are investigated. For interactions in- 
volving two or more nucleons it is useful to divide pions into three classes: potential, 
radiation, and soft. The momentum threshold for the production of radiation pions is 
Qr = \^MN^m^. It is shown that radiation pions can be included systematically with 
a power counting in Qr- The leading order radiation pion graphs which contribute 
to NN scattering are evaluated using the the PDS and OS renormalization schemes 
and are found to give a small contribution which vanishes as the singlet and triplet 
scattering lengths go to infinity. The power counting for soft pion contributions is 
also discussed. 

Chapter 7 shows that in the limit where the NN ^Sq and ^Si scattering lengths, 
a*^^"^") and a^'^^'^\ go to infinity, the leading terms in the effective field theory for 
strong A^A^ interactions are invariant under Wigner's SU(4) spin-isospin symmetry. 
This explains why the leading effects of radiation pions on the S-wave A^A^ scattering 
amplitudes vanish as a^^"^"^ and a^^^'^^ go to infinity. The implications of Wigner 
symmetry for A^A^ NN axion and '-fd—^np are also considered. 

Finally, in chapter 8 results for the ^Si —^Di mixing parameter, ei are presented. 
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This observable provides an example where the theory with perturbative pions gives 
better agreement with the data at p ~ m^r with fewer parameters than the theory 
without pions. 
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